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WALTER E. LILJESTRAND 


Chief Research Engineer, 
Mission Manufacturing Company, 
Houston, Texas 


Simplified Approach to Understanding the 
Suction Problems for Reciprocating Mud Pumps 


There are no standards or generally approved formulas for designing or analyzing mud- 


pump suction systems such as already worked out for centrifugal pumps. 


This paper 


analyzes pump-suction and suction-piping characteristics showing responsibility areas 


of users and manufacturers. 


plained. 


Empirical formulas and verifying test procedures are ex- 


This material should be used as a basis for preparation of standards on pump 


suction systems as used in the oil field. It will reduce pumping troubles and costs and 


| design of mud-pump suction systems has ap- 
parently been done empirically or by experience. The results 
have been reasonably satisfactory, but many instances of trouble 
have raised questions about the corrections needed. Corrective 
steps often did not help because the true cause of the trouble was 
not known. There are a few people who know pump-suction re- 
quirements and understand the hydraulic effects of the suction- 
piping system, but their knowledge has not been made broadly 
available. Standards have not been prepared and distributed. 
Increased competition in drilling requires that pumping costs be 
kept as low as possible and that all pumps perform at maximum 
output and efficiency. 

A good understanding of the mud-pump suction problem ap- 
plied to suction systems can give the following advantages: 


1 Less replacement of fluid-end parts. 

2 Longer life of power-end drive system. 

3 Increased pump capacity. 

4 More dependability of pump under adverse conditions. 
5 Possibly less costly suction-piping system. 


~ 


Less concern and worry by the crew. 
7 Effective use of suction dampeners and prechargers. 


P Contributed by the Petroleum Division and presented at the 

Petroleum Engineering Conference, Houston, Texas, September 20- 

23, 1959, of THE AMERICAN Society OF MECHANICAL ENGINEERS. 
Note: Statements and opinions advanced in papers are to be 


understood as individual expressions of their authors and not those 
Manuscript received at ASME Headquarters, July 
Paper No. 59 


of the Society. 


10, 1959. Pet-23, 


Nomenclature= 


permit increased pump output. 


Centrifugal-Pump Approach 


The chemical industry has taken advantage of studies and test- 
ing of suction systems of centrifugal pumps. Manufacturers 
have supplied accurate data on centrifugal-pump suction charac- 
teristics using water. Users have studied and designed suction 
systems so that good pump performance is assured even when the 
suction systems are built at minimum height and therefore at 
lowest cost. Cavitation problems are in the minority and, in case 
one occurs, simple direct corrections are easily made. The flow 
through the suction system of a centrifugal pump is steady for a 
normal condition so the conditions are easily measured. The re- 
quirements for a different installation are easily specified. Per- 
formance in the field is the duplicate of test-stand performance. 


Reciprocating-Pump Approach 

The flow is cyclic in the suction-piping system of a reciprocating 
mud pump because it follows the start and stop motion of the 
pump pistons. Each pump revolution the flow pattern is re- 
peated. The pattern depends upon the number of cylinders, 
single or double-acting, piston-rod area, connecting-rod crank- 
length ratio, and relative angular positions of cranks. Flow rate 
at any instant in the pattern or revolution also depends upon pis- 
ton area, piston-stroke length, and pump-crankshaft speed. De- 
tailed explanations of the flow pattern can be found in other 
writings. 

The cyclic flow pattern makes the analysis of the conditions in 
a suction system more complex than for centrifugal pumps. The 
effect of the cyclic flow must be understood in order to design a 


NPSHieq’'a = net positive suction head 


required, psia 


flow velocity in suction- 
piping system, psia 


taneous flow divided by 
pipe flow area, fps 


NPSHavair = net positive suction head Py. = deceleration head  pro- V?/2g = velocity head, ft of mud 
available, psia duced during decrease P. = reading of a suction gage 
P, = atmospheric pressure, psia of mud-flow velocity in in psig with P, as datum 
Z = elevation of free surface of suction-piping system, A, = acceleration of pump pis- 
mud above datum eleva- psia ton, fps per sec 
tion, ft (or —Z if below P., = vapor pressure, psia N = speed of pump crankshaft, 
datum) sp gr = specific gravity of fluid or rpm 
h, = head lost at entrance of mud weight, lb per gal/ S = stroke of pump piston, in. 
suction-piping system, 8.33 L = length of suction pipe, ft 
ft of mud g = acceleration due to gravity G = mud weight, lb per gal 
h, = head lost in friction in 32.17 fps per sec D, = piston diameter, in. 
suction-piping system V = velocity in suction aver- D, = suction pipe inside diame- 
from entrance to pump aged for cross section of ter, in. 
suction flange, ft of mud pipe flow area deter- R = length of pump crank, in. 
Pye = acceleration head used mined at point of gage y = elevation of suction gage 


during increase of mud- 
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above or below datum, ft 
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suction systeni. In fact, it will be shown later in this paper that 
in the suction-piping systems for conventional mud-pump in- 
stallations (not reciprocating pumps in general) the physical 
dimensions are such that the cyclic-flow-pattern effect is dominant 
and is a factor indirectly limiting the capacity of the pump. It 
controls to such a degree that understanding and control of it will 
solve most of the problems for mud pumps. In this respect, it is 
reasonable to say that the understanding of suction problems for 
mud pumps is simplified by approaching it through the study of 
the cyclic flow. 


The Suction System 


The complete suction system of a mud-pump installation, Fig. 
1, will be referred to as consisting of two parts divided at the 
pump-suction flange. The part from the suction flange to the 
faces of the pistons and to the discharge valve decks will be called 
the ‘pump-suction manifold.”’ The part from the suction flange 
to the free mud surface will be called the “suction-piping sys- 
tem.’”’ The object of the complete suction system is to fill the 
pump cylinders with fluid. 


| 
| 


Pump Suction | 


~— —Pt—Suction Piping System 
ManiFoLo | 


Fig. 1 Physical characteristics of a pumping suction system 


Pump-Suction Manifold 

Failure to fill the pump at the same rate that the piston strokes 
will cause “fluid knocking” which is detrimental to pump parts 
and limits the pump capacity. Fluid knocking is a loud sound 
in the pump caused when the fluid catches up with the piston 
after lagging behind at the beginning of the suction stroke. Its 
Collier! 
A characteristic of the pump-suction manifold is that it requires a 


characteristics and causes were explained by S. L. 


small absolute total pressure or head to: 
1 Provide the velocity of flow in the pump-suction manifold. 
2 Overcome friction loss in the pump-suction manifold. 
3 Operate the suction valves. 
4 Provide when needed the acceleration of mud in the pump- 
suction manifold. 
5 Overcome vapor pressure of the mud. 


Friction and velocity head increase with the flow rate. Ac- 
celeration of mud increases with the pump speed and piston 
diameter. Acceleration of mud or pressure needed for accelera- 
tion head is greatest for the leading head-end piston just starting 
its suction stroke and it will have four high values each pump 

18. L. Collier, “Know Your Mud Pump—It Can Save You Money, 
Part 1—Introduction,’’ World Oil, vol. 147, July, 1958, p. 100. 

2S. L. Collier, ‘‘Know Your Mud Pump—It Can Save You 
Money, Part 2—Suction Pipe and Acceleration,’’ World Oil, vol. 147, 
August, 1958, p. 81. 
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revolution for a duplex double-acting pump. The absolute pres- 
sure value of the fluid at the suction flange (corrected to datum) 
that is just sufficient to overcome in the pump-suction manifold 
the friction loss, valve-opening loss, velocity-head increase, and 
the highest acceleration head in the cycle after subtracting the 
absolute vapor pressure is called the ‘‘net positive suction head”’ 
required or NPSH req’'a- The pressures needed for the other three 
high values are smaller because of effect of piston-rod area, con- 
necting-rod effect, and flow pattern in the manifold, so are of little 
interest because the highest value is the first one to cause fluid 
knocking. To prevent fluid knocking, or to insure filling the 
pump at the same rate the piston strokes, the suction-piping 
system external to the pump-suction flange must always supply 
fluid at a pressure equal to or greater than the NPSH req’a. 


Suction-Piping System 

The total pressure of mud (absolute pressure corrected to 
datum plane plus the velocity head) at the pump-suction flange 
delivered by the suction-piping system less the absolute vapor 
pressure of the mud is the “net positive suction head” available 
or NPSHavait. It consists of head values such as atmospheric 
pressure, the height of free surface above or below datum plane 
and the specific-gravity effect less losses such as entrance loss, 
friction loss, acceleration head, and less the absolute vapor pres- 
sure. In symbols this is expressed as: 
sp gr sp gr 


— — (h 


NPSH 


The fluid-head values in Equation (1) are all a part of the suction- 
piping system. They are dependent on the: 


1 Rate of mud flow. 
2 Rate of change of mud flow (acceleration). 
3 Physical size and condition of the suction-piping system. 


They can be calculated or measured. Considering a conven- 
tional-sized suction-piping system for the average mud pump, the 
scale values of pressures and losses shown in Fig. 2 are approxi- 
mately correct. The atmospheric pressure P, and the static head 
Z are easily determined so the only other major portions to de- 
termine are acceleration head Py. and vapor pressure Py». In 
this respect the suction problem for mud pumps is greatly simpli- 
fied because acceleration head is dependent on the physical size of 
the suction-piping system. The acceleration head Py. can be 
calculated. One expression for this was used by 8. L. Collier and 
it is shown in the Appendix. 


Cyclic-Flow Aspects 


There is a simultaneous need for acceleration of mud in the 
pump-suction manifold and the suction-piping system during the 
start of a given piston stroke. In a duplex double-acting pump 
there are four accelerations each revolution. 
the mud-flow velocity diagram, Fig. 3. 


These are shown in 


The effect of the acceleration in the pump-suction manifold is 
to increase the N PSH req’a because this acceleration is an included 
factor. In contrast, the effect of the acceleration in the suction 
piping system P,- is to decrease the NPSHavaii because this ac- 
celeration is an excluded factor. 

The value of NPSHavaii can be determined for any point on the 
revolution shown in Fig. 3. At point A, the highest velocity, 
the entrance loss A, and friction loss h; are maximum but still 
very small and the acceleration head P,, is zero. The NPSHavait 
will therefore be a large value at point A. At point B on the other 
hand the entrance loss h, and the friction loss h, are somewhat 
reduced and still very small, but the acceleration head Pye is 


Ibid. 
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large and at its maximum value. The NPSHavaii will be minimum 
at B for the revolution. Fig. 2 is an illustration of conditions at 
point B in Fig. 3. Here the slope of the velocity curve shows 
maximum mud acceleration. 

The scale values shown in Figs. 2 and 3 are dependent upon the 
speed of the pump. Fig. 4 illustrates the character of NPSHavait 
and N PSH req’a Versus speed of pump. 

The NPSHavaii at zero speed will be 


sp gr 
2.31 


NPSHawvail = P, = (2) 


The NPSHavaii is represented in Fig. 4.as a band or region. Ata 
given pump speed the width of the band is produced by alternate 
accelerations and decelerations of the mud. The deceleration 
just preceding point B of Fig. 3 represented by the expression Pde 
. is found experimentally to be equal or greater than the accelera- 
tion Px. following point B. The pressures caused by deceleration 
are of no direct value for fluid-knock prevention because they do 


hy 


Fac 


$1 
=| 
excess 


NPSH- avail. 
NPSH-reqd. 


Pie 


Fig. 2. Pressure diagram of pressures and heads (all in psi) showing 
scale values for suction system shown in Fig. 1 
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Fig. 3 Typical suction flow velocity diagram for a duplex double-acting 
pump 
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PUMP SPEED 


Fig. 4 Ranges of NPSH values with changing speed. The width of 
NPSH region at a given speed represents the maximum cyclic effect of 
acceleration (and deceleration). 


not last long and are gone before the following acceleration takes 
place. The term Pa. was purposely omitted from Equation (1). 
The minimum values of NPSHavaii obtained at time of maximum 
acceleration are of chief concern. 

The NPSHreq’a by the pump at zero speed is the pressure re- 
quired to open a suction valve. The band width of the NPSHreq’a 
for pump speeds is also caused by the cycle of the pump. The 
maximum value at a speed occurs the same instant in the revolu- 
tion that the minimum NPSHavair occurs. The difference be- 
tween them or ‘“‘excess’’ of NPSHavair over NPSHreq’a is the re- 
serve available to take care of unexpected adverse changes. Such 
changes may be increased vapor pressure, partially clogged line, 
loss of head of mud, and so forth. The maximum speed M that 
the pump can run without fluid knocking is found at the inter- 
section of the NPSHavaii minimum line with the NPSHreq’a 
maximum line. These two lines of course represent the condition 
for the leading head-end cylinder. The acceleration-head peaks 
for the other cylinders were smaller so the pump speed at which 
fluid knocking begins in them will be a little higher. This is one 
of the reasons that fluid knocking is very irregular when it is first 
encountered with increasing pump speed. Plunger pumps with 
symmetrical suction-velocity diagrams will also show fluid-knock- 
ing irregularities with speed increases. 

The position of the maximum NPSHreq’a curve in Fig. 4 for 
a given pump and piston diameter is rather fixed. Its position 
and shape are not dependent on suction-piping-system conditions. 
The position and shape of the minimum NPSHavaii curve how- 
ever are dependent on dimensions and suction conditions of the 
suction-piping system. This provides the possibility of collecting 
field test data over a period of time to plot the NPSHreqa for the 
pump. Each time a new operating condition occurs such as a 
change of static head, mud weight, mud temperature, or pipe 
dimensions, a test should be run by raising the pump speed 
until fluid knocking develops. This is the speed at which the two 
curves cross in Fig. 4. Values of NPSHavaii can be computed if 
the suction-piping system is normal, that is, not clogged or leaking 
air. Each change in conditions moves the NPSHavyaii minimum 
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curve and provides a new crossing point. A series of such points 
determines the NPSHreq’a curve. 

Each computed or measured value used in Equation (1) to 
obtain NPSHayaii and to determine the NPSHreq’a will have a 
tolerance or probable error. These unknowns should not be over- 
looked but it is suggested that the allowance made on each deter- 
mined value be as small as possible so that the excess of NPSHavaii 
over NPSHreq’a be as large as possible. This will give a clearer 
picture of how much allowance is available for unknowns and the 
unexpected. 


Determination of NPSH...,-. 

There is no way to calculate the NPSHreq’a for centrifugal 
pumps, so their manufacturers have worked out very accurate 
testing methods. It may be possible to write equations for 
NPSHereq’a for piston pumps after sufficient experimental back- 
ground is obtained. Manufacturers of piston pumps can also 
measure NPSHreq’a quite accurately. Fig. 4 shows that 
NPSHereq’a is equal to NPSHayair at the speed where the pump 
just begins to knock. NPSHayaii can be very accurately deter- 
mined at this speed on a special test stand Fig. 5, designed so that 
all the difficult-to-measure terms of Equation (1) are zero. The 
terms of h,,h;, and P.- will be zero for the suction-piping system if 
the pipe length is zero, Fig- 6. In other words the mud tank 
should be connected to the pump with a short well-rounded noz- 
zle. The tank must be closed and capable of sustaining full 
vacuum. The surface of the mud should be only a few feet above 
the datum plane of the pump and kept constant by returning the 
pump discharge to it. Also the mud should be deaerated. Con- 
tro] of the air pressure over the mud will give successive test points 
for various pump speeds as shown in Fig. 7: 


At knocking 
NPSHreq’ta = NPSHavaii 


NPSH,..11 = P, — Px — Pop —— 2 3) 


Such data should be obtained on each pump style and size for 
various speeds and piston sizes. Water should be the first fluid 
tested but it should be followed by muds of various types and 
weights. 

Electronic suction-gage readings will faithfully reproduce the 
eyeles but will be difficuit to interpret. It may eventually be 
possible to use these data to determine NPSH values. The 
problems of interpretation of these data will be partially evident 
by referring to the Appendix where interpretation of a gage read- 
ing is shown with several actual oscillograph readings of suction 
pressure. 


VACUUM 
= 


required by a pump 


90 MAY 1960 


Fig. 5 Schematic diagram of laboratory test stand for measuring NPSH 


Work to Be Done 


Thus far in this paper the general nature of the suction problem 
and relation of its components has been presented. If the user 
wishes to obtain the benefits offered by having a correctly designed 
system, further work on this approach must be done to refine it 
and to agree on details of symbols, terms, computations, tables, 
curves, and so on. The available papers relating to suction 
problems show considérable variance in symbols, terms, and 
definitions. Early agreement must be reached on these by the in- 
terested parties. 

More data for correct evaluation of terms must be obtained, 
For example, the vapor pressure of various muds should be 
checked to see if water vapor pressures are sufficiently correct as 


Re =O 


Poi 


NPSH-avail. 


Fig.6 Pressure diagram for suction system using the test stand of Fig. 5. 
If the vacuum P., is increased at the test speed, the NPSH,,.\; will be 
reduced to equal NPSH,,...;. Fluid knocking develops when this test 
condition is reached. 
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Fig. 7. NPSH,....4 curve determined by successive laboratory pump suc- 
tion tests using test stand of Fig. 5 


Transactions of the ASME 


he = O 
Psg lan 
he 
Z 
ve 
: 
t 
= 
at 
> Test 
| 
ol : | 
| | 
slid Ww = 
| 


substitutes. This may be quite important for the higher mud 
temperatures and different mud types and weights. If water 
vapor pressures give conservative results by being substantially 
higher than mud vapor pressures, the suction-piping systems 
mav be expensive being overdesigned. Acceleration head of the 
suction-piping system being a dominant term, agreement should be 
reached on a correct expression for computing it. If possible, the 
results should be shown in curves or tables to simplify using it in 
design work. It might even be possible within a given user com- 
pany to work out standard suction-piping systems. 

NPSH ean be expressed in either feet of water (as in centrifugal 
pumps) or in psia. Psia may be the best for oil-field work con- 
sidering that mud is a normal fluid and is not of standard weight. 

No emphasis was placed in this paper on the importance of a 
standard reference plane or datum plane for static head Z and for 
correcting pressure-gage Curves and charts 
could be quite confusing without agreement on the location of this 
datum plane in the pump. Locations might be the piston-rod 
center line, the top of the discharge valve deck, or the center line 
of the suction flange. The height of the plane above the feet of 
the pump is an indirect factor affecting the static head Z. 


measurements. 


Conclusion 


Several steps can be taken to obtain full benefit of suction-pip- 
ing svstems designed to eliminate or control fluid knocking. The 
first step can be taken immediately. 


1 Use the formula for acceleration head and design systems 
for minimum acceleration. In general the shortest, straightest, 
and largest suction pipe will produce the largest NPSHayaii. The 
excess over NPSH,.<’s will not be known but can be large 
enough to take care of most unexpected changes. 

2 Determine the NPSHqa on pumps by field testing at 
knocking speeds. This burdens the pump with all the short cuts 
of this approach but it will then be possible to match a suction 
system to a pump. 

3 Determine the NPSHreq’a by laboratory tests to give pump 
all the credit it is due. 

Refine equations for caleulating NPSHavait in light of ex- 
perience with laboratory tests to make field-test data agree. 
Standardize the svmbols and methods. 
data, 


Publish curves and 
The last step may not be completed for several vears. Begin- 
ning the first step will reduce system costs, operating costs, and 
increase the pump capacity. 

Some arrangements of piping systems may have to have 
suction-pulsation dampeners to reduce the acceleration head and 
increase the NPSHayaii. In some instances where suction pipes 
must be long or lift high or gas is entrained in the mud, a pre- 
charging centrifugal pump when properly selected and installed 
will produce very beneficial results. Poor results on precharging 
in the past can be attributed to a lack of understanding of the 
problem. 


Increased speeds of pumps, especially plunger pumps, make it 
imperative that suction problems be completely defined. When 
completely defined, such knowledge would be of great value to 
related reciprocating-pump problems such as salt-water disposal 
systems and well-servicing pumps. 

Use of suction characteristics can be worked into in easy 
stages, but the most progressive and cost conscious will go to work 
on it at once. Until more experience develops, there should be 
left a definite excess of NPSHavair over NPSHyeq'a to take 
care of tolerances and unexpected changes. 
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APPENDIX 

Suction Lifts 

In many cases suction pits are used instead of tanks. Equation 
(1) expresses these conditions using the minus sign for the static 
lift Z. Fig. 8 shows scale values of pressures and losses. Even if 
all other factors remain the same, the loss in static head with pits 
compared to tanks will reduce NPSHavaii substantially. The 
excess over NPSH req’a will be small and may even require 
running the pump slower. The entrance head loss h, may in- 
crease unexpectedly by clogging of the inlet. 


PSI 


excess 
yt 

29 
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Fig. 8 Pressure diag for system with a suction lift from pit 
ortanks. The reduction of NPSH,,.\) is substantial compared to Fig. 2. 


Acceleration Head — P.. 


A theoretical equation (5) for acceleration head P,. can be 
derived considering: (a) The acceleration of the piston, and (6) 
the force required to accelerate the mud in the suction-piping 
system. 

General form equation for acceleration of 

# a body (crankpin) moving in a circular 

? path of radius r with speed v. Units: 
v in fps, and r in ft, a in fps per sec. 


2 Substituting terms used in a suction sys- 
TN? 
~ 10,800 
~ oe Change radius R to stroke S/2. 


= 0.000456.V2S Apply connecting-rod-effect factor = 
1.20 for connecting-rod length = 5 crank 
lengths. (4) 
General form equation for force in Ib 
which accelerate mass (weight divided by 
acceleration of gravity). 


A,, = 0.000548.N2S 


P 
F = ma 
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Fig. 9 Pressure diagram for suction system (Fig. 1) showing relation- 
ship of a vacuum gage reading (Fig. 11) to other values. Its cyclic 
changes are similar to those of acceleration (and deceleration) head. 
Readings from a mechanical gage are difficult to make. 


VAGUUM 


PRESSURE 
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Fig. 10 Relationship of suction gage reading (used in Fig. 11) to other 
elements (converted to psi) entering into an NPSH,,,i; measurement. 
The velocity head is a portion of the suction gage reading P., when the 
gage is reading a vacuum. 


Considering the volume and specific gravity of the mud in the 
suction-piping system moving with acceleration proportional to 
ratio of piston area and suction-pipe inside area: 


Py X = — 4 
231x3217.. + ~ 
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Fig. 11 Installation of suction gage on suction-piping system near pump 
manifold suction flange 


Fig. 12 Suction manifold pressures for different speeds and liner sizes. 
One revolution of pump lapsed for each trace shown. Pressure scale is 
170 psi perin. Atmospheric line is shown. The effect of four piston 

lerati lution with r is seen in these complex 


per rev 


substituting terms used in a suction system 


D 2 
Pye = 0.00161 X LX GX (22) (5) 


These equations show that fluid knocking if caused by too 
much acceleration head P.- loss can be reduced by: (a) Reducing 
pump speed N, (b) reducing suction piping-system length L, (c) 
reducing mud weight G, (d) reducing piston diameter D, or in- 
creasing suction-pipe diameter D,. 


The Suction-Gage Reading 


Much misunderstanding of suction-gage readings is possible. 
The location of the gage must be considered in order to be able 
to select the factors making up its reading. Figs. 9 and 10 show 
the pressures and losses to scale in bar form and on the gage face 
for an installation shown in Fig. 11. The values are instantaneous 
for maximum acceleration head Pa. The standard vacuum gage 
is not recommended because its hand will swing widely from 2 to 
20 eps. 

Oscillograph records of suction pressures are shown in Fig. 12. 
These are samples taken from hundreds and were selected to 
show the variety possible. Many conclusions can be drawn by 
intensive study of all the suction records. The important one is 
that the minimum suction pressure occurs when the need for 
pressure is greatest. 

The four accelerations per revolution (of a duplex double-acting 
pump) may appear in simple form or in multiples up to 16 per 
revolution. Fluid knocking superimposes and may obliterate 
some of the pattern causing an odd number of surges per revolu- 
tion to occur. 
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gested. 


cracking, which is defined as the 
spontaneous fracturing of steel subjected simultaneously to a 
corrosive hydrogen-sulfide aqueous medium and a static stress 
less than the ultimate strength of the metal, became of major im- 
portance to the oil industry with the discovery of sour-condensate 
production in the late 1940’s. Sour-condensate reservoirs, which 
have been discovered in Alberta, Canada; Arkansas, Wyoming, 
Mississippi, and Texas; and in Lacq, France; are produced from 
limestone formations that are located at depths ranging from 
10,000 ft to more than 15,000 ft. Bottom hole pressures and 
temperatures are in the range of 5000 to 9500 psi and 180 to 280 F, 
respectively. The hydrogen-sulfide content of the well fluids 
varies from 1!/2 per cent in the Fashing Field, Tex., to 46 per cent 
at Neiber Dome, Wyo. Carbon dioxide is always present in 
lesser and unrelated amounts ranging from 2 per cent (20 per 
cent H.S) at Lacq, France, to 9 per cent (19 per cent H,S) at 
Federal-Stockham, Wyo. 

In the early development of the fields at Pincher Creek and 
Jumping Pound, Alta., Canada; Federal-Stockham, Wyo.; and 
Lacq, France; serious difficulties from sulfide-stress cracking 
failures in critical items of well production equipment were ex- 
perienced. Because of the hazards and the great expense in- 
curred by such failures, the sulfide-stress cracking problem was 
studied extensively by each of the operating companies, and on a 
co-operative basis under the auspices of the National Association 
of Corrosion Engineers. The test methods that have been used by 
the various researchers and the results that were developed have 
been published [1-10]! and will not be discussed in detail. How- 
ever, a brief review of the factors that are now believed to bear 
importantly on the sulfide-stress cracking phenomenon will be 
included for ready reference. 

Factors Influencing Sulfide-Stress Cracking. Despite the extensive 
investigatory work that has been done, a completely satisfying 
explanation of the sulfide-stress cracking mechanism has not been 
evolved. The various theories that have been advanced to ex- 
plain the failure mechanism involve: (a) The effect of nascent 
hydrogen absorbed by the steel; (b) true stress-corrosion crack- 
ing; and (c) a combination of both. In instances where the 
cracking occurs very rapidly, in a matter of hours or a few days, 
absorbed hydrogen is believed to be responsible for the failure. 
When the time to failure is prolonged to months or years, true 
stress-corrosion cracking appears to predominate. Failure at 
intermediate time intervals may be a consequence of the com- 


1 Numbers in brackets designate References at end of paper. 
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Prevention of Sulfide-Stress Cracking in 
Sour-Condensate Well Equipment 


Several case histories of equipment failures in sour-condensate wells are presented, 
reasons for these failures are discussed, and means for avoiding such failures are sug- 


bined effects of absorbed hydrogen and stress-corrosion cracking. 

The problem of hydrogen in steel is not new and it has been 
known for many years that hydrogen embrittlement, which is 
manifested by loss in ductility without a corresponding loss in 
yield strength or ultimate strength, will develop in steels after 
exposure to a corroding hydrogen-sulfide environment. While it 
is conceded that hydrogen embrittlement will precipitate 
failure under circumstances where dynamic loading is involved, it 
is not considered a basic cause of the failures that have developed 
in sour-condensate service. The mere fact that the material has 
been embrittled does not mean that it will fail spontaneously at 
applied stresses that are within the elastic range. 

The current theories [5, 8, 11, 12, 13] on the manner in which 
absorbed hydrogen influences sulfide-stress cracking are based on 
solid state reactions between absorbed nascent hydrogen and 
highly stressed steel. Each of the theories incorporates the con- 
cept that the interaction of the absorbed hydrogen with disloca- 
tions and/or lattice imperfections is responsible for a local tensile 
stress, which when added to the applied and internal stresses is 
sufficient to initiate cracking. Inthe authors’ opinion, stress-cor- 
rosion plays a minor role, since virtually all of the sulfide-stress 
cracking failures have developed rapidly and apparently cannot 
be stopped by cathodic protection. Corrosion appears to be of 
importance only in so far as it is the source of nascent hydrogen. 
This observation has been supported by the results of cathodic 
charging experiments [8] in which elastically stressed steel speci- 
mens failed in essentially the same manner as experienced under 
sulfide conditions. 

Despite the disagreement on the details of the failure mechan- 
ism, it is generally accepted that the cracking tendency of ferritic 
steels in sour-condensate service is related to: (a) Absorption of a 
sufficient quantity of nascent hydrogen by the steel; (b) a total 
tensile stress of critical magnitude; (c) a susceptible metallurgical 
condition; and (d) time. If any one of these factors is absent or 
is of insufficient magnitude, sulfide-stress cracking will not de- 
velop. Absorption of nascent hydrogen is a factor that cannot be 
eliminated, since, in most instances, it is impractical to shield 
major items of equipment from the sour well fluids. The fact that 
the failures that have been experienced have developed after short 
exposures to the sour well fluids precludes replacement of equip- 
ment after a finite life as a means of preventing failure. 

Susceptibility of ferritic steels has been ascribed, at least in 
part, to the presence of martensite and some of the structures re- 
sulting from tempering martensite. Since the mechanical proper- 
ties are related to the microstructures that are developed, they 
have been used as a measure of susceptibility. For convenience, 
hardness, which at best is a rough parameter of the mechanical 
properties, has been used as a criterion of susceptibility. Steels 
with hardnesses of Rockwell-C 20 or lower are highly insensitive 
to sulfide-stress cracking, provided they have not been subjected 
to an excessively high degree of plastic deformation. In this 
hardness range, the structures are usually mixtures of pearlite and 
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ferrite, or spheroidized carbide and ferrite. At hardnesses of 
Rockwell-C 30 and higher, all ferritic steels are susceptible and 
become increasingly sensitive to cracking as hardness increases 
even in the absence of plastic deformation. In the range Rock- 
well-C 20 to 30, the cracking tendency is highly variable and is 
strongly influenced by the amount of plastic deformation and, 
more importantly, by the type of microstructure that has been 
developed. For a given hardness in this range, steels may ex- 
hibit the same or decidedly different structures 
variable mechanical properties 


and therefore 
depending upon the chemical 
composition and the manner of processing. In some steels mixed 
structures containing appreciable quantities of white martensite 
are developed by normalizing; in others the structures may be 
uniform  pearlite-ferrite mixtures. In properly quenched 
and tempered materials, a uniform structure commonly called 
tempered martensite is developed. The presence of damaging 
martensitic structures can be avoided by appropriate heat-treat- 
ment, which for alloy steels will involve tempering after cooling 
from the austenitizing temperatures. 

The critical stress necessary for initiation of sulfide-stress crack- 
ing is the total tensile stress to which the metal is subjected and in- 
cludes the internal stress present in the metal as well as the ap- 
plied service stresses. Although the magnitude of the applied 
stress can be precisely determined, it is virtually impossible to 
determine with any degree of accuracy the magnitude of the in- 
ternal stresses that are induced as a consequence of heat-treat- 
ment and/or plastic deformation. Since the sulfide-stress crack- 
ing phenomenon is stress dependent, the magnitude of the ap- 
plied tensile stress should be minimized by appropriate design 
in so far as practicality will permit. 


Service Failures—Cause and Prevention 


A review of the spontaneous brittle failures that were ex- 
perienced in the development of the Pincher Creek Field will serve 
to show the nature of the failures and the means that were de- 
veloped to avoid a recurrence of them. The case histories that 
will be presented are typical of those that have been exp rienced 
in other sour-condensate fields. 

Subsurface Equipment. Tubing. uring the production testing of 
the discovery well at Pincher Creek, spontaneous brittle failure 
developed in a used string of API Grade N-80 tubing, hardened- 
steel fishing tools, and hard-drawn stainless-steel wire lines. 
While the laboratory studies of these failures were still in progress, 
two failures developed in a new string of 9 per cent nickel-steel 
tubing after 6 days’ exposure to the well fluids during the produc- 
tion testing of the second well. 

The first failure, Figs. 1 and 2, occurred at a depth of 26 ft be- 
low the surface while the well was being killed for a cement plug 
back to shut off water. During subsequent fishing operations 
the tubing again parted, this time at a depth of 1201 ft under a 
tensile load of 80,000 to 85,000 lb. This failure is shown in Fig. 3. 
A very comprehensive appraisal of the conditions attending the 
failures showed that the induced stresses were well below the 
yield strength of the material and it appeared reasonably certain 
that the failures were not mechanical in origin. It was originally 
postulated that the failures were ascribable to a form of stress- 
corrosion cracking, which were duplicated in the laboratory by 
exposing beam specimens, Fig. 4, of the 9 per cent nickel-steel 
tubing to samples of the well fluids that had been recombined in 
the laboratory to approximate reservoir conditions—5000 psi and 
190 F. Under these conditions beam specimens failed within 
1'/2 to 4 days when loaded to an outer fiber stress of 120,000 psi 
(yield strength 121,100 psi). At an applied fiber stress of 60,000 
psi, the time to failure was prolonged to 67 days. 


The sensitivity of the 9 per cent nickel-steel tubing was asso- 
ciated with the presence of retained austenite at the grain bound- 


94 / may 


1960 


a 


aries. The combined effect of plastic deformation from cold 
straightening of the pipe at the mill and the applied service 
stresses was considered to be responsible for the eventual de- 
composition of the retained austenite to martensite. The highly 
strained condition of the martensite lattice was considered the 
reason for the susceptibility of 9 per cent nickel steel. The im- 
portance of the composition of the retained austenite and its 
transformation product, martensite, on susceptibility was 
demonstrated by experimental heat-treatments involving double 
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Fig. 3 


normalizing at. 1650 and 1450 F followed by tempering at 1100 F 
for 2 hr [1,2]. The purpose of double normalizing was to mini- 
mize chemical inhomogeneities of the steel; extended tempering 
at 1100 F was intended to stabilize the retained austenite by in- 
creasing the solution of carbide. Success of this heat-treatment 
in desensitizing 4 per cent nickel steel was demonstrated by the 
fact that a beam specimen of the steel so treated did not fail after 
more than four years of exposure to a sulfide environment in the 
laboratory. However, the heat-treatment was considered imprac- 
tical from a mill-production standpoint as it required very precise 
control of tempering temperature and time. In the light of the 
above, it appeared quite definite that 9 per cent nickel steel was 
not a suitable material for sour-condensate well tubing. 
Concurrent with the failure investigation, numerous laboratory 
cracking studies were made for the purpose Of screening potential 
tubing materials prior to field testing in the flowline of one of the 
Pincher Creek wells. For the laboratory studies, an artificial 
corrodent was used, consisting of tap water saturated with hydro- 
gen sulfide and carbon dioxide. Fortunately, the results of the 


laboratory tests were in excellent agreement with those that were 
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later obtained in the field test. The results of this work gave 
further evidence that the cracking of ferritic steels was related to 
the presence of martensite in the structure. API grade J-55 
material, which usually is made from plain carbon steel (0.35-0.40 
per cent C; 1.20-1.30 per cent Mn) exhibits a pearlite-ferrite 
microstructure in the hot-rolled condition, and is completely in- 
sensitive to cracking even when it is plastically deformed to an 
appreciable extent. 

Normalized grade N-80 material exhibits a highly variable 
response to cracking, depending upon the metallurgical condition 
of the steel. This material is usually made from a carbon-man- 
ganese-molybdenum steel (0.40-0.45 per cent C; 1.5-1.75 per 
cent Mn; 0.15-0.20 per cent Mo) which is subject to segregation 
of manganese. The high hardenability imparted by manganese 
permits hardening of the segregated zones when the material is 
air cooled from the normalizing temperature. An example of 
segregated martensite, typical of that observed in many normal- 
ized N-80 products, is shown in Fig. 5. To satisfy the engineering 
requirements for the Pincher Creek wells it has been necessary to 
use N-80 tubing. For the early development wells, the tubing was 
selected from heats of steel which exhibited maximum yield 
strengths of 90,000 psi. The purpose of the limitation on 
yield strength was to more nearly insure that the chemical com- 
position of the steel would be on the low side of the usual range, so 
that the possibility of manganese segregation and potential sus- 
ceptibility would be minimized. A special heat-treatment has 
since been developed that has proved effective in rendering car- 
bon-manganese-molybdenum insensitive to cracking. 
This treatment involves tempering at 1150 F minimum after 
normalizing to develop a maximum yield strength of 90,000 
psi. However, to allow an adequate manufacturing tolerance, 
the minimum yield strength requirement for the special tubing 
was reduced to 75,000 psi. 

In addition to the regular N-80 steels, this treatment has also 


steels 


been effective in desensitizing tubing made from AISI 4140 and 
4145 steels. To date, 16 strings of special N-80 tubing, eight each 
of the carbon-manganese-molybdenum, and 4140 or 4145 steels, 
have been procured from three different pipe manufacturers. 
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The chemical and mechanical properties of all heats of steel in- 
volved are shown in Tables 1 and 2. Laboratory sulfide-stress 
cracking control tests have been made on samples representing the 
beginning, middle, and end of the tempering cycle for each heat of 
steel. All of the beam specimens were prestrained by stressing to 
the tensile strength and loaded to a maximum fiber stress of 80 
per cent of the yield strength. Under these conditions, none of 
the beam specimens have failed after exposure periods ranging 
from 125 to 235 days. 


Table 1 Range of chemical and mechanical properties of 21 heats of 
C-Mn-Mo steel for special N-80 tubing; normalized 1550 F, tempered 
1150-1175 F 


Carbon, per : 
cent 0.40-0.45 Tensile strength, 102,200-121,900 
psi 
Manganese, 


per cent 1.42-1.68 Yield strength, 


psi 75, 600-89 , 400 
Silicon, per 
cent 0.18-0.25 Elongation, per 
cent in 2 in. 21-3 
Molybdenum, 
per cent 0.15-0.23 Hardness, Rock- 18-21 
well C 


Table 2 Range of chemical and mechanical properties of 29 heats of 
AISI 4140 or AISI 4145 steel; normalized 1550 F, tempered! 150~1175 F 


Carbon, per 


cent 0.38-0.48 Tensile strength, 
psi 104, 500-124, 000 
Manganese, 
per cent 0.82-0.99 Yield strength, 
psi 77, 200-90, 000 
Silicon, per 
cent 0.25-0.30 Elongation, per 
cent in 2 in. 18-34 
Chromium, 
per cent 0.93-1.04 Hardness, 
Rockwell C 18-22 
Molybdenum, 
per cent 0.17-0.20 


Until recently, all of the tubing procured for the Pincher 
Creek wells was internally lined with a baked phenolic plastic. 
Initially, plastic linings were considered advisable primarily for 
protection of the metal during acid treatment of the wells. The 
acidizing solution contained an arsenic inhibitor, which was not 
completely effective, so that steels with questionable susceptibility 
might be damaged or conditioned for rapid failure during acidiz- 
ing. Plastic linings were never intended to provide permanent 
protection against general corrosion, since mechanical abrasion 
by bottom hole equipment and wire lines that are run through the 
tubing during well completion operations would damage the 
lining and expose bare metal to the well fluids. In addition, the 
plastics that have proved suitable for linings will crack and spall 
because of the sudden release of absorbed gas when the tubing is 
removed from the well. The success of the special heat-treat- 
ment in desensitizing the steels used for N-80 tubing indicates 
that the metal will not be damaged by exposure to the acidizing 
solution and plastic-lined tubing is no longer considered essential 
for sour-condensate wells. 

Casing. Specially heat-treated N-80 casing is not available and 
to avoid sulfide-stress cracking failures in regular N-80 casing, the 
sour well fluids must be excluded from the annulus. To minimize 
the possibility of joint-leaks in the tubing string, integral joints 
of the shouldering type that provide a metal-to-metal seal are 
recommended. The integral joint is considered superior to the 
API coupled joint in that it provides a more positive seal and 
the number of actual pipe connections are reduced by a factor of 
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two. In addition, the shouldering type joint offers the ad- 
vantage of hoop stresses that are significantly lower than those 
associated with the tapered API coupled joint. 

Bottom hole packers are used to isolate the well bore from the 
annulus. In fact it is considered good practice to use dual packers, 
so that in the event one of the packers fail, the annulus will still 
be protected. The practice of filling the annulus with a high pH 
mud has been followed as a means of affording at least partial 
protection of the casing should the well fluids infiltrate the 
annulus, because of a packer or tubing joint leak. The buffering 
action of the alkaline mud might provide adequate protection of 
the casing until such time that the leak can be repaired. Freedom 
from casing difficulties at Pincher Creek indicates that the afore- 
mentioned precautions are effective in protecting the casing string. 

Surface Equipment. The integrity of valves is of paramount im- 
portance since it is essential that they operate on demand should 
it ever become necessary to close in a well for any reason. The 
high-pressure production valves that are presently available as 
stock items all contain critical components that are made from 
susceptible materials that will fail by sulfide-stress cracking under 
circumstances where the induced service stress is of a critical 
magnitude. 

Valve Stems. In the case of gate valves, one of the most impor- 
tant problems has been providing valve stems that are immune 
not only to sulfide-stress cracking, but also to mechanical failure. 
The difficulties that have been experienced from the brittle 
spontaneous failure of hardened-ferritic-steel valve stems can be 
illustrated by reviewing one of the recent failures that developed 
at Pincher Creek. The valve stem involved was from a 4-in. 
series-2900 gate valve that had been located on the top of the 
wellhead to facilitate wire line work. After 15-days exposure to 
the well fluids, the stem failure, shown in Fig. 6, developed. The 
stem had been made from AISI 4140 steel that had been heat- 
treated to Rockwell-C 32, which is within the range Rockwell-C 
30 to 38 customarily used for high-strength valve stems. In ad- 
dition to the main fracture, numerous radial cracks were evident 
on the top flanks of many of the threads. On subsequent metallo- 
graphic examination, it was found that the radial cracks were 
associated with a shallow layer of untempered martensite as is 
evident in Fig. 7. The light etching constituent is the un- 
tempered martensite on the flanks of the threads; the darker 
etching constituent is the normal tempered martensite that de- 
velops during heat-treatment. 

The untempered martensite apparently resulted from localized 
overheating and subsequent rapid cooling at the time the threads 
were ground. In this case, the failure was ascribable to the pres- 
ence of untempered martensite in combination with high hardness 
of the 4140 steel. If, as may have been possible, the radial crack- 
ing had developed during grinding, a brittle torsional failure 
may have been precipitated under otherwise normal operating 
conditions because of subsequent embrittlement of the steel by 
absorbed hydrogen. However, if, as was most likely the case, the 
radial cracking was not due to grinding, it could have been 
initiated by sulfide-stress cracking with little or no applied stress 
because of the highly strained condition of the martensite lattice. 
Even in the absence of untempered martensite or radial cracking, 
4140-steel valve stems that are treated to hardnesses greater than 
Rockwell-C 30 would exhibit a high propensity for sulfide-stress 
cracking since they are subject to high tensile loads during normal 
valve operation and, under certain conditions, to high-impact 
tensile loads. 

Experience has shown that type 410 stainless steel, which has 
been used extensively for valve stem application, offers no ad- 
vantages over AISI 4140 steel in providing improved resistance 
to cracking. When treated to similar mechanical strength levels, 
the response of 4140 and type 410 stainless steel to sulfide-stress 
cracking is essentially the same; this is evident from the data in 
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Table 3  Sulfide-stress cracking data for materials used for valve stems 


Yield 
strength, 
psi 


Tensile 
strength, Hardness Prestress, 
psi R. psi 


Applied Time to 
stress, failure, 


psi days Heat-treatment 


AISI 4140 


147,500 36.0 151,000 


147,500 158,000 


118,000 


118,000 Normalized 1650 F; 
quenched 1500 F; 


pered 1050 F 


Type 410 


161, 400 173,200 


161,400 173,200 


K-Monel 


137, 400 
137 , 400 


38.0 
37.0 


* No failure, removed from test. 


173,700 
154,600 


129,100 


Quenched 1750 F; tempered 
750 F 
113,000 


110,000 
110,000 


1402* Cold rolled; aged for 10 hr at 
950 F 
1402¢ 


Norte: Corrodent—tap water saturated with H2S and CO:. 


Table 3 which have been selected to show the relative cracking 
resistance of materials that have been used for valve stems. A 
beam specimen of 4140 steel treated to a hardness of Rockwell-C 
36 failed in 14 days when it was prestrained by stressing to the 
‘tensile strength and loaded to a maximum fiber stress of 80 per 
cent of the yield strength. When the prestress was reduced to a 
point midway between the yield strength and ultimate strength, 
the time to failure was 17 days under the same conditions of load- 
ing. A beam specimen of type 410 stainless steel failed within 
three days when tested under comparable conditions; the time 
to failure was increased to 24 days when the applied stress was 
decreased to 70 per cent of the yield strength. 

The results of corrosion studies that were made in the flowlines 
of one of the Pincher Creek wells showed that the corrosion rate 
of type 410 stainless steel was approximately 7 mils per yr; the 
corrosion rate for carbon and low-alloy steels was in the range 0.3- 
0.4 mil per yr. In the light of these data, it is clearly evident that 
type 410 stainless steel is unsuitable for sour-condensate service 
since it does not provide any improvement in resistance to crack- 
ing and is actually less resistant to corrosion than low-alloy con- 
structural steels. 

The valve-stem problem has been resolved by substituting pre- 
cipitation-hardened K-Monel for the ferritic steels. In the fully 
hardened condition, K-Monel exhibits a hardness in the range 
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Rockwell-C 32 to 38 and is highly resistant to cracking in hydro- 
gen-sulfide environments. The data in Table 3 show that beam 
specimens of K-Monel did not fail after 1402 days’ exposure when 
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they are prestrained and loaded in the usual fashion. Despite the 
high resistance to cracking under conditions of static stress, there 
is evidence that K-Monel valve stems are susceptible to mechani- 
cal failure under circumstances where impact loads are involved. 
In this respect, the embrittlement of K-Monel as a consequence 
of absorbed hydrogen may be an important contributory factor. 
Laboratory experiments, wherein K-Monel was cathodically 
charged at a current density of 80 ma/sq in. in a 5 per cent solu- 
tion of sulfuric acid that contained a small concentration of so- 
dium arsenite, have indicated that the load carrying capacity is 
reduced by approximately 20 per cent after an exposure time of 
10 days when the applied stress is 80 per cent of the yield strength. 
Under similar conditions hardened ferritic steels fail spontane- 
ously in a matter of minutes, even when the applied stress is as low 
as 50 per cent of the yield strength. 

To minimize the possibility of tensile impact failures with K- 
Monel valve stems, the cross-sectional area should be increased 
as much as design limits will permit in order to reduce the magni- 
tude of the service stresses. When requested, one valve manu- 
facturer has provided high-pressure valves with oversized K- 
Monel stems for sour-condensate service. For example, the sub- 
stitution of a 3-in. stem for a 2-in. stem has increased the mini- 
mum area of the stem 38 per cent and decreased the stress ap- 
proximately 26 per cent. Under normal conditions this change 
has eliminated mechanical failure of K-Monel stems. However, 
under circumstances where valve operation has been impaired by 
hydrates that form in the valve body, overstressing has been re- 
sponsible for a few failures of the oversized K-Monel stems. A 
recurrence of failures of this type obviously can be and has been 
averted by heating the valves to eliminate the hydrates, so that 
the valve can be operated without difficulty. 

Valve Bodies. Although the stresses that are induced in valve 
bodies by the contained pressure are relatively low, certain design 
features and fabrication methods may be responsible for high 
local tensile stresses and damaging metallurgical structures which 
may promote sulfide-stress cracking. The difficulties that may be 
experienced under sulfide conditions from the indiscriminate use 
of valves that are made by regular production methods can be 
illustrated by the ease history of a 3-in. series-2900 wing valve 
that developed a leak in the valve body after 106 days’ exposure to 
the sour Pincher Creek well fluids. The leak occurred at a crack 
that extended into the threads of the lubrication plug hole, Fig. 8. 
It was evident by metallographic examination that the crack was 
confined to a repair weld that apparently had been made to cor- 
rect a casting defect. The nature of the repair weld is shown in 
Fig. 9, which is a photomacrograph of an etched section that was 
taken axially through the lubrication plug hole. The dark etch- 
ing portion outlining the weld is the heat-affected zone which 
exhibited a hardness of Rockwell-C 32 to 38. The weld metal ex- 
hibited a hardness in the range Rockwell-C 23 to 27, while the 
metal remotely removed from the lubrication plug hole had a 
hardness of Rockwell-C 17. The structures that were developed 
in the weldment were martensitic and there was no evidence that 
the valve body had been treated after welding to condition the 
weld. The analysis of the steel involved corresponded to a modi- 
fied AISI 8625 steel, which, because of its alloy content, possessed 
sufficient hardenability to permit hardening of the weld during 
cooling. Because of the tapered thread of the lubrication plug, 
the weldment would have been subjected to a local hoop stress 
which when added to the internal stress of the weld was sufficient 
to initiate sulfide-stress cracking of the hardened weldment. 

The special requirements that have been developed for valves 
for sour-condensate service stipulate that, in the event repair 
welding is necessary, valve bodies should be tempered so that the 
hardness of the weldment does not exceed 225 Bhn. For carbon 
steels tempering at 1150 F minimum has proved very satisfactory 
and hardnesses that are well below the maximum requirement 
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will invariably develop. In the case of alloy steels, the response 
to tempering is variable depending upon the composition of the 
steel involved. Since various manufacturers use different grades 
of steel for valve bodies, it is not possible to specify a single tem- 
pering temperature that will be effective for all bodies. It is the 
responsibility of the valve producer to select a temperature that 
will properly condition repair welds. Since the sulfide-stress 
cracking phenomenon is so vitally dependent on stress, considera- 
tion should be given to means by which highly localized tensile 
stresses in valve bodies, such as those resulting from the make-up 
of tapered lubrication plugs, can be averted. 

Valves of the lip-seal construction, which ordinarily are in- 
tended for high-temperature-steam service, have proved unsatis- 
factory in high-pressure sour-condensate service. During the 
first six weeks of operation of the gas-cleaning plant at Pincher 
Creek, nine of approximately 120 lip-seal valves, that had been 
installed in the sour-gas section of the plant, failed. In each in- 
stance the failure originated in the seal weld between the bonnet 
and the valve body, Fig. 10. The length of the visible cracking 
varied from '/s in. in some of the failures up to a maximum of 
2'/. in. in others. The body and bonnet material was 1.25 per 
cent chromium-0.5 per cent molybdenum steel which as ordi- 
narily processed developed a hardness in the range 141 to 175 
Bhn. The seal welds had been made by the Heliare process. In 
some valves the weld had been made without the use of a filler 
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rod, while in others a filler rod of 1.25 per cent chromium steel 
had been used. The metallurgical examination indicated that the 
weld failures were ascribable to the combined effect of me- 
chanically deficient welds and sulfide-stress cracking. The 
mechanical strength of the weld would have been depreciated by 
extensive porosity which had developed in the weld metal as a 
consequence of thermal decomposition of a graphitic lubricant 
that had been used to facilitate assembly of the valves. This was 
indicated by the presence of graphite particles on the shoulders of 
the bonnet and valve body, Fig. 11. The voids, which in some 
instances penetrated approximately two-thirds of the weld thick- 
ness, were filled with decomposed lubricant as evident in Fig. 12. 
In each of the failures the weld metal and heat-affected zone ex- 
hibited hardnesses in the range Rockwell-C 30 to 39 and as such 
were in a susceptible metallurgical condition for sulfide-stress 
cracking. The stress concentrating effect of the weld porosity 
intensified the service stresses and probably contributed to the 
rapid failures that developed. However, under the circumstances 
involved, it is believed that sulfide-stress cracking would have 
eventually developed even if the welds had been mechanically 
sound. Because of the manner in which the lip-seal valves were 
constructed, it was not possible to thermally process the welds to 
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render them insensitive to sulfide-stress cracking; the heat at- 
tending tempering would have seriously impaired the operating 
characteristics of the valves. The questionable susceptibility of 
untempered welds, particularly in alloy steels, indicates that 
valves of the lip-seal weld construction are unsatisfactory for 
high-pressure sour-condensate service. The corrective measures 
employed in the already discussed situation, involved replacement 
of the questionable lip-seal valves with plain carbon-steel valves 
having bolted bonnets. These valves have been in service for 
over two years and have given trouble-free service. 


Conclusions 


Despite the fact that the present knowledge of the sulfide- 
stress cracking problem has not provided a satisfactory explana- 
tion of the failure mechanism, the information that has been ac- 
cumulated by experimentation and actual service experience has 
provided means of averting spontaneous failure in sour-conden- 
sate well equipment. The discussion that has been presented in- 
dicates that freedom from service failures can be achieved by 
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providing materials that are metallurgically satisfactory and by 
appropriate design, maintaining the applied stresses at as low a 
level as it is practical. The present concepts of the requirements 
for sour-condensate well equipment indicate the following ma- 
terials considerations: 


1 When well conditions permit, grade J-55 tubing made from 
plain carbon steel should be used. If the engineering require- 
ments are such that J-55 is inadequate, specially processed N-80 
is indicated. Experience has shown that the carbon-manganese- 
molybdenum steels and AISI 4140 or 4145 can be rendered in- 
sensitive to sulfide-stress cracking by tempering at a minimum 
temperature of 1150 F to develop yield strengths in the range 75,- 
000 psi minimum-90,000 psi maximum. Specially treated N-80 
casing is not now available and to avert sulfide-stress cracking 
failures in potentially susceptible N-80 materials, it is necessary 
to exclude the well fluids from the casing-tubing annulus. This 
can best be accomplished by using bottom hole packers and by 
insuring against leaks in the tubing joints. 

2 Critical components of wellhead equipment that require 
materials with hardnesses of Rockwell-C 30 or greater should be 
made of K-Monel, which under the conditions encountered in 
sour-condensste well service is immune to cracking. 

3 Equipment that is fabricated by welding or which may be 
repaired by welding should be thermally processed after welding 
to limit hardness to 225 Bhn maximum. For carbon steels, this 
can be accomplished by tempering at 1150 F; because of the 
variable response of alloy steels to tempering, it is not possible to 
specify the temperature at which the maximum hardness re- 
quirement will be achieved. Therefore the equipment supplier 
should establish the minimum tempering temperature required 
to properly condition the particular type of steel that may be 
involved. 
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DISCUSSION 
J. P. Fraser? 


The authors have done an excellent job of presenting the 
methods which are now available for avoiding cracking-type fail- 
ures in sour-gas well environments, and I agree with their general 
approach. If equipment for use in sour gas wells is specified in 
accordance with their suggestions, there is little likelihood of ob- 
serving sulfide-stress cracking.* Nevertheless, I would like to 
differ with the authors on a few minor points and to ask several 
questions. 

First, I would like to ask how the authors determine when it is 
necessary to use the precautions and special materials specifica- 
tions which they describe. Such precautions and specifications 
are useful if there is risk of sulfide-stress cracking, but they do cost 
money. 

Second, what do the authors mean by “critical stress necessary 
for sulfide-stress cracking?’ This is a term which we have used 
(Reference [9]) to denote the stress for 50 per cent failure 
probability in the given corrosive medium. I presume that the 
authors mean something which could be termed a threshold 
stress, below which failure will not occur. Is there any assurance 
that such a threshold stress exists for a steel (e.g., API grade N-80 
steel) which is susceptible to cracking? Our own work indicates 
that there may not be such a threshold (see Reference [9] ), al- 
though the probability of cracking diminishes rapidly as the stress 
is decreased below our critical stress for cracking. 

The authors have stated that the load carrying ability of K- 
Monel is lower after cathodic charging than before. Does this 
mean that they observed cracking-type failures or was the ulti- 
mate tensile strength reduced? The authors further state that 
K-Monel suffers from loss of impact strength, and they attribute 
a part of this to occluded hydrogen. I would be very much sur- 
prised to find that occluded hydrogen affects the impact strength 
of any alloy. For example, it has been demonstrated‘ with hy- 
drogen charged steel that no loss of ductility is apparent at very 
high strain rates; embrittlement only shows up at moderate and 
slow strain rates. This is apparently a function of the diffusion 
rate of hydrogen in steel which although fast is still finite. In 
common with the authors, we regard K-Monel as one of the most 
sulfide-stress cracking resistant of all engineering materials; no 
failures have occurred in laboratory or field tests using stressed, 
hardened K-Monel exposed to hydrogen sulfide bearing environ- 
ments. The only conditions under which we have observed 
cracking failures have been those described by Bergman and 
Bjorkman,’ which involved a mixture of moist crvstals of NH,Cl 
+ NaS. 

2 Shell Development Company, Emeryville, Calif. 

3 Note: I prefer to use the term “sulfide corrosion cracking’’ which 
is the nomenclature used in most of the reports of NACE Committee 
T-1G, which studied this problem for several years. My preference 
for this term, rather than “‘sulfide-stress cracking’ as used by Mec- 
Guire and Vollmer, stems from my belief that the mechanism of 
cracking-type failures in the presence of hydrogen sulfide is basically 
stress corrosion cracking. It is my belief that occluded hydrogen 
plays a secondary role in crack propagation, not a primary one. My 
reasons for this belief are set forth in Reference [2] of this paper. 

4T. Toh and W. M. Baldwin, ‘‘Ductility of Steel With Varying Con- 
centrations of Hydrogen,” Stress Corrosion Cracking and Embrittle- 
ment (Robertson, Editor), John Wiley & Sons, Inc., New York, N. Y., 
1956, pp. 176-186. 

6 D. J. Bergman and 8. A. Bjorkman, ‘‘Some Typical Cases of 
Corrosion in Platforming Units,’’ API Division Refining, Phila- 
delphia, Pennsylvania (May 13, 1957). 
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J. F. Mason, 


The experiences related by the authors should prove very help- 
ful to other operators in the Pincher Creek as well as other sour- 
condensate well areas, and they should be complimented on the 
excellent and concise manner in which they described the metal- 
lurgical methods and design changes used in combating sulfide- 
stress cracking of equipment. 

Other than the theoretical discussion of the metallurgy of 
9 per cent nickel steel which is open to further clarification at 
some later date when current extensive laboratory studies are 
completed, there are several other methods that might be con- 
sidered for the mitigation of problems encountered in sour gas 
condensate wells. 

Mention was made of the spontaneous brittle failure of hard- 
ened steel fishing tools and hard drawn stainless steel wire lines. 
In the first instance, it is apparent that failure was induced by 
sulfide-stress cracking, while in the latter case the indications 
are that chloride-stress cracking was involved. Under either 
circumstance it should be noted that the situation can be cor- 
rected by the use of Monel and K-Monel. Both of these materials 
possess excellent resistance to sulfide and chloride stress corrosion 
cracking and as a matter of fact they have been adopted on a 
limited scale for fishing tools and wire lines. Because of its 
higher strength and hardness due to thermal treatment, K-Monel 
has been preferred for the fishing tools while regular Monel has 
proved adequate for wire line service. 

It is fortunate that a special heat-treatment has been de- 
veloped that is effective in rendering carbon-manganese-molyb- 
denum steels insensitive to cracking at yield strengths below 90,- 
000 psi. However, under circumstances where higher mechanical 
properties are required, the use of Gray Monel-lined tubing would 
appear to be a practical solution to the problem of sulfide-stress 
cracking. This product was developed by the Gray Tool Com- 
pany, Houston, Texas, and involves the hydrostatic expansion of 
a seamless, light gage Monel tube inside a joint of steel tubing. 
During this operation the Monel yields slightly. As the pressure 
is relieved, contraction of the steel actually leaves the Monel in 
compression, resulting in a highly efficient mechanical bond to 
the steel. This hydrostatic-forming operation extends beyond the 
tapered sealing surface on each end. When the string is made up 
with a Graloc joint including a K-Monel sealing ring, this pro- 
vides a continuous corrosion resistant internal surface in the 
tubing. 

An experimental section of this duplex product made up of two 
30-ft lengths of 3!/2 in. OD N-80 oil well tubing equipped with 
Graloe joints and lined with 0.035-in. wall Monel tubing was 
installed in a flow line of a well in the Pincher Creek gas field 
handling the following conditions: 


Line pressure........... 1500 psi 

Flow through tubing... . 13,000,000 cu ft/day 
Hydrogen sulfide. . : 10 per cent 

Carbon dioxide...... 6 per cent 
Temperature of gas...... 135 F 

Water in gas...... 6 bbl/1 million cu ft 
Depth of well........ 13,000 ft 


After six-months’ service the section was removed from the 
flow line and visually examined. The Monel liner was in excel- 
lent condition, being free of significant corrosion and free of bulg- 
ing, blistering, or cracking. This test is not conclusive because of 
its relatively short duration and the tubing has been returned to 
service for further testing. However, the results to date do indi- 
cate that this type of product may provide a successful means for 
handling sour gas condensate wells where sulfide-stress cracking 
of high-strength steels may be a serious problem. 


6 Supervisor, Development and Research Division, The Inter- 
national Nickel Company, New York, N. Y. Mem. ASME. 
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R. S. Treseder’ 


The authors are to be complimented on the clear and concise 
way in which they have presented a complex subject, one which 
has a number of controversial aspects. The work we have done 
in the laboratory and our field experience have led us to con- 
clusions which are in substantial agreement to those of the 
authors. The following comments are concerned with a few of 
the secondary points in the paper. 

No mention was made of the use of inhibitors, and it would be 
of interest to have the authors’ current views on this subject. We 
have thought it desirable to use inhibitors in sour-condensate wells 
even though all materials used were specified to be those with 
adequate resistance to cracking in the particular well environ- 
ment. In our experience such inhibitors are not required from 
the standpoint of preventing metal loss, but only as insurance 
against cracking failures arising from off-specification materials 
that might be present in the tubing string or wellhead owing to 
the inadequacies of inspection procedures. 

In several places the authors mention that certain materials, 
e.g., J-55, tempered N-80, etc., are insensitive to sulfide-stress 
cracking. I believe this statement should be modified in order to 
be strictly correct. I would agree with this statement only if the 
added qualification, ‘under environmental conditions so far en- 
countered in sour-condensate wells” is stated or implied. Under 
laboratory conditions more severe than any well conditions we 
have yet encountered, we have been able to obtain sulfide-stress 
cracking of J-55 and tempered N-80; in both cases the material 
was cold worked. This has been reported in paper [10]. An 
added point to be made here is that condensate well environments 
do vary in severity as regards sulfide-stress cracking. From the 
information reported I would estimate that Pincher Creek was of 
the same order of severity as in Jumping Pound Field in Alberta; 
the latter has represented the most severe well environment that 
we have encountered. 

As a final point, it would be of interest to know more about the 
nature of the hydrogen attack experienced with K-Monel. For 
example, on page 98, it is mentioned that the load carrying ability 
of K-Monel was reduced 20 per cent following cathodic charging; 
was this reduction determined by a standard tensile test? 


Author’s Closure 


In reply to Mr. Treseder’s question concerning the use of inhibi- 
tors, I would like to state that we are not using chemical in- 
hibition as a means of preventing cracking, but have relied upon 
the adequacy of the special heat-treatment to render the steels in- 
sensitive. The effectiveness of inhibitors is vitally dependent on 
uninterrupted treatment to maintain the concentration of the 
inhibitor at an effective level. Occasional interruption of the in- 
hibition treatment cannot be avoided, particularly in Western 
Canada where severe winter weather conditions make access to 
the wells impossible at times. Under such circumstances sulfide- 
stress cracking may, and in fact did, develop when the concen- 
tration of inhibitor falls below the minimum required for complete 
protection. We are hopeful that the practice we have followed, 
involving sulfide-stress cracking control tests on samples of the 
tubing from each heat of steel representing the beginning, center, 
and end of the tempering cycle, will minimize the possibility of 
susceptible material being placed in service. 

We are in complete agreement with Mr. Treseder’s comment 
that the preventive measures presented in the paper should apply 
only to the environmental conditions that so far have been en- 
countered. At this time we consider each new sour-condensate 
field potentially dangerous and have recommended the use of 
special materials until such time as experience dictates otherwise. 
In our opinion the severity of unknown environments can best be 
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established by field tests involving the exposure of susceptible 
materials to the sour well fluids. In the more severe environments 
that We have had to cope with, the precautions indicated have 
been adequate. In situations where the hydrogen sulfide con- 
tent is low, for example, 1'/2 per cent in the Fashing Field, the 
initial wells have been equipped with special tubing and acces- 
sory equipment asa precautionary measure. We agree with Mr. 
Fraser that the special precautions are expensive, but in the ab- 
sence of a reliable means of predicting the material requirements 
in unknown situations, the expense is not unwarranted. 

The cathodic charging test on K-Monel involved the exposure 
of 0.125-in. diameter specimens at a current density of 80 ma/sq 
in. under a static tensile load equivalent to 80 per cent of the yield 
strength. The high current density used to accelerate the labo- 
ratory studies was substantially greater than that associated 
with severe hydrogen sulfide corrosion which would approximate 
0.25 ma/sq in. The electrolyte used for these tests was a 5 per 
cent H2SQO, solution which contained 20 mg/liter of AsoO,. After 
10 days’ exposure without cracking, the specimens were pulled to 
destruction. The reduction in load carrying ability questioned 
by Mr. Treseder and Mr. Fraser represents the reduction in the 
tensile strength of the material. We did not, as Mr. Fraser 
indicates, state that K-Monel suffers from loss of impact strength 
as a consequence of occluded hydrogen. The question apparently 
is directed to that . there is evidence 
that K-Monel valve stems are susceptible to mechanical failure 
The 


complex stress system and the moderately fast strain rates in- 


the statement 
under circumstances where impact loads are involved.” 


duced in the stem when the valve is snapped into the full open 
or closed position are significantly different than those induced 
in the high speed tensile impact and notch-bar impact. test 
specimens referred to by Mr. Fraser. 

Preliminary cathodic charging studies of stressed K-Monel 
show that the normal tensile properties are not recovered when 
the material is aged to eliminate the effect of absorbed nascent 
hydrogen. When tested under the charging conditions referred 
to above, the tensile strength of K-Monel was reduced from 
165,000 to 122,500 psi after an exposure period of 10 days. In 
this case the specimen was pulled to destruction without removal 
from the environment. Another specimen, which had been aged 
for one week at room temperature after charging for 10 days, 
failed at 130,000 psi. This indicates that the loss in load carrying 
capacity may be, at least in part, permanent, and, if this is true, 
it is quite conceivable that the impact characteristics of K-Monel 
may be influenced by occluded hydrogen. 


More experimental 
vork.will be necessary to establish whether or not this is so. 
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The critical stress necessary for sulfide-stress cracking may be 
considered similar to the endurance limit observed in fatigue test- 
ing. It should be emphasized that the critical stress is the total 
effective tensile stress to which the material is subjected. Under 
the test conditions that we have elected to use; that is, tap water 
saturated with hydrogen sulfide and carbon dioxide, steels that 
are heat-treated to yield strengths in the range of 90,000—110,000 
psi can sustain applied loads in the range of 60-65 per cent of the 
yield strength. When these materials are treated to the yield 
strengths in the range 90,000 psi or less, failure does not develop 
when the applied stress is 80 per cent of the yield strength. | Per- 
haps the apparent absence of the threshold (critical) stress re- 
ferred to by Mr. Fraser is ascribable to the presence of acetic 
acid in the corrodent used for his studies. 

Mr Mason has indicated that the failure of the hard drawn, 
stainless steel wire lines in sour-condensate wells involved chlo- 
ride-stress cracking. While we agree the presence of the chloride 
ion will certainly promote stress-corrosion cracking of stainless 
steel, we are not at all certain that the failures experienced at 
Pincher Creek by us were of this nature. In fact, failures were 
experienced in lines made of Types 316 and 304 stainless steel 
under circumstances where the presence of the chloride ion was 
very remote; i.e., in flowing wells that had not been acidized 
with hydrochloric acid. failures to 
Much of the wire line work that has 
been subsequently done at Pincher Creek involved the use of 
carbon steel wire. To prevent failure, the tubing was filled with 
inhibited diesel oil. 


We have ascribed these 


sulfide-stress cracking. 


We agree that under circumstances where 
the wire line is in contact with the sour well fluids, regular Monel, 
which presumably had been cold drawn to provide the desired 
mechanical strength, would be much more satisfactory than the 
carbon steel or stainless steel. 

In conclusion, we would like to supplement the information on 
the performance of the Gray Monel-lined tubing reported by Mr. 
Mason. Several years ago a sample length of this tubing was 
installed in the flowline in the Stockham-Federal Wellin Wyoming 
and handled the production of the sour well fluids for approxi- 
mately six months. A very careful examination showed the sam- 
ple to be in excellent condition without anv evidence of corrosion, 
either of the general or pitting type and there was no indication 
that the Monel lining had blistered from entrapment of hydrogen 
at the junction of the Monel and steel. We always have ex- 
pressed the belief that the Gray Monel-lined tubing would be 
suitable under severe conditions, provided the additional expense 
of the Monel lining can be justified. 
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Bones [1]! and by Marin and Rimrott [2]. 


from both torsion and tensile test specimens. 


theories. 


properties throughout and as nearly perfect as possible. 


shearing of steel slabs, 6 to 10"/2 in. thick. 


coincides to a transverse direction in the ingot. 


sections were forged to 2-in. diameter round bars. 


of commercial forged bars of these steels, 


from which the cylinders were machined. 


Petroleum-Mechanical Engineering Conference, 


ENGINEERS. 


Society. 
1959. 


Paper No. 59—Pet-1. 


Nomenclature 


based on both tension and torsion data. 
sures 1s described and checked against the experimental data. 


Sa BEHAVIOR of thick-walled cylinders subjected to 
stresses in the plastic range has lately received renewed attention 
and various new working theories have been suggested, of which 
excellent reviews can be found in recent papers by Crossland and 
In reference [3] a 

theory was proposed based on the stress-strain relation of the 
; material both in the elastic and the plastic ranges as derived 
A number of 
cylinders were pressurized to destruction and compared with the 
The materials of the cylinders were English steels 
and extreme care was taken to obtain materials of uniform 


On the basis of the very encouraging results of these tests it 
was decided to continue the tests using, for test material, four 
American commercial steels to specifications described in Table 1. 

Sufficient stock of each of the four steels was obtained by hot 
The stock was sheared 
off in 6-in. lengths, so that the long dimension of the stock 
The slab 
In view of the 
fact that the axes of the forged bars do not coincide with the 
axes of the ingots, the bars cannot be considered representative 


All four materials were tested in the as-received condition. 
Both torsion and tension specimens were positioned at stations 
along the bars, between and alternating with the sections, 
Previous tests showed 
no significant differences in results between torsion specimens 
cut in the longitudinal and transverse directions in the bar, so 
! Numbers in brackets designate References at end of paper. 
Contributed by the Petroleum Division and presented at the 


Houston, Texas, 
September 20-23, 1959, of THe AMERICAN Society OF MECHANICAL 


Note: Statements and opinions advanced in papers are to be un- 
derstood as individual expressions of their authors and not those of the 
Manuscript received at ASME Headquarters, June 16, 


thick-walled, closed-end cylinders of four different steels. 
through the plastic and strain-hardening ranges to destruction, at pressures up to 100,000 
A theoretical method is given for computing expansion and bursting strength, 
Finally, a simplified formula for ultimate pres- 


Table 1 


USS T-1 


Grade 


Heat-treatment Quench 
and 
temper 
Austenitizing tempera- 
ture, deg C:...... 900 
Tempering tempera- 
ture, deg C. 610 
Carbon, per cent. 0.15 
Manganese, percent... 0.87 
Phosphorus, per cent.. 0.015 
Sulfur, per cent O15 
Silicon, per cent... . 0.24 
Nickel, per cent....... 0.91 
Chromium, per cent... 0.50 
Molybdenum, per cent. 0.47 
Vanadium, per cent... 0.04 
Boron, per cent ... 0.003 
Copper, per cent... . 0.38 
90 ,000 
Ultimate, psi 105,000 


Overstrain Tests on Thick-Walled Cylinders 


The experimental work described in the following covers pressure-expanstion tests on 


The tests covered stress levels 


Materials® and position of specimen 


HY-80 A-302-B A-212-B 
Normal- Normal- Normal- 
izeand izeand ize and 
temper temper temper 
850 850 850 
630 610 600 
0.15 0.23 0.27 
0.25 1.00 0.82 
0.014 0.018 0.022 
0.020 0.015 0.032 
0.20 0.18 0.20 
2.50 0.17 0.03 
1.19 0.15 
0.30 0.48 
0.06 0.07 0.05 
80,000 50,000 38,000 
95,000 80,000 70,000 


The manufacturer's specifications are given. 


transverse torsion specimens only were tested, except for the 


A-212 material, which was tested in both directions. 
Tension specimens in both directions were produced and 


tested. 


line of the bars. 


across the cross section of the bar. 
In reference 


high. 
the following. 


For both tension and torsion all transverse specimens 
were located 0.32 in. and all longitudinal 0.37 in. from the center- 


Several hardness pieces were also cut from the bars and the 
transverse surface ground to enable the hardness to be measured 


[3], the test results, when correlated with the 
theories, showed the theory based on torsion superior to the 
theory based on tensile test specimens, which were consistently 
An empirical correction factor is therefore introduced in 
In the present paper, the equivalent stress at ra- 
dius r, where the strain equals u/r, is assumed to equal the stress 
on a conventional stress-strain curve, corresponding to strain 


initial inner radius I 


b = initial outer radius , 
r = initial intermediate radius 
r, = initial radius at yield . 
r, = initial radius at end of per- ; 
fectly plastic zone 
Tmy?, = initial intermediate radii 

P = internal pressure 
P,, P; = pressure summation Fy 
AP,, AP, = increments of pressure sum- ou 
mation € 
P, = pressure summation for €, 
elastic zone € 
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pressure summation for Ug 
perfectly plastic zone Uy 

equivalent true stress 

equivalent true stress at Um Un 
outside radius 

equivalent true inter- 
mediate stresses Ty 

true yield stress ps 

true ultimate tensile stress K, 

actual radial strain 

hoop strain m 

axial strain E 


radial deformation at bore 

radial deformation at radius 
lp 

intermediate 
formation 

shear stress 

yield shear stress 

shear strain 

initial-diameter ratio 

actual or strained-diameter 
ratio 

Poisson’s ratio 

Young’s modulus 


radial  de- 
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u/r. In reference [3], the corresponding strain is taken from an 
equation which, for small strains, gives the strain value of (2/+/3) 
(u/r). When the test results in feference [3] were rechecked 
on the basis of using strain u/r, much better agreement was 
found with the test results, and this approach is therefore used 
in the present paper. 

Due to this change it became necessary to change the tabulation 
which is used to facilitate the arithmetical integration of the 
pressure-expansion curve. The new table, as shown in the 
appendix, has also been revised in case of the torsion theory 
so as to facilitate its use and cut down the number of calculations. 

The development of the basic methods was given in references 
[1] and [3], and only the main formulas are given in the ap- 
pendix. 


Apparatus and Techniques of Testing 


The apparatus used for tension, torsion, and pressure tests 
is fully described in references [1] and [3], and consequently 
no details are given here, except a brief description. 

In the torsion and tension tests the machines used have a 
stiff load bar placed in tandem with the specimen. With such 
machines it has been found, for steels, that the load and strain 
readings in the plastic region rapidly reach steady values after 
an increment of strain has been applied. Consequently, the 
test procedure, after the elastic behavior and yield point 
had been determined, was to apply increments of strain and 
measure the load and strain after a delay of three minutes. 

For the pressure tests small increments of pressure were applied 
in the elastic region, and readings of pressure and strain noted 
after three minutes to allow the slight temperature increase, 
caused by the compression of the oil, to be dissipated. In the 
post-yield and plastic regions, pressure was maintained constant 
for a short period and then the intensifier was switched off and 
the readings taken after a further three minutes. 

The tests were continued until the ultimate pressure was 


180,000 


160,000 


passed and some of the specimens were then burst without 
further measurements. 


Results 


Tension Tests. The results of the tension tests are shown in 
Fig. 1, which shows the nominal s‘ress-strain curves of the four 
steels, and Figs. 2, 3, 4, and 5, ‘u which the true stress at the 
minimum diameter is plotted against the strain at the minimum 
diameter. At large strains the true stress at the neck has been 
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Fig. 2 T-1 tension tests 
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Fig.3 HY-80 tension tests 
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determined using the analysis due to Bridgman [4], projecting 
an image of the neck and measuring the radii of the neck profile. 
All the stress-strain curves have scatter bands in the post-yield 
and plastic ranges, but it should be noted, that for each single 
specimen the experimental points lie on a smooth curve. 

Average values of the ultimate strength and yield stresses 
are given in Table 2. Values of Young’s modulus were found 
but are not very reliable as the specimens tested were small with a 
short gage length, so that the accuracy of strain measurement 
in the elastic region was not high. However, reliable values for 
Young’s modulus were found from the pressure tests. 
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Fig. 4 A-302 tension tests 


Table 2 Average values of ultimate strength and yield stresses 
T-1 HY-80 A-302 A-212 
Elastic limit, psi. ..... 106,000 67,000 
Upper yield, psi....... 55,000 37,000 
Lower yield, psi....... 53,000 32,300 
Ultimate tensile, psi... 126,000 91,000 78,000 63,000 


Torsion Tests. The results of the torsion tests are shown in 
Figs. 6, 7, 8, and 9, as indicated by the solid-line curves, from 
which the shear stress-strain curves have been derived using the 
method given by Nadai [5]. As in the case of the tension tests 
there is a considerable scatter of points in the post-yield and 
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Fig.7 HY-80 torsion test and derived shear-strain curves 
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Fig. 8 A-302 torsion test and derived shear-strain curves 
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Fig.9 A-212 torsion tests and derived shear-strain curves 


plastic ranges. The average values of the physical properties 
for the four steels are given in Table 3 


Table 3 Average values of physical properties 


T-1 HY-80 A- 302 

106 106 
60,500 38,000 


Material, 
Shear modulus, psi. 


Elastic limit, psi... .. 
Upper yield, psi. 


23-26 000 
Lower vield, psi. 


16,600 
Pressure Tests. The theoretical pressure/expansion curves at 
small strains, based on both tension and torsion data are shown in 
Figs. 10, 12, 14, and 16, which also show the experimental 
points. It should be noted that, after initial vielding had oe- 
curred at the inside diameter, and the internal pressure continued 
to rise, there were large differences between the diametral strains 
measured on two mutually perpendicular diameters in the 
central plane, and in the graphs the average of these has been 
plotted. 

In previous tests [3], it was noted that unsymmetrical ex- 
pansion occurred in the partially-plastic range only for a soft 
steel having a drop of stress at yield. That it also occurred in 
these tests for steels like T-1 and HY-80, which show no drop in 
yield stress, must be attributed to other causes, most likely 
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anisotropy due to the method of fabrication of the steel bars. 
It was noted too, when the bars were cut up and machined, that 
they contained an unusual high amount of slag inclusions, voids 
and other faults. 

From the elastic portion of the curves the elastic constants 
may be calculated and mean values of these quantities are given 
in Table 4 


Table 4 Mean values of the elastic constants 

HY-SO A-302 

29:6 30.0 x 
106 106 

288 0.282 


Material . 


Young’s modulus, psi. 


0.286 


Poisson’s ratio 

The theoretica! pressure/expansion curves at large strains, 
based on tension and torsion data are shown in Figs. 11, 13, 15, 
and 17 together with the experimental points. The large dif- 
ferences between the diametral strains on mutually perpendicular 
diameters, which were noted at small strains, were maintained at 
large strains, and local bulging occurred well before the peak 
pressure was reached. In many of the cylinders the bulging did 
not occur in the central plane in which the measurements were 
made. This is contrary to findings in previous tests, references 
{1] and [3], in which the cylinders expanded symmetrically 
around the axes and uniformly along the gage length until the 
ultimate pressure was reached. This is best explained by as- 
107 


MAY 1960 / 


| | | 
| | 
100,000 
| 
= 
| To 007 | | 
ox 
RES 
20,000 
| | 
4 / | | | | | | | Baer 
‘ 0,02 0.04 0.05 0.06 
be 
iy | | | | 
| 
j | 
0.1 0. 0. 04 Os 13 
A-212 
106 
A-212 
30.4 X 
| Be 


100,000 


| } = 


1K21.25 


| 
— —|—-2 — — 


PRESSURE LBS /IN@ 


20,000 
THEORETICAL CURVES: 

| ———— TENSION DATA 
----- TORSION DATA 
© EXPERIMENTAL POINTS 


| 
| 


0.001 oo002 0003 0004 0005 


DIAMETRAL STRAIN . AT EXT. SURFACE 
Fig. 10 T-1 pressure/expansion curves at small strains 


100,000 


60,000 


PRESSURE LBS/IN® 


Kel.25 


THEORETICAL CURVES: 


————— TENSION DATA 

o EXPERIMENTAL POINTS 
| | 


° 0.01 002 0.03 004 


DIAMETRAL STRAIN €, AT EXT. SURFACE 


Fig. 11 T-1 pressure/expansion curves at large strains 


suming that an abnormally high degree of anisotropy existed in 


the materials, and this also accounts for the comparatively 
larger deviations between theoretical and experimental values 
than were found in previous tests, where unusual care was 


taken to procure test specimens as nearly perfect as practically 
possible. 


As explained hereafter, three additional cylinders were tested, 
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Fig. 13 HY-80 pressure/expansion curves at large strains 


two of HY-80 and one of A-302 steel. Their pressure-expansion 
curves are not included because the diametral ratios were close 
to ratio of the heaviest cylinder of each of the two materials, 
with the result that the test points were very close together and 
the theoretical curves overlapping. 

Normally the expansion curves for large strains when plotted 
for all pressures up to the ultimate, will be shorter as the cylinders 
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Fig. 14 A-302 pressure/expansion curves at small strains 
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Fig. 15 A-302 pressure/expansion curves at large strains 


get heavier. This is very pronounced for the A-212 steel, Fig. 
17. The reason for this is that, at ultimate pressure, the stresses 
throughout a thin-walled cylinder will be close to the ultimate 
tensile stress. For very heavy cylinders only the stresses in 
the inner layers will be close to the ultimate, the stresses in the 
outer layers will be much lower, which accounts for the relatively 
shorter outside strains at ultimate pressure. 
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Fig. 16 A-212 pressure/expansion curves at small strains 
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Fig. 17 A-212 pressure/expansion curves at large strains 


It will be noticed that this effect is not very obvious for the 
other steels in these tests. As a matter of fact T-1 and A-302 
on Figs. 11 and 15 seem to show that just the opposite is true. 
This again must be attributed to anisotropy. These two ma- 
terials show a very large range in ductility, as indicated by the 
reduction in area of the tensile specimens as follows: 
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Ultimate Pressures. The ultimate pressures plotted against the 
diameter ratios are shown in Figs. 18, 19, 20, and 21, including 
the theoretical curves based on both tension and torsion data. 


In all cases, the curves based on tension data lie consistently 
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Fig. 18 T-1 ultimate pressure tests 
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Fig. 19 HY-80 ultimate pressure tests 


110 | may 1960 


higher than those based on torsion data, and apparently seem to 
be in better agreement with the test data. However, this must 
be considered as entirely fortuitous, in view of the anisotropy 
and extreme variations in properties of all four materials. 

When it was found that the test results for the cylinders of 
HY-80 and A-302 seemed particularly high, two further cylinders, 
one of HY-80, and two of A-302, were tested as there was some 
slight doubt as to whether there was an error in pressure measure- 
ment. However, these tests removed any doubt, as they con- 


— 


80,000 } —— 


° 
8 
| 


ULTIMATE PRESSURE LBS /IN@ 


+— = 
THEORETICAL CURVES: 
TENSION DATA 
------ TORSION DATA 
© EXPERIMENTAL POINTS 


2 3 


DIAMETER RATIO K 
Fig. 20 A-302 ultimate pressure tests 
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Fig. 21 A-212 ultimate pressure tests 
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Table 6 Pressure summation for T-1 steel (stresses and pressures in 1000 Ib) 


Tension method 


- Low ————-High Low High 
1 2 3 4 5 6 7 8 9 10 11 12 13 
r o APs re o APs Ps T AP; P; T AP; P; 
2.8 140 279 149 306 73 276 80 308 
8.3 8.5 ca 8.2 
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firmed the previous results.’ They are included in the graphs of 
the ultimate tests, but not in the graphs of the expansion curves, 
as already mentioned, as they overlapped previous tests and 
their inclusions would have made the graphs less legible. 

Hardness Tests. Hardness checks were made at various stations 
along each bar. At each «station the hardness along three 
equally placed diameters were measured. These hardness 
checks confirmed the variation in properties along the axis of 
the bar established by torsion and tension tests. They also 
showed for HY-80 and A-302 that the properties were not sym- 
metrical about the axis of the bars, one diameter showing noticea- 
bly lower hardness figures than the other two 

Simplified Formula. Substituting stresses and strains from stress- 
strain curves, Figs. 2, 3, 4, and 5, into formulas (12), (13), and 
(14), the ultimate pressures have been estimated for each test 
cylinder, and compared with the test results, as shown in Table 


7. 
It will be noted that the ultimate pressures derived from the 
simplified formulas are almost identical with those derived from 
the pressure-expansion summation curves, as indicated by a 
comparision with Figs. 18, 19, 20, and 21. 
Form of Fractures. Typical fractures are shown on photographs 
of cylinder cross sections, Figs. 24, 25, 26, and 27. The T-1 and 
HY-80 steels show shear fracture to be predominant, whereas the 
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Torsion method 


80 


Table 7 Ultimate pressures, estimated and test 


Diameter --——Ultimate pressures, psi-—-— 

ratio, -——Theoretical—— Experi- 

Material K High Low mental 
T-1 1.25 31,200 29,800 29,400 
1.50 57 ,000 53 , 500 53 , 800 

Lito 79,400 74,000 75,000 

2.00 98 , 500 91,800 93 , 500 

HY-80 1.25 22,200 20 ,000 20 200 
1.80 58 ,000 52,600 53,000 

2.40 86 ,000 78,600 81,200 

2.76 100 , 400 90 , 000 97 , 100 

2.81 102,200 90 , 400 99 , 600 

3.00 107 , 200 97 ,600 100 , 500 

A-302 1.25 18,100 16,300 13,500 
2.00 57 ,400 51,800 51,700 

2.76 83 , 500 77,500 78,200 

3.29 96 , OVO 87 ,500 87,700 

3.44 100 , 200 91,300 98 , 500 

A-212 1.50 26,100 24,400 23,800 
2.75 61,900 56 , 900 59 , 000 

3.925 79, 700 74,000 79,500 

5.00 93 , 400 86,700 90 , 500 


fractures of A-302 and A-212 seem to be more like radial cleav- 
ages. No conclusions on the mechanism of failure have been 
drawn. 
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APPENDIX 
Basic Formulas 


The formulas used in this paper and the simplifying assumptions 
used as basis for their derivations have been described in previous 
reports, see [1] and [3], therefore only a short résumé is given 
here of the formulas needed in constructing the graphs of pressure 
expansion and pressure-summation. 


Tension Method 


The pressure-summation or radial-pressure differential across a 
cylinder wall 


2 o 
1 
P= 2/2 dr ( ) 


For small strains, as found in the elastic and plastic strain 
zones 


constant 


and in the elastic zone 


or? = constant 
V/3P b? 


For a partially plastic cylinder, where the initial yield has 
penetrated to radius r,, the pressure-summation for the elastic 
zone 


(5) 


and for a perfectly plastic zone between radii r, and rp 


9 


2 r 
P, = —~ log,“ (6) 
aro 
In the strain-hardening zone 
u(2r + u) = &2 + €)r? = constant (7) 
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and between radii r, and r = a equation (1) takes the form 


2 Tp + up 

Po =x 
V3 atue T+U 
This integration is carried out arithmetically by cutting the 
part of the cylinder wall in the strain-hardening range between 
radii r = a andr = r, into a suitable number of divisions. As- 
sume two consecutive radii r,, and r, and their strained values 

Tm + Un, andr, u,. The equivalent simple strains €, equal 


dir + u) (8) 


(9) 


The equivalent true stresses o corresponding to strains € are 
taken from a true tensile stress-strain curve. 

Between r,, and r, the equivalent stress is assumed to be 
constant and equal to 1/2 (¢,, + ¢,), and the radial stress gradi- 
ent across the subdivision per equation (6) 

_ + Tn 


Tr + Un 
= VJ/3 log. Tn + Un 


If the radial displacement is known or assumed for one radius, 
the displacement at all other radii may be calculated from 
equation (7) and the integration of equation (8) completed in 
accordance with 


(10) 


r=p 


1 Tr + Un 
v3 r=a 


Ta Un 


(11) 


A simplified method for computing the pressure-expansion 
curve is given in the following. A complete and exact pressure- 
expansion curve is necessary for certain types of work, for 
instance in case of autofrettage calculations, but in many cases 
the sole problem is to estimate the bursting strength of a cylinder. 
As already pointed out [6], the tension method leads to a simple 
formula for the bursting pressure 


+ % 


> = 2 
ge V3 log, K, (12) 
where 
ii 


The outside wall stress o, is taken from a true stress-strain. 
The corresponding strain is 


kK, 
= zx 


(14) 


Torsion Method 


The development is similar to and runs parallel with that 
of the tension method. The pressure summation 


2f2 dr 
r 


P, = 


(15) 


The expression for a constant cross-sectional area remains 
unchanged 


ur = er? = constant 


and in the elastic zone 


Tr? = constant (16) 
P b2 
(17) 
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u 
Tn 
4 
- 
: 
| 
u=e? = (2) 
| 
: 
| 
‘ 
(1-2 = 4 
(2) 
= 


When the initial yield has penetrated to radius r, 


In the strain-hardening zone 


u(2r + uw) = e(2 + €)r? = constant (7) 


and between radii r, and r = a equation (15) can be written thus 


—d(r + u) (19) 
at+ua r + 


From the known or assumed radial displacement at any radius, 
the displacement for all radii throughout the wall thickness can 
be computed from equation (7). The shear strain equivalent 
to a displacement of u at radius r 


= (20) 


From a torsion test a shear stress-strain curve may be derived 
using the method described by Nadai [5] and corresponding 
values of y and 7 can be found from the derived curve. 

Following the procedure outlined for the tension method, the 
radial-stress gradiant from radial increment can be written 
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Fig. 22 Radial stress /log r based on tension data 
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Tu + Un 


AP, = (Tm ») log, 
(T 


and the integration of equation (21) is carried out arithmetically 


r=p 
Tr Un 
> (7m + Tq) log, (22) 
Tn Un 


Strains at the Outer Surface 


At the outer surface for a closed-ended and wholly elastic 
cylinder the hoop strain and axial strain 


P 
= (2 - 2: 
= (1 — 2p) (24) 
E( K? — 


When the inner portion of the cylinder is plastic up to r, and 
the outer position is elastic 


2 
£ 
Ty 25) 


When all or part of the cylinder is in the strain-hardening zone, 
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Fig. 23 Radial stress/log r base: 
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(K = 1.25) 


Fig. 24 T-1 steel cylinders 


Fig. 25 HY-80 steel cylinders 


the strain can be found from the pressure-expansion curves, as 
described in the following. 


The Pressure-Expansion Curve at Large Strains 


The pressure-expansion curve is plotted on semilogarithmic 
paper following the method first proposed by Manning [7]. 
Table 5 shows the application of the method to a cylinder of 
diameter ratio K = 28 and having a bore strain vu = 1. The 
magnitude of strains and expanded radii throughout the entire 
wall from r = 1 to r = 28 is then found from equation (7). 

The equivalent strain for each value of r is given in Table 5, 
in column 4 for the tension method and in column 7 for the 
torsion method. The corresponding stress values o and 7 for 
the two methods are taken from the stress-strain curves. 

For each increment the sum of the stresses at the end points 
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is multiplied by the corresponding figure in column 6, and the 
summation is then carried out in accordance with equations (11) 
and (22). 

As an example, a complete computation for T-1 steel is given 
in Table 6. 

The pressure-expansion curves for the four steels tested are 
shown in Figs. 22 and 23. 

For any cylinder of given diameter ratio, the pressure-expansion 
curve is then obtained by considering a pair of ordinates a distance 
equivalent to log K apart and moving along the r-axis from right 
to left. At any position the strain at the outer surface will be 
equal to € = u/r in column 4, Table 5, corresponding to the r 
value of the right-hand ordinate, and the corresponding internal 
pressure is given by the vertical distance between intersections 
of the ordinates and the pressure-expansion curve. 
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Fig. 27 A-212 steel cylinders 
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DISCUSSION 
Irwin Berman? 


The overstrain tests on thick-walled cylinders are a welcome 
addition to the literature on this subject. In addition to the 
test results, the author has presented for comparison two cal- 
culated sets of results. One of these is obtained from tension 
tests by the use of the method proposed by the author in his 
reference [6]. The other calculated set of results is obtained from 
torsion tests by means of the method first proposed by Manning 
(reference {7]). This discussion will first consider the author’s 
use of each of these methods of calculation to determine the rela- 
tionship between the internal pressure and the tangential strain 
at the outer surface. 

It is regrettable that the author repeats the inconsistent and 
frequently incorrect determination of the basic tension method 
formulas which he also used in his reference [6]. This is par- 
ticularly unfortunate because it makes meaningful discussion 
rather difficult and tends to obscure the possible value of the 
method. Rather than make an attempt to unravel the author’s 
presentation, the basic equations of the author’s tension method 
will be derived in a manner consistent with the theory of elas- 
ticity. Thus the necessary consistent assumptions will be known 
and the results may be evaluated in a proper light. 

Before proceeding with the derivation, however, one very 
questionable step which is presented in the paper and which will 
later be discussed in more detail should be pointed out. If the 
author’s Equations (2) and (3) are combined with the definition 
of o which the author presented in his reference [6], the results 
may be expressed as 


Ee = (V/3/2)(o, — 6,) 


where € is the radial strain, o, the tangential stress, and o, the 
radial stress. This result which the author obtains in the elastic 
range may only partly be explained by the author’s incorrect use 
of «. He uses € = u/r whereas € is actually equal to du/dr and 
the tangential strain €, is equal to u/r. 

A consistent derivation of the author’s Equations (4), (5), (6), 
and (8) will now be presented. Where possible, the author’s 
notation will be used. The assumptions are listed: 


(a) The material is isotropic and homogeneous in its strength 
properties. 

(b) As the pressure is increased, an element of the material 
may exhibit three types of behavior: elastic, perfectly plastic, 
and strain-hardening, in the order listed. 

(c) The stresses are dependent upon the strains and not upon 
the history of loading. 

(d) The material is both elastically and plastically incom- 
pressible. 

(e) The axial strain vanishes throughout the cylinder. 

(f) Mises’ yield condition is applicable for initial yielding and 
for the perfectly plastic material. 

(yg) When the material has been strain-hardened, the shape of 
the vield surface does not change. 

(h) In the plastic region, the true stress a plotted against the 
conventional strain €, is identical to the relationship in the eyl- 
inder between the effective stress (1/3/2) — o,) and the 
tangential strain €, which is defined by u/r. 


The stresses and strains are independent of the axial and cir- 
cumferential co-ordinates and the principal directions of stress are 
radial, circumferential, and axial. The only equation of equi- 
librium which remains to be satisfied is 


o, — o, — r(do,/dr) = 0 (26) 


2 Research Division, 
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Incompressibility, the zero axial strain, and the use of the defor- 


mation expressions € = du/dr and e, = u/r result for small 


deformations in the equations 
é, = —e = C/r? (27) 


where C is a constant. 
Consider first a fully elastic cylinder. From elastic incom- 
pressibility and Hooke’s law for axial strain, the relationship 


20, =0,+ (28) 


may be obtained. By means of Equations (27) and (28) and 
Hooke’s law for circumferential strain the expression 


o, — 0, = 4EC/3r? (29) 


is obtained. The value of C is determined from Equations (26) 
and (29) and the use of the boundary conditions, at r = b, ¢, = 0 
and atr = a,o0, = —P. Thus 


C = (3/2E)[Pb?/(K? — 1)] (30) 


The zero pressure boundary condition at the outer surface, to- 

gether with Equations (26), (29), and (30), lead to the well-known 

Lamé equations for the radial and circumferential stresses. 
Define an equivalent stress o as 


20 = V/3(0, — (31) 
where in the elastic range o may be considered to be the stress in 
a tensile specimen for which the distortion strain energy of the 
tensile specimen and the cylinder are equal. By the use of Lamé’s 
equations and Equation (31), the author's Equation (4) is ob- 
tained. 

Consider now the partially plastic cylinder under internal pres- 
sure. Referred to the initial co-ordinates, the portion r, <r < 6b 
is in the elastic condition, the portion r, < r < r, is in the per- 
fectly plastic condition, and the portion a < r < r, is in the strain 
hardened condition. Equation (31) with modification of the 
meaning of o may be made applicable to all three regions. For 
the perfectly plastic portion, Mises’ yield condition reduces to 
Equation (31) if ¢ = a, is the yield stress of the tensile specimen. 
This is so because for the incompressible cylinder in plane strain, 
Equation (28) is valid for a perfectly plastic material. If it is also 
assumed that Equation (28) is valid in the strain hardened por- 
tion of the cylinder, a similar comparison between the stress in 
the tensile specimen and the stress in the cylinder reduces to 
Equation (31) in which o is the current yield stress of the tensile 
specimen. 

At the elastic-plastic boundary r = r,, Equations (29) and (31) 
are both valid. Thus the value of C for the partially plastic 
cylinder may be written as 


C= (V/3/2E)o,r,2 (32) 


By the use of Equations (26), (29), and (32) and the boundary 
conditions that at r = b, ¢, = Oandatr = r,, 0, = —P,, the 
author’s Equation (5) for the elastic region may obtained. Simi- 
larly, Equations (26) and (31) and the boundary conditions that 
atr =r,,o, = —P,andatr =r,, 0, = —(P, + P,) lead to the 
author’s Equation (6) for the perfectly plastic region. 

In the elastic and perfectly plastic regions, the radial deforma- 
tions are small in comparison with the radius so that the original 
radial co-ordinates are used. In the strain-hardened region, the 
deformations may be large so that the current radius is used in the 
integration. Also, Equation (27) should be replaced by the 
author’s Equation (7) which simply states that the cross-sectional 
area of the cylinder is constant. Equations (26) and (31) and 
the boundary conditions that at r = r, + u,, 0, = —(P, + P,) 
and atr = a+u,,0, = —(P, + P, + Pz) lead to the author’s 
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Equation (8). With the use of assumption (/) of this discussion 
and a true tensile stress-strain diagram, the author’s Equation (8) 
may be numerically integrated. 

The author in the section entitled “Strain at the Outer Surface” 
puts down three strain expressions, Equations (23), (24), and 
(25). These expressions do not follow from any previous work 
and are not used by the author in the tension method calculation. 
The author then proceeds to determine the pressure expansion 
curves. He uses the equations which pertain to the strain 
hardened region in all the regions of a partially plastic cylinder. 
Even in the case of a wholly elastic cylinder, he uses Equations 
(10) and (11) of his paper. It is important to know what effect 
this has upon the results. 

Fora wholly elastic cylinder, the tangential strain €, may be ob- 
tained by a combination of Equations (27) :nd (30). Atr = b 
the result may be written as 


(€:) = (1.5/2)(F — 1)j. (33) 


Consider now the author’s implicit assumption. Namely, that 
this discusser’s assumption (h) is true for the elastic cylinder. 
As pointed out at the beginning of the discussion, this means that 
the arbitrary assumption is made that for the wholly elastic 
cylinder 


In order to compare this result to Equation (33) substitute the 
values for o; and o, into Equation (34). The result may be 
written as 


(€,)rmep = (1.732/E)[(P/(K2 — 1)). (35) 


Thus it is shown that Equation (35) may be obtained from the 
correct Equation (33) by means of a multiplication of the right 
side of Equation (33) by 2/+/3. 

For comparison, the tangential strain at the outside surface of 
an elastic closed end cylinder in which the material is compressible 
may be written as 


= [(2 — w)/E][P/(K? — 1)}. (36) 


If steel with Poisson’s ratio 0.29 is used, the strain of Equation 
35) is within 2 per cent of that of Equation (36). Thus the good 
correlation between the author’s results and the test results in 
the elastic renge are due to the fact that the inconsistent factor 
of 2/+/3 by which the tangential strain is multiplied results in 
a value which is very close to that obtained from the most gen- 
eral elasticity expressions for a closed end cylinder. 

For a partially plastic cylinder, in plane strain, the stress and 
strain in any outer annulus of the tube is only dependent upon the 
radial pressure transmitted at the interface. Thus the discussion 
of the relationship between strain and pressure just given is ap- 
plicable to the elastic region of the partially plastic cylinder. 

It has been suggested (reference [3]) that the octahedral shear 
stress and strain should relate the tensile specimens to the 
cylinder. To accomplish this, assumption (h) of this discussion 
may be replaced by: 

h') In the plastic region, the relationship between the true 
octahedral shear stress and the octahedral shear strain is the same 
for the tensile specimen as it is for the cylinder. 

If the validity of this assumption is extended to the elastic 
region of the cylinder, the result may be written as 


Ee, = (3/4)(0, — G;) (37) 


By the substitution of the Lamé values for o, and o, into Equa- 
tion (37), Equation (33) is obtained. Thus the assumption of 
equivalent octahedral shear stress and strain in the tensile speci- 


men and in the cylinder is consistent with the theory of elasticity. 
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This, of course, does not mean that the relationship is necessarily 
valid or even preferable in the plastic range. It simply means 
that, if, as the author does, the single graphical method is used in 
the elastic as well as the plastic regimes, the result does not 
violate the basic expressions of the theory of elasticity. However, 
because of the assumption of elastic incompressibility, the indi- 
cated relationship between pressure and tangential strain at the 
outer surface would not correspond to that of an actual cylinder. 
It should also be noted that, for large strains, the use of the rela- 
tionship (1 + €,)(1 + €) = 1 is necessary in order to relate €, 
to the strain of the tensile test. This added difficulty does not 
arise in the author’s method since in his case €, is related directly 
to the strain of the tensile test. 

There is an obvious discrepancy in the results which were ob- 
tained from the torsion methods of calculation as shown for the 
T-1 material. In the graphical method (Table 6) the author 
again used a single expression (Equation 19) for both the elastic 
and plastic regions of the cylinder. In this case the use of this 
expression is consistent with the elastic analysis of the discusser. 
It leads to Equation (33) as the relationship between the outer 
tangential strain and the internal pressure. The author’s Tables 
5 and 6 in the elastic range are in agreement with Equation (33) 
in that, for a particular value of €, of the outer surface, the internal 
pressure as calculated by the torsion method is higher than the 
result for a closed end cylinder. For some unexplained reason, 
the graphs shown by the author in his Fig. 10 indicate that for 
identical values of €, the pressure as calculated by the torsion 
method is identical with that of the closed end cylinder. This 
discusser believes that the curves for the elastic region of the tor- 
sion method may have been obtained from a separate calculation 
by the use of Equations (23) and (25). This in effect means that 
Equation (33) was arbitrarily multiplied by 2(2 — y)/3. 

Before some concluding remarks the discusser would like to 
point out two important nomenclature errors which occur in this 
paper. The symbol € is frequently used instead of €,, i.e., Equa- 
tions (2) and (7) and the heading of Column 4, Table 5. Log r 
is used instead of “r plotted on a log seale,’’ i.e., the abscissa of 
Figs. 22 and 23. 

The tension method of calculation which is used by the author 
has merit for the pressure vessel designer because of its simplicity 
of calculation and its automatie ‘ 


” 


in the elastic re- 
gion to compensate for the assumption of elastic incompressibil- 
ity. Also, tension methods in general have the advantage of the 
greater availability of tensile stress-strain curves. 


‘correction 


Although torsion specimen stress-strain curves are more dif- 
ficult for the designer to obtain, their stress-strain relationship 
may resemble that of the overstrained cylinder more than do the 
tension specimen curves. The torsion graphical method as 
shown by the author, however, does not ‘‘correct’’ the pressure- 
strain curve for the assumption of incompressibility in the elastic 
component. 

The use of the methods of analysis considered in this paper 
must be limited to those cases in which they yield results within 
the required accuracy. For the materials tested, at large 
strains the change of the internal pressure with increase in the 
strain at the outside surface is a fairly small percentage of 
the total pressure. Thus these graphical methods, with cor- 
rection factors if necessary, could apply to the design deterraina- 
tion of the pressure at large strains. On the other hand, if the 
axial stress and strain are desired, it must be obtained by some 
other method. 

One should be very hesitant about the use of the methods of 
analysis considered in this paper in order to determine the magni- 
tude of the internal pressure which is needed for autofrettage and 
the subsequent stresses which are set up. The autofrettage 
strains are generally small so that most of the cylinder is either 
elastic or perfectly plastic. Thus both elastic compressibility and 
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a consistent application of the theory of elasticity may be very 
important in the determination of the subsequent stresses which 
are set up in the eylinder. 


B. Crossland’ 


As noted by the author all the tests which are the subject of 
this paper were carried out at the request of the Foster Wheeler 
Corporation of New York, in the Mechanical Engineering Labora- 
tories of the University of Bristol. My main anxiety when these 
tests had heen completed was that the results might be used to 
support one particular theory or another. It must be emphasized 
that the material provided was in the form of very roughly 
forged bar, and in the case of one of the steels, A-212, large cavi- 
ties were found throughout most of the bars. The results from 
numerous tension and torsion tests showed considerable and 
random variations in properties, much larger than would be 
normally found, or accepted, in practice. As requested, the tests 
were carried out on specimens in the as-machined condition, that 
is, without vacuum stress relieving, so in addition there was the 
possibility of residual stresses. 

No theoretical treatment could allow for these completely ran- 
dom variations in properties, nor for the unknown residual stress 
distribution in the virgin cylinders. Consequently, the results 
of these tests should only be used to demonstrate qualitative 
agreement with one theory or another, and they cannot be quoted 
to differentiate between, for instance, theoretical treatments 
based on tension and torsion. 

In his treatment of the analysis based on tension results, the 
author has departed from the analysis given in reference [3]; in the 
present paper he has assumed that the equivalent tensile strain 
in the evlinder is u/r. Now in a thick-walled cylinder it has 
been shown experimentally that there is no axial plastic flow, and 
consequently the plastic strain in the radial and tangential direc- 
tion must be equal and opposite if there is constancy of volume. 
However, in a tensile specimen there is plastic flow along the 
three principal axis, equal compressive plastic straining in the 
two transverse directions and a tensile plastic strain in the axial 
direction. How then is it possible to state that the equivalent 
tensile strength is equal to w/r unless the author advocates the 
use of the maximum strain criterion, which has been adequately 
disproved by numerous experimenters in the past sixty years? 
In reference [3] it was assumed that there was 4 unique relation- 
ship between the Mises-Hencky function of stresses and strain or 
between the octahedral stress and strain, whether this was ob- 
tained from a test on a torsion, tension, or thick-walled cylinder 
specimen. Using this assumption, for which there is some ex- 
perimental evidence (ef., J. L. M. Morrison, Proceedings, I. 
Mech. E., vol. 168) the analysis given in reference [3] is perfectly 
correct. 

Finally, I would like to comment once again on the reasons for 
preferring shear stress-strain data for this work. As has already 
been indicated, the state of stress in a thick-walled cylinder is one 
in which there is no axial stress deviation and the stress devia- 
tions in the radial and tangential directions are equal and op- 
posite. Consequently, the state of stress reduces down to a 
simple strain stress acting along 45-deg helices, with a super- 
imposed volumetric tensile stress equal to the axial stress. Now I 
have demonstrated experimentally (ef., B. Crossland, ‘The 
Effect. of Fluid Pressure on the Shear Properties of Metals,’’ Pro- 
ceedings, I. Mech. E., vol. 168) that hydrostatic pressure has no 
effect on the vield and plastic flow of ductile metals except to 
increase the strain to failure, so it is not unreasonable to assume 
that reasonable values of hydrostatic tension will likewise not ef- 


3 Formerly, Senior Lecturer in Mechanical Engineering at the 
University of Bristol; at present, Professor of Mechanical Engineering, 
The Queen’s University, Belfast, Northern Ireland. 
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‘fect the yield and plastic flow. If this is accepted, then the state 


of stress in a thick-walled cylinder is identical to simple shear as 
obtained in a torsion test. This then avoids all arguments about 
the correct criterion for yield and plastic flow, which simply 
does not enter the problem at all. The only point of dispute is 
then the question of complexity of the torsion test versus the 
tension test. From my experience there is no doubt that to 
obtain the necessary large strain data it is altogether more simple 
to carry out a torsion test. Admittedly having obtained a 
torque twist curve it is necessary to use Nadai’s construction to de- 
duce the shear stress-strain curve, but the validity of this con- 
struction has been verified on several occasions. In tension, 
however, even if the diameter and radius of curvature of the 
neck are measured as it develops the true stress has to be de- 
rived using an analysis suggested by Bridgman, which even 
though it may be correct has never been verified experimentally. 
Again, there is no autographic method of recording the diam- 
eter of the neck and the radius of curvature, and indeed unless 
a tension machine is used which has a very stiff load bar it is 
extremely difficult to obtain any information after the neck has 
just begun to form. 

In conclusion, I would say that for this particular problem, at 
least, there is little doubt that the shear-strain data are to be pre- 
ferred, and the steel companies should provide shear stress-strain 
data as a matter of routine rather than by special request. 


W. R. D. Manning’ 


In essence this paper follows up the previous work of this 
author in conjunction with B. Crossland and J. A. Bones® which 
was presented to the Society at Denver, Colo., in 1958. Four 
different steels were used, and tensile, torsion, and cylindrical 
bursting specimens were machined from the same material in 
each case. The true shear stress versus shear strain was extracted 
from both the tensile and torsion test results, and by their use the 
relationship between diametrical expansion and internal pressure 
for the cylindrical specimens was derived and compared with 
actual tests. The theory used to calculate the cylinder behavior 
from the shear stress data is a modification of that in the previous 
paper. 

On the whole the agreement is tolerably good although it is 
clear that none of the materials was as uniform as could be de- 
sired. On the other hand, the photographs of sections through 
some of the burst cylinders which are reproduced in the paper do 
not indicate any serious flaws in the case of T-1 and HY-80, al- 
though for the A-302 and A-212 the fissures appear to proceed 
along paths which are more nearly in the radial direction than 
would be expected from failure by shear, as admirably demon- 
strated by the examples shown for T-1 and HY-80. The A-302 
cylinder of K ratio 1.25 appears to have two cracks diametrically 
opposite one another. The photograph is not very clear, but if 
this is so it would suggest that one crack may be the result of some 
internal defect. 

The author notes that during the early stages of overstrain the 
cylinders did not expand uniformly in all directions and he sug- 
gests that anisotropy is the cause of this. It could have been ex- 
plained by the phenomenon of upper and lower yield stress, which 
causes the cylinder wall to yield in a series of discontinuous wedges, 
as Steele and Young® have shown for low carbon steels. Mr. 
Jorgensen, however, states categorically that his: specimens did 
not exhibit this phenomenon; nevertheless, the present reviewer 
cannot help wondering whether this may not still have been 

4 Imperial Chemical Industries, Ltd., Welwyn Garden City, Eng- 
land. 

5B. Crossland, 8. M. Jorgensen, and J. A. Bones, “The Strength 
of Thick-Walled Cylinders,” Trans. ASME, series B, JouRNAL oF 
ENGINEERING FOR INpusTRY, vol. 81, 1959, p. 95. 

6 M.C. Steele and John Young, Trans. ASME, vol. 74, 1952, p. 355. 
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present, at any rate in the cases of A-212 and A-302 where the 
stress-strain curves of Fig. 1 reveal considerable horizontal por- 
tions in the yield zone. One wonders whether his load-bar tech- 
nique was sensitive enough to resolve any possible fall in stress 
during this process. It would certainly appear from the photo- 
graphs of the burst sections that the two harder materials were 
much the cleaner, since they show perfect spiral shear fissures 
compared with the jagged more or less radial cracks of the others 
which, apart from this, one would expect to be much more duc- 
tile. 

The theory, which is described in some detail in the appendix, 
calls for little comment; essentially it depends on the fact that the 
stress system in a thick-walled cylinder ean be resolved at all 
times into a pure shear with a superimposed hydrostatic tension. 
It is then assumed that the latter plays no part in the process of 
overstrain, and the correlation between shear stress and shear 
strain is assumed to be the same in the cylinder as it is in either 
the tensile or the torsion tests. The author suggests that the 
hydrostatic tension should be disregarded in the plastic range, 
and by so doing he obtains a better agreement with his test results. 
In the opinion of this reviewer, however, this must be fortuitous 
since the removal of the hydrostatic tension component is in- 
consistent with the general assumption that plain sections re- 
main plain and perpendicular to the axis of the cylinder. If the 
fundamental theory does not fit these experimental facts so well, 
surely the explanation is more likely to be in the admitted lack of 
uniformity of the experimental material. 

The paper is presented in the clear style one has come to expect 
from this author, and his arguments are logically developed and 
easy to follow. It seems a pity therefore that this should be 
used to introduce some modifications to the more generally 
accepted theory when the experimental results are inevitably in 
doubt on account of the suspect material used. 


Author’s Closure 


It is noticeable that the discussions are mainly about the theory 
based on the tension method, and that much of the discussion has 
been supported by the strain correlations. 
Rather than being a clarification, this seems to have complicated 
the issue. 


introduction of 


The author agrees that a correlation must be estab- 
lished between the strain in the test specimen and some signifi- 
It will be seen that two of the dis- 
cussers, Messrs. Crossland and Berman, not only disagree with the 


cant strain in the evlinder. 


author but also with each other on this subject. 

This much discussed function is € = w/r, in which wu is the 
radial expansion as it is measured with a micrometer across the 
diameter of the evlinder. It is 
then postulated that the equivalent strain in the test specimen 
must equal e. 


The above function defines e. 


Further reasons for this assumption will be given 
in the following. 

On the subject. of strains, the author wishes to point out that 
equation (2) in the Appendix is a carry-over from reference ma- 
terial. It should not have been included in the paper, as it is not 
valid throughout the whole range of experiments, and it has not 
been used in developing the theory, nor has it been used in the 
interpretation of the test results. 

Dr. Berman’s discussion covers primarily one particular part 
Both in 
the subject paper and in the reference papers describing the deri- 
vation of the formulas, the author indicated that, in the elastic 
range, the cylinder material is considered to be compressible, 
and also, that the Lamé formulas for the component. stresses 
should be used to calculate the expansion of the cylinder. Only 
the cylinder parts in the plastic and strain-hardening ranges are 
considered incompressible. This admittedly results in a state of 
discontinuity at the interface between the elastic and plastic 


of the paper: The elastic range of the tension method. 


zones due to the abrupt change in Poisson’s ratio at this point. 
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However, the discontinuity in the theoretical curves from this 
source is generally within the deviations due to imperfections in 
the material, as shown by the pressure-expansion test results. 

It might be more consistent, as Dr. Berman maintains, to use 
a constant Poisson’s ratio = 0.5 throughout the whole stress- 
strain range, but it can hardly be called the correct approach for 
the elastic range. 

As mentioned earlier, € is defined by € = u/r and not by du/dr. 

The author has used two methods in calculating the pressure- 
expansion curve for a wholly elastic cylinder. One is equation 
(23) which is derived directly from Hooke’s law and the Lamé 
values for o, and o,. The other is derived by continuing the 
stepwise tabulation of Table 5, using equations (3) and (4). For 
all practical purposes the results from the two methods are identi- 
eal. 

On this basis, the author finds Dr. Berman’s equation (33) un- 
acceptable. A total of 32 test cylinders covering 6 different ma- 
terials and tested over a wide range of thicknesses and pressures 
up to 100,000 psi showed a remarkable accuracy when checked 
against the author’s method. If checked against Dr. Berman’s 
equation (33), the deviation would be around 15 per cent. 

For similar reasons Dr. Berman’s equation (37) is unacceptable, 
as this equation apparently is based on assuming the cylinder 
material incompressible. 

It is evident that the difference between Dr. Berman’s and the 
author’s method can be expressed by the ratio of 2- for incom- 
pressible and compressible material. Dr. Berman has found 
this ratio to be close to 2/4/3, which is the same as the ratio of the 
von Mises-Hencky functions for test specimen and cylinder. 
This is entirely fortuitous, as there is no connection whatsoever 
between Poisson’s ratio and the von Mises-Hencky ratio. 2//3 
equals 1.15, which accounts for the 15 per cent deviation of Dr. 
Berman’s equation (33), as mentioned before. The author 
therefore contends that the factor 2/+/3 is not inconsistent but is 
in the equation for pressure-summation with good, theoretical 


“backing. 


There is no discrepancy in the results obtained from the tor- 
sion method. The calculations run parallel with those for the 
tension method. The integration is carried out arithmetically 
in the strain-hardening range, and the expansion in the elastic 
range is found from equation (23). 

With regard to using the methods for autofrettage, it should be 
pointed out that autofrettage is most likely to be used for vessels 
operating at. very high pressures. The material used for such 
vessels will be high tensile, heat-treated steels similar to T-1 
and HY-80 steels, which do not have a perfectly plastic range in 
their stress-strain curves. 

The methods generally used at present for autofrettage calcula- 
The 
methods proposed in this paper permit calculations for the 
strain-hardening ranges also and should therefore prove far supe- 
rior to the older methods. 

The author agrees with Professor Crossland that the results of 
these tests cannot be used to demonstrate that of the two methods 
described, based on either tension or torsion tests, one is superior 


tions cover only the elastic and perfectly plastic ranges. 


to the other, and no such claim has been made in the paper. 

The problem of whether the equivalent strain in the tensile 
test. specimen should equal u/r or (2/+/3):(u/r), has been the 
subject of extensive discussions between Professor Crossland and 
the author, but no agreement has been reached. The author is 
of the opinion that the strain correlation equation (32) in paper 
reference [3], is based on assumptions that seem very questiona- 
ble. In that correlation, the strains in the tensile specimen are 
separated into a principal strain and two lateral strains, whereas 
the cylinder strains are total strains. This strain correlation 
therefore is not symmetrical with the stress correlation. 

The argument is limited to stresses and strains outside the 
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elastic range. The author suggests that the disputed function is 
supported by the following argument: 


1 At a radius r in the cylinder, the tangential and radial 
strains are of opposite sign, but equal in size to u/r. 

2 There is a unique relationship between u/r and stresses 
o, and o, throughout the cylinder. 

3 We now form the correlations 


 (2/r/3)0 
u/r € 


4 It is postulated that € = u/r. If this were not so, then at 
some radius r the strain would be € = ku/r; we would then have: 


o, + 0, — 
But that would conflict with the basic stress correlation 


0, +4, = (2//3)o 


The value of any argument whether it supports one theory 
or another is, at best, of dubious value considering our limited 
knowledge of materials in the plastic and strain-hardening state. 
Until more data are available any theory must be evaluated by 
comparison with experimental data. 

Those who Professor Crossland’s work will 
know he is a firm believer in the superiority of the torsion test 
method and he always presents his ease well. 


have followed 


The author will- 
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ingly concedes that this method has an advantage in simplicity. 
However, on the basis of experimental work so far, both methods 
seem equal in their close agreement with the experimental results. 

Mr. Manning is correct in assuming that, in the case of A-212 
and A-302, both of these steels exhibited an upper and lower 
yield stress, which certainly would account for the uneven ex- 
pansion of the cylinders from these materials. 

It was unexpected that uneven expansion also occurred in the 
‘ases of HY-80 and T-1, as these steels do not have an upper and 
lower yield stress. 

Mr. Manning also comn-ents on the proposed suggestion to 
consider the hydrostatic component negligible. However, this 
idea was not developed and is not used in the theory, and the 
misleading sentence was deleted. 

The change in theory is not due to any modification of the 
generally accepted theory. Acknowledging the sharp difference 
of opinion, the change has been described as an ‘“‘empirical cor- 
rection factor.’ It is, however, the author’s opinion that the 
present version is the more correct. 

In conclusion, the author should like to point out that be- 
havior of materials when overstrained is a comparatively new 
field, and that every new investigation generally has resulted in a 
new theory. With the hope that these tests will contribute to- 
ward better understanding of the problem, the author wishes to 
thank the discussers for their very generous and valuable con- 
tributions. 
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§ ALMOST any investigation of vibrations in drill 
strings, it is necessary to consider their natural frequencies. 
The calculation of these frequencies may not be intrinsically 
difficult, depending on the assumptions made, but for any but 
the simplest drill-string configurations it is generally very tedi- 
ous. In this work, a quick method of making these calculations 
is demonstrated, and the effects of some parameters on the natu- 
ral frequencies are considered. 

The drill string considered consists of drill collars of one or 
more different sizes and of drill pipe. The string rests on the 
bit at the bottom of the hole and hangs from the traveling block 
at the top. The boundary conditions that result at the ends of 
the string are important in the calculations, but they are not 
at all well understood at present. They will be discussed later. 
The effects of any driving or damping of the vibrations that 
may occur in practice are neglected. Furthermore, only lon- 
gitudinal and torsional vibrations of the string are considered. 
Lateral motions of the string (due to bending, buckling, whirling, 
whipping, and so on) are neglected, and it is assumed that the 
torsional and longitudinal motions considered are independent. 


Vibration of Stepped Drill Strings 


The case of longitudinal vibration in a drill string as shown in 
Fig. 1 will be considered in terms of 


— = the longitudinal displacement 

x =a length co-ordin: . of the string 

t = time 

w = angular frequency 

c = speed of propagation of longitudinal (sound) waves in the 
material of the drill string 

F = axial force in the string 

= Young’s modulus 


= cross-sectional area of the string 


m 


The symbols F and € refer only to the unsteady quantities as- 
sociated with vibration. 

Longitudinal and torsional vibrations in drill strings are ex- 
actly analogous and, if only the symbols are changed, the same 
equations govern both. Hence the equations and results that 
follow apply immediately to torsional vibrations if the following 
substitutions are made: 


£ becomes the angular displacement 


Contributed by the Petroleum Division and presented at the 
Petroleum Mechanical Engineering Conference, Houston, Texas, 
September 20-23, 1959, of THe AMERICAN Society oF MECHANICAL 
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An Analytical Study of Drill-String Vibration 


Longitudinal and torsional vibrations in a drill string consisting of drill pipe and 
collars, and the boundary conditions at the ends of the string, are discussed. 

A trial-and-error method of solution for the natural frequenctes 1s demonstrated. 
Since it is very lengthy, certain useful charts are then derived, and a graphical method of 
solution for the natural frequencies is developed. By use of this method the natural fre- 
quencies can be quickly and easily obtained. The influences of various parameters on 
the natural frequencies are considered. 


c becomes the speed of propagation of torsional (shear) waves, Cg 
A becomes the polar moment of inertia of the string, J 

F becomes the torque in the string 

E becomes the shear modulus 


Also, the boundary conditions at the ends of the string will be 
different in the two cases. 
The equation governing the longitudinal motion is 
= — (1) 


and the general solution to this equation is 


&(z, = | A sin r+ Beos r)(Dsinwt + E cos at), 


(2) 
where A, B, D, and E are arbitrary constants, and wa/c is a di- 
mensionless length parameter describing the 2 co-ordinate in 


As 


7 EROSS SECTIONAL AREA A, 
™~POLAR MOMENT OF INERTIA I, 


Fig. 1 Stepped drill string nomenclature 
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terms of the wave length of the vibration. 
of as a phase angle and is also referred to as Q. 

The natural frequencies of the system are found by finding 
values of w that fit the required boundary conditions. zx is de- 
fined as zero at the bottom of the string and increasing upward. 
If we assume that the bottom of the string is fixed, ie, & = 
0 at x = 0, and choose the time scale arbitrarily so that & = 
0 for t = 0, the solution in the lowest section of drill collars can be 
written 


It may be thought 


&, = C, sin = x sin wt, (3) 
where C, is the amplitude of the vibration. 

Consider the motion in the next section of the stepped string, 
which begins at z = 22. The general solution for the motion in 
that part of the system can be written 


cs) sin (wt + 


where €;. and €' are constants chosen to suit the boundary condi- 
tions at the join. These boundary conditions will be taken as 


€12'), (4) 


= 22) (5) 


and 
F (a2) = F222), (6) 


i.e., the string is continuous at the join and there is no net force 
there. 

The boundary conditions in equations (3) and (4), and the 
result F = EA d&/dz2, can be used to show that 


tan ( Lie + «:) 
Ay 


A, 
tan — 22 


and 


sin Tie t 


Thus at the change of section A; to A» at ty, there is both a 
phase discontinuity, €2, and a change in amplitude of the wave- 
form (C,/C; # 1, in general), though of course the amplitude at 21, 
is continuous as required by equation (5). Fig. 2 shows the 
nature of the change at the discontinuity. 

Now consider the application of these results to the actual 
calculation of the natural frequencies of a drill string. The 
boundary condition & = 0 has already been assumed at x = 0. 
Values of A2/A; are known at each change of section of the drill 
string. Ifthe boundary condition at the top of the pipe is known 
(or at least assumed), solutions can be made. 
numerically by trial and error. Choosing some frequency w and 
starting at one end of the drill string, we solve transcendental 
equations like (7) at each change of section, and the boundary 
condition that is calculated for the other end of the pipe is then 
compared with the required condition. Further trials are made 
until a sufficiently accurate solution is found. The procedure 
can then be repeated for other modes. In order to find the 
relative magnitude of the fluctuating stresses in the string, it is 
necessary to solve equation (8) at each change of section for each 
mode of interest 
If many solutions are to be made, the tedious numerical work 


This can be done 
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involved in this method of solution can be largely eliminated by 
the use of graphical solutions to equations (7) and (8). 
It will be seen from equation (7) that 


A c 


C; C; (4 c 
Values of €;2 and C2/C, in terms of these parameters are shown in 
Figs. 3 and 4. If J2/I, is used instead of A2/A; and w2/ce in- 
stead of wr/c, then € represents the phase jump and C2/C; the 
amplitude ratio in the torsional case. Instead of plotting over 
some range of A,/A,, it is convenient to take values correspond- 
ing to specific tubular-goods sizes of interest. Thus in Figs. 3 
and 4 the curves are given such identifications as 113/,-in. DC 
to 7%/,-in. DC (longitudinal), and so forth, and the numbers 
shown are the values of A2/A; or J2/J; as the case may be. For 
other combinations of tubular goods, the curves can readily be 
calculated or they can be plotted by interpolation when we know 
the required values of A2/A,; or J:/J;. A sample calculation 
demonstrating the method and the use of Figs. 3 and 4 is shown 
later. 

The drill-collar string is generally much shorter than the drill 
pipe. Also, the boundary conditions at the bottom of the string 
seem to be less uncertain than those at the top. Thus calculat- 
ing from the bottom of the string upward instead of from the 
top downward is likely to result in less error in w2j2/c and w2o:/c. 
It can be seen in Fig. 3 that quite a small change in wz,2/e can in 
some circumstances give a relatively large change in €. Hence it 
is generally more accurate to start the calculations at the bot- 
tom of the string. 


€i2 = 


and hence 
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Fig. 2. Predicted second mode of longitudinal vibration 
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Fig. 3 Phase-angle discontinuity, «°, versus phase angle, {2°, at change 
of section (K = Ao/A; or I2/ii, etc) 


Sample Calculation of Natural Frequency and Relative Displacements 
and Stresses. ‘To illustrate the method previously described, con- 


sider longitudinal vibrations in a drill string consisting of 


95 ft 
321 ft 
3175 ft 
3591 ft 


11%/,-in.-DC 
73/.-in. DC 
6-in. DP 
Total 
*Note that this 6-in. DP is not regular drill pipe. It has the 
dimensions 6.000 in. OD, 5.424 in. ID, 18.4 lb/ft without tool 
joints, and 19.64 lb/ft including tool joints. 


In calculating the natural frequency of the second longitudinal 
mode, we assume both the bottom and top ends of the string 
to be fixed. The relative amplitudes of the displacements and 
the stresses in various parts of the string will be found. 

For longitudinal vibration, the angular phase change along a 
section of steel pipe of Jength L is 


- L = 2.195 deg, 


where n is the frequency in eps; i.e., the phase changes in the 
different parts of the string are as follows: 


113/,-in. DC 
73/,-in. DC 
6-in. DP 


2.09 n deg 
7.05 n 
69.6 


A table as will be shown is then constructed. 
represents a trial solution. 


Each column 
The table is essentially a table of 
phase angles (or changes), given to the nearest 1/. deg. The 
lines are: 
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Fig. 4 Amplitude ratio C./C), etc, versus phase angle, ‘)°, at change of 
section (K = A:/A; or ete) 


1 The frequency in eps that is being tried. 
2 The bottom boundary condition, 22 = 0 in this case. 
3. The phase change in the 114/,-in. DC, AQ = 2.09 x row 1. 
{4 The phase angle at the bottom of the 11%/,in. DC, row 
4 = row 2 + row 3 [in this case where the bottom end is taken 
as fixed, i.e., 2 = 0, rows 2 and 4 are trivial, but in general they 
must be included]. 
5 The phase jump from the 11%/-in. DC to the 73/,-in. DC, 
€2; this result is read from Fig. 3. 
6 The phase angle at the bottom of the 7/,-in. DC, row 6 = 
row 4 + row 5. 
Proceeding similarly, 
= 7.05 X row 1, 
= row 6 + row 7, 
is read from Fig. 3, 
= row 8 + row 9, 
= 69.6 & row 1, and 
= row 10 + row LI, the resulting boundary condition at 
the top of the string. In this case row 12 should be 360 deg in the 
desired solution (from the original assumption). 
5.0 5.05 5.07 
0 0 
10!/., 
~ 


bottom 
11/,-in. DC 
10'/, top 


bottom 
73/,-in. DC 
top 


5 

3851/2 
401/, 

—33 


71/, 


—26}/2 


bottom 
6-in. DP 


top 


279 


284!/5 


351"/> 
"359 
That is, the natural frequency of the second longitudinal mode is 
5.07 eps. ‘Thus an accuracy of better than 1/2 per cent has been 
rapidly attained. In practice, the calculation is complicated 
by the fact that the required top boundary condition is a func- 
tion of the frequency. <A graphical method of solution that fol- 
lows takes this into account and also eliminates the need for suc- 
cessive approximations, except for very accurate calculation. 
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Suppose the displacement amplitudes in the three parts of the 
string are C,, C2, and C;, as defined in equation (3), and so forth. 
Then from Fig. 4 


and 


giving 


C3/C, = 11.35. 

Now the displacement amplitude at any point is C sin Q, and the 
amplitude of the resulting stress component is proportional to 
CcosQ. The following table shows the relative displacement and 
stress amplitudes in the string. If the actual value of stress or 
displacement were known at any one point, a table of actual 
values could be made. 


113/,-in. 73/,-in. 
DC 


DC 6-in. DP 
Relative Bottom 0 0. 2002 0.148 
displacement Top 0.182 0.149 0 
Max. value 0.182 0.200 11.35 
Relative Bottom 1.00 2.29 8.75 
stress Top 0.980 1.72 11.35 
Max. value 1.00 2.29 11.35 


* This apparent inconsistency (9 per cent) is a result of working 
to the nearest '/, deg with small angles and is not considered im- 
portant for the present purposes. 


The relative stress amplitudes above are those due to longi- 
tudinal vibration only. The stresses due to the weight of the 
string and torsional vibrations (and indeed to bending and any 
other stress-producing mechanisms) must also be considered in 


order to estimate the actual stress at any point. 


Boundary Conditions at the Ends of the String 


For accurate calculation of the natural frequencies, the bound- 
ary conditions at the ends of the drill string must be known with 
precision. Unfortunately, this is not the case. 

Consider first the conditions at the top of the string. If 
measurements of both strain and displacement were made at 
the top of the string for both the longitudinal and torsional os- 
cillations, it should be possible, by comparing the amplitudes and 
phasing of the strain and displacement, to estimate the boundary 
conditions. In a series of tests run recently and deseribed in an 
accompanying paper,! such measurements were made, but they 
did not prove satisfactory for the calculation of the boundary 
conditions. 

As an alternative it is possible to estimate the boundary condi- 
tion at the top of the string for the longitudinal case by con- 
sideration of the mechanical system there. The drill pipe hangs 
from the traveling block, which is effectively a concentrated 
mass. The block then hangs from the derrick. The drilling 
lines and derrick are treated as a spring. From the data availa- 
ble on the components, an estimate can be made for the stiff- 
ness of the spring and the behavior of the spring-mass system. 
The spring constant of the system can also be measured directly 
by setting the pipe in the slips and taking a strain on it with 
the elevators. The length of line reeled onto the drum for a 
known change of load (without moving the pipe in the slips) is 
measured, 

Perhaps a more directly meaningful experiment would be to 
hang the drill string on the hook and then suddenly change the 
load on the hook. This might be done, for example, by at- 


1], Finnie and J. J. Bailey, ‘‘An Experimental Study of Drill-String 
Vibration,” printed in this issue, pp. 129-135. 
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taching a rope to the hook and to the rig structure and lifting 
until the rope parted. By observing the resulting free oscilla- 
tions of the system we could determine the boundary condition. 
Such an experiment might also be adapted for the investigation of 
torsional vibrations, 

Although, for longitudinal] motions, a spring-mass system seems 
avery plausible analog for the equipment at the top of the string, 
it is a much less satisfactory conception for the torsional motions. 
Perhaps the kelly bushing might be considered as a concentrated 
moment of inertia and the chain drive as providing a spring 
element. At least it appears likely that the natural frequency of 
the system driving the drill pipe is low and that it is felt by the 
drill pipe as a more-or-less fixed end. There seems to be some 
experimental evidence for this, but at best the torsional end con- 
ditions are rather uncertain. 

At the bottom of the string, the conditions are perhaps some- 
what more certain. It seems unlikely that the bit operates 
under conditions such that it moves up and down and jumps 
clear from the hole bottom. Assuming it does not jump, and 
neglecting the up-and-down motion of the bit due to the rolling 
of the teeth, we can treat the bottom of the string as a more-or- 
less fixed end, since the amplitude of 0&/dt at z = 0 cannot then 
exceed the penetration rate. This argument would not hold at 
frequencies as high as the tooth frequencies. These also seems 
to be the possibility of driving at lower frequencies, correspond- 
ing to the cone frequency or to the difference between tooth fre- 
quencies on different cones, for example. Nevertheless, the bot- 
tom of the string has been assumed to act as a fixed end for the 
present purposes. 

Tests have shown that the torque of a bit is not very much 
affected by the rotary speed over normal speed ranges. Thus 
for torsional oscillations the bottom end of the string can be 
treated as more or less free. 

fecapitulating, then, values at the bottom of the hole for the 
dimensionless length parameter Q will be taken as 0 deg for the 
longitudinal case and 90 deg for the torsional case. At the top of 
the string in the longitudinal case the parameter may take values 
from 0 to 180 deg within the frequency range of interest. In 
the torsional case 2 is thought to be fairly small within the 
range of interest. 


Calculation of Boundary Conditions at the Top of the Drill 
String 

Consider the idealized system shown in Fig. 5. The derrick 
and drilling line are represented by the light spring A, with'spring 
constant k. The traveling block is represented by mass M, B in 
the figure. The rod C represents the drill string. 

The motion of the top end of the rod, z = L, can be written as 


oL , 
= asin — sin wt, 
c 


The tensile load in the rod is 


and if we take into account the equilibrium of the mass, 


M 
atz = L, 
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: 
2 
Cy C, = 2.30 
C3/Cz = 4.95, 
: 
* 
An —, 
ox 
AG 
which gives the result 
wL _,| ¢ {Mw k 
c AE \ g 


In the torsional case, the model may have little validity. The 
moment of inertia of the drilling table is the analog of the mass B. 
The chain drive from the compound might act somewhat as the 
spring A. It seems likely that the effect of the mass would pre- 
dominate over the range of interest, giving a more-or-less fixed 
end at the top of the string (i.e., #/cg small). 
perimental evidence for this. 

S ple Calculati 


rhere is some ex- 


Boundary conditions at the top of the 
drill string previously described. 

For the 6-in. DP, A = 6.22 in.?, giving c/AE = 8.80 x 107% 
ft/lb sec. 

The total weight of the hook traveling block and swivel was 
23,000 Ib. An estimated correction was applied for the fact 
that the ‘spring’ in the system was not “‘light,’’ giving M/g = 
8.82 lb-sec?/ft. 

The effective spring constant for the derrick and drilling line 
system was estimated from the dimensions of the components as 
Kes. = 8.00 X 105 lb/ft. Later, an experimental value for k was 
found by measuring the length of drilling line reeled in to effect a 
known change of load on the hook. This procedure gave kexp = 
3.29 & 105 lb/ft. 

Putting these values into the result of the previous section 
gives the required top boundary conditions 


(0.503 n — 4.62/n) 


(0.503 n — 1.122/n), 
where n is the frequency in cps. The corresponding natural fre- 
quencies of the system in the absence of the drill string are 3.04 
and 4.73 eps, respectively. 


Graphical Solution for the Natural Frequencies 


The complicated form and the general uncertainty of the 
boundary conditions make numerical solutions for the natural 


bJ 
Fig. 5 Idealized model of equipment at top end of string. (A) Spring 


representing derrick and drilling lines. (B) Mass representing traveling 
block, swivel, and hook. (C) Rod representing drill string. 
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frequencies rather unsatisfactory. It is laborious to find the solu- 
tions numerically and even then, unless a large number of 
solutions are made, it is hard to see the effects of changes of the 
various uncertain factors. A graphical method of 
can be used which obviates these difficulties. 

The graphical method for the longitudinal case is shown in 
Fig. 6. Starting at the bottom end of the string, we assume a 
boundary condition corresponding to a fixed end. 


solution 


Now the re- 
sulting phase angleQ at the top end of the string is plotted against 


4 5 6 
Nn cps 
Fig. 6 Top-end condition, ‘)°, versus frequency, n cps, for longitudinal 
vibration 


frequency. For convenience, values 0 < 2 < 180 deg are plot- 


ted, giving a series of diagonal lines. 


Points on these lines are 
calculated by use of the numerical method shown earlier. It 
has been shown that the boundary condition required at the top 
of the string is frequency-dependent, 


Q = cot™ (“we 
AE g 


This curve is plotted on the same diagram, and where it inter- 
sects the diagonal lines previously drawn, the two parts of the 
system are so matched that there are no net forces at the top of 
the string; that is, the frequency indicated is a natural frequency 
of the combined system. In this way, good accuracy (within 
1/10 eps in the range 0 to 10 cps shown, for example) is obtained 
quickly and without successive approximations. 

Furthermore, the effects of some changes in the system can be 
seen from Fig. 6. If the string were of constant cross section, 
the diagonal lines would be straight with slope 2.195 x 107? L 
deg/cps. The deviation is due to the stepped nature of the 
Also, it can be seen from Fig. 3.that this deviation changes 
most rapidly with frequency when the value of 2 at any change 
of section is near 90 deg, 270 deg, and so forth. 

The diagonal curves show at n = 0 eps the value of the bottom 
boundary condition chosen. 


string. 


Errors in the choice of this bound- 
ary condition will in effect move all these lines parallel to the n 
axis. Thus the effect of errors can be roughly estimated. In the 
case shown, an error of 20 deg, say, in the choice of the bottom 
end as fixed would give an error of about 0.25 eps in the calculated 
natural frequencies, and all the natural frequencies would be 
affected. This figure is only an approximation, for changes in 
the bottom boundary condition will, in fact, somewhat change the 
shape of the diagonal curves. Nevertheless, it may be useful 
for small changes, and it does give an idea of the sort of errors 
due to the uncertainty of the boundary condition. 

The natural frequencies are also dependent on the boundary 
condition at the top of the string, which is represented by an S- 
shaped curve. At low frequencies, the spring predominates and 
Q approaches 180 deg, whereas at high frequencies the mass pre- 
dominates and 2 approaches 0 deg. For 2 = 90 deg, the fre- 
quency is the natural frequency of the spring-mass system, 
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The slope of this line atQ = 90 deg is proportional to Mc/EAg, 
i.e., larger M or smaller A makes for a more rapid transition of 
the top boundary condition from large to small values of Q for 
constant values of M/k and vice versa. Two curves for the top 
boundary condition are shown in Fig. 6 corresponding to two 
values of the effective spring constant k, differing by a factor of 
2.43 (the values are the experimental and estimated values men- 
tioned previously). It will be seen that the change in k has 
little effect on the natural frequencies except near the natural fre- 
quency of the system supporting the drill pipe (derrick, drilling 
lines, traveling block, and so forth). This is a convenient re- 
sult, for in practice the value of k will vary appreciably (by a 
factor of 1.5, say) as the traveling block moves between its ex- 
treme positions. 


DISCUSSION 
J. E. Goldberg? and J. L. Bogdanoff* 


It is common knowledge that the dynamical behavior of a drill 
pipe in rotary drilling is extremely complicated, being determined 
by the method of drive, the forces acting on the bit, and by the 
forces arising from contact of the pipe with the sides of the hole 
or with pieces of rock between the pipe and the sides of the hole. 
In particular, the time behavior of the forces at the bit and on 
the side of the pipe cannot be described in any deterministic sense 
because of their random or chance nature. These points suggest 
that a statistical description of drill pipe dynamics might be 
fruitful. 

We presented a paper entitled ““A New Analytical Approach to 
Drill Pipe Breakage” at the September, 1958, Denver meeting 
of the ASME Petroleum Mechanical Engineering Conference. 
The central purpose of the paper was to exhibit the advantages of 
a statistical approach to the very complex problem of drill pipe 
dynamics and to recommend an experimental program to secure 
necessary statistical information or input data with which this 
new approach could be put on a quantitative basis. A certain 
aspect of drill pipe dynamics is discussed by the present authors 
from a deterministic point of view. It is therefore of considerable 
interest to contrast the two methods of approach. It is necessary, 
however, to first describe the system considered by us inasmuch 
as our paper is not readily available. 


COUPLING 
PRIME DRILL PIPE BIT 


OVER tL Q x 


° 


Fig. 7 


The mode] considered in our paper is shown in Fig. 7. For 
purposes of preliminary discussion, only longitudinal displace- 
ment u(z, t) and torsional displacement O(2, t) were allowed. 
Gravity was neglected, viscous damping was assumed, and the 
pipe was taken to be uniform in cross-section;* this leads to the 
equations of motion® 

2 Professor of Civil Engineering, Purdue University, Consultant, 
Midwest Applied Sciences Corporation, Lafayette, Ind. 

3 Professor of Engineering Sciences, Purdue University, Consultant, 


Midwest Applied Sciences Corporation, Lafayette, Ind. Mem. 
ASME. 

4Nonuniform pipe complicates the computations but not the 
basic methods. 

6 Notation is that of our paper. 


Journal of Engineering for Industry 


Ou 
—- + 2b =a,’ 
ot? 4 ot dr? 
(9) 
00 
— + 2bp — = a7? 
ot? ol Ox? 


The boundary conditions at the left end, i.e., at the top of 
the hole, depend upon the method of drive, method of support. 


ete., and in each problem may be estimated fairly easily. We 
took for purposes of illustration only 
ou 
AE = k,u(0, ¢) | 
Or 
(10) 
00 | 
— = k,O(0, t) 
or zr=0 


which represent the conditions for an elastic connection. 

We assumed that a random axial force F(l, t) and a random 
couple T(J, t) to be second-order mean square continuous, weakly 
stationary random functions with zero means 


E{F(l,i)} =0, =0 (11) 
and covariance 
— t) = E{T(l, t)} (12) 
The boundary conditions at the lower end were therefore 
! \ 
AE—| | 
or r= ! 
> (13) 
00 | 
=T(l,t) | 
Or 


It might seem at first sight that w(x, t) and O(z, t) do not de- 
pend upon one another. This would appear to follow from the 
fact that the equations ot motion and boundary conditions for the 
two displacements are separate. However, F(l, t) and T(l, t) 
are related statistically through (12). Hence u(z ,t) and O(z, t) 
are not independent of one another in our analysis. 

Our views on the boundary conditions at the top of the hole 
coincide with the authors’; i.e., they may be estimated with fair 
accuracy for each type of drive. While these conditions may be 
complex in certain circumstances, they can be handled without 
undue difficulty. 

We have assumed an axial force and torque (random, however) 
acting at the bottom of the pipe. This assumption leads to 
natural frequencies, when damping is neglected, which correspond 
to a freeend. The authors, on the other hand, assume that the 
string is fixed longitudinally and free torsionally at the bottom. 
The assumption that longitudinal displacement is zero implies the 
presence ot a rigid constraint and neglects elasticity of the rock, 
of the bit, and of the possibility of the bit bouncing up from the 
rock which it is supposed to be cutting. If one really wants an 
accurate estimate of natural frequencies, such factors should not 
be ignored. 

The natural frequencies in undamped linear systems provide 
only partial information about the dynamical behavior of such 
systems. If, for example, there are simple harmonic components 
in the exciting force system which coincide with one or more of 
the natural frequencies and do not act at nodes of the correspond- 
ing modes, the phenomenon of resonance occurs. Amplitudes 
may not be estimated, however, unless damping is assumed. 
Moreover, if the actual damping present is large, undamped 
natural frequencies are at best only rough guides to frequencies of 
oscillation. As is well known, a complete description of dynamical 
behavior of a damped linear system is contained in the frequency 
response (responses if there is more than one point of excitation). 
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The frequency response in longitudinal motion which we ob- 
tained is 


cosh + sinh 


14) 
WwW) = —— (1 
A AE y* . | 
-sinh y + 7 cosh Y 
a 
where 
k,l 
a 
AE 
(15) 
il / 
Y= — Vw? — 
a4 
'G,(2, iw)| becomes large when w is such that 
| 
sinh y + cosh ¥)} (16) 
a 
is small compared to the absolute value of the numerator. For 
b, = 0, ie., for the undamped case, this occurs when w satisfies 
the equation 
wl k,l 1 ae 
tun — = = (17) 
4 AE ol 
a4 


and these w are the undamped longitudinal natural frequencies of 
our model; they do not correspond to those obtained by the 
authors. As by increases from zero, the values of w satisfying (17) 
depart further and further from the values of w which make 
1G4(a, iw)| large. Because of the possibility of rubbing between 
the pipe and wall, the presence of internal damping, the presence 
of slipping in tool joints, ete., it is reasonable to expect sub- 
stantial damping in drill strings. In addition, the length of a 
string changes as cutting progresses. Hence undamped natural 
frequencies begin to assume a somewhat academic aspect even if 
a deterministic point of view is adopted in studying drill pipe 
dynamics. 

These remarks also imply that a visual search for a set of 
natural frequencies based on the assumption of no damping and 
definite boundary conditions in records such as the authors 
present, which appear to be random, is hardly adequate. What 
is needed is a power spectral density analysis of these records. 
Sharp peaks in these densities would indicate either dominant 
exciting frequencies or natural frequencies, if any of the latter 
are present. 

Many of the quantities of interest in our paper are of the form 


|H(x, iw)|2P(1, w)dw 
0 


(18) 


H,(2, iw) = AE 
Ox \ (19) 


P,(l, w) = power spectral density of axial bit force 
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The peaks of |H4(2, iw)|? as a function of w determine the re- 
gions of resonance. The frequencies at which these peaks occur 
are determined by the boundary conditions, physical constants, 
and damping. Their magnitudes, if large, are primarily controlled 
by damping. If P4(l, w) is a slowly varying function of w, the 
exact positions of these peaks will not appreciably influence the 
value of (18), provided the magnitudes of the peaks are about the 
same. If the peaks are not large, the same remark is valid even if 
P,(l, @) varies somewhat more rapidly as a function of w. Hence 
(18) is not unduly influenced by changes in boundary conditions 
and damping, if the latter is reasonably large. Thus our analysis 
is much less sensitive on the whole to doubtful or uncertain 
boundary conditions than a deterministic analysis. 

Since the very essence of a statistical analysis is the accommoda- 
tion of uncertainty (in loads in our analysis, but uncertainty in 
boundary conditions and damping could also be included), its 
range of application in highly complex problems involving ran- 
dom loads, etc., is substantial and the advantages it offers in such 
problems is worthy of consideration. 


Authors’ Closure 


The authors wish to thank Professors Goldberg and Bogdanoff 
for their description of their paper and their stimulating discus- 
sion. Statistical approaches have been used with success to pre- 
dict results based on past measurement in a number of technologi- 
cal areas. In particular, the theory of random vibrations has 
received considerable attention in aircraft and missile applica- 
tions, where vibration caused by excitation due, for example, to jet 
and rocket engine noise, is logically treated as a statistical prob- 
lem. 

Clearly the objectives of our paper are quite different from 
those of the discussers, for the present paper concentrates on the 
determination of the lower natural frequencies of an idealized 
model of a drill string and the manner in which certain changes in 
the system alter the natural frequencies. 

We do not wish to underestimate the practical difficulties of 
our assumptions, indeed we have already pointed many out. 
However, the degree to which the assumptions are reasonable will 
vary greatly from hole to hole, depending on such factors as 
the hole depth, the crookedness of the hole, the type of drilling 
fluid used, the nature of the buckling in the collars, and so forth. 
All these factors will also complicate a statistical analysis. 

The discussers point out that their equation (17) does not agree 
with our own results. The reason, of course, is that they have 
assumed different boundary conditions from those that we used. 
In particular, they have assumed the bottom end of the string 
effectively free to longitudinal vibrations. Our best present 
analyses do not support such a result; nor yet do we feel it likely 
that the bit, carrying mean loads approaching 100,000 Ib, should, 
in any ordinary circumstances, continually bounce clear of the 
rock. 

In summary, then, we agree with Professors Goldberg and Bog- 
danoff that a treatment of the drill-string vibration problem by the 
techniques of random vibration theory might well prove interest- 
ing. We must also agree with them that (as they pointed out in 
their original paper) significant calculations by these techniques 
are out of the question at the moment, as the basic experimental 
information necessary is not available. 


Transactions of the ASME 


> 
a 
where, for example, 
| « 
; ) 
|| 


1. FINNIE 


Shell Development Company, 
Emeryville, Calif. 
Assoc. Mem. ASME 


J. J. BAILEY 


Shell Development Company, 
Houston, Texas 


An Experimental Study of 
Drill-String Vibration 


This paper describes equipment which has been developed for the measurement of 
axtal force, torque, and axial and rotational motions at the top of a drill string. Meas- 
urements made in two wells with this equipment represent the first time these quantities 


have been recorded during drilling. 

A large number of frequencies of vibration were found in the recordings. Some 
of these, particularly torsional readings, correlated with predicted natural frequencies, 
but many did not. Several explanations have been proposed for these ‘‘extraneous’’ fre- 
quencies, but no completely satisfactory solution is available. In addition, some in- 


A... mechanical systems possessing mass and elas- 
ticity have a number of natural frequencies of free vibration. 
Associated with each of these frequencies is a natural mode (i.e., 
shape) of vibration. If a mechanical system, such as a drill 
string is acted upon by fluctuating forces whose frequency is close 
to a natural frequency of the system, vibrations of very large 
amplitude may result. When the exciting frequency is equal to a 
natural frequency, the amplitude of vibration is limited only by 
damping (i.e., energy dissipation) in the system.! If the damp- 
ing is insufficient, failure occurs. Many records of this type of 
failure can be found in the engineering literature (diesel-engine 
torsional vibrations, steam-turbine blade vibration, propeller- 
shaft longitudinal vibration, and so forth). 

This type of failure is brought about by external application of 
forces to the system. Another broad class of vibration problems 
is self-excited vibration. Here, the exciting forces are produced 
by motion of the system itself. No external fluctuating forces are 
required to maintain such a vibration once it has started. This 
type of vibration has been found in aircraft wings, turbine 
blades, suspension bridges, electric transmission lines, and the 
like. 

Since so many vibration problems are encountered in engineer- 
ing, it is rather surprising that they seem to have been of rela- 
tively little importance in drilling. The answer, in part, must be 
that the drill string, although possessing an infinite number of 
natural frequencies, is generally sufficiently highly damped by the 
surrounding fluid and by contact with the formation. Under cer- 
tain drilling conditions, however, very noticeable motion of the 
drill string is produced, and drill-string vibrations do seem to be a 
possible cause of drill-string failures. Certainly, this aspect of 
drilling mechanics appears to deserve investigation. 

To obtain greater insight into the vibration behavior of drill 
strings, a series of tests was made on two operating drilling rigs. 
This paper describes the equipment developed for the tests, the 
test procedure, and the interpretation of the results. An ac- 


1 Strictly speaking, the relationship of the exciting forces to the in- 
ertial forces corresponding to the natural modes must also be con- 
sidered. 
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teresting interrelations between axial and torsional vibrations were observed. 


companying paper? describes the analytical procedures which 
were developed for the prediction of natural frequencies. 


Test Equipment 


In studying the behavior of a drill string in longitudinal vibra- 
tion, for example, it is desirable to measure the longitudinal dis- 
placements at all points of the string as functions of time. If this 
is not possible, it is desirable to measure at least the displacement 
E and the force at both ends of the drill string as functions of time. 
The derivative of the displacement £ with respect to the distance z 
along the string, 0&/0z, is the longitudinal strain in the drill string 
and is hence proportional to the axial force. Instrumentation 
which. enables the continuous recording of downhole-vibration 
data is not available. Since the object of the tests was the meas- 
urement of drill-string vibration and not the development of 
radically new equipment, vibration observations were confined to 
the top end of the drill string. 

Torque, axial force, and rotational and axial displacements were 
measured as functions of time. Wire-resistance strain gages 
were used for torque and axial-load measurements, which were 
made below the rotary table. For displacement recording, the 
choice of a transducer was less clear cut. Either a velocity-sensi- 
tive pickup with an integrated output or an accelerometer signal 
integrated twice will yield displacement. Direct-reading displace- 
ment pickups are not convenient because of their low natural fre- 
quency and large size. Generally speaking, an accelerometer re- 
quires more elaborate instrumentation than does the self-generat- 
ing type of velocity pickup, particularly if displacement records 
are required. However, accelerometers have a flat response from 
zero frequency to very high frequencies, whereas velocity pickups 
drop off rapidly in output at low frequencies. In the present 
equipment, velocity pickups were selected for axial and rota- 
tional displacement measurements, principally because the 
pickups and their associated electronic equipment were available 
in the laboratory. The equipment was designed and constructed 
on a rush basis for testing certain wells. This led to a choice of 
units and design of equipment suited to availability and rapid 
construction rather than to strictly optimum performance. 

The test equipment consisted of four main sections, Fig. 1: 

1 A sub (6!/s-in OD by 26-in length) which attached to the 
bottom end of the kelly. The sub contained strain-gage bridges 
for the measurement of axial load and torque. The strain gages 
were connected to a cable about 50 ft long which passed up 
through the kelly. 


2 J. J. Bailey and I. Finnie, ‘‘An Analytical Study of Drill-String 
Vibration,’’ published in this issue, pp. 122-128. 
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Fig. 1 Schematic layout of test equipment 


2 Asub (73/y-in OD by 36-in length) which attached to the 
top end of the kelly. A pressure-tight gland at the bottom end of 
this sub led out the cable from the bottom sub. A heavy flange 
was welded on the sub as a mounting for the slip-ring assembly 
and velocity pickups 

3 The slip-ring assembly, Fig. 2, which had an inner section 
that bolted to a flange on the top sub. There were leads to the 
inner section from the velocity pickups and from the strain-gage 
bridges on the bottom sub. The outer section of the slip-ring as- 
sembly remained stationary, and the inner section rotated. 
Electrical connections between the two sections were made by 
brushes which rubbed on the slip rings. A 50-ft cable led from a 
plug on the stationary part of the slip-ring assembly to the re- 
corders. 

1 Conventional equipment which was used to record the 
strain-gage and velocity-pickup signals. The rccorders were 
placed in the driller’s doghouse. 


In practice, this equipment performed very satisfactorily. 


Test Procedure 


The procedure used was relatively simple. When the equip- 
ment was designed, the sizes of the subs were chosen to match a 
spare kelly. This kelly was moved to the first test location and 
was placed on the pipewalk. The subs were then made hand- 
tight to the kelly, and the gland in the top sub was sealed around 
the cable from the bottom sub. (The two subs and the kelly were 
intended to be handled as a single unit, but a downhole fire at 
the first test location destroyed the strain gages on the lower 
sub, which had to be removed and replaced by a spare lower sub. ) 
Drilling continued while the subs and spare kelly were being as- 
sembled so that no rig time would be lost. 

On the next trip out of the hole to change bits, the regular kelly 
was laid down and the kelly-subs combination was picked up. 
The subs were tightened on the kelly and the kelly was then held 
in the kelly bushing while the slip-ring assembly was lowered by 
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hand over the top sub. After the slip-ring assembly was bolted 
to the top sub and the velocity pickups were attached, the whole 

unit was picked up and was placed in the rathole; thereafter, it 

was treated as a conventional kelly. Some rig time was lost dur- 

ing this phase of the operation, mainly because of the time re- 

quired to change the kelly. At the second test location, this 

procedure was simplified by attaching the slip-ring assembly be- 

fore the kelly was changed, see Fig. 2. 

In subsequent drilling operations, the slip-ring housing was left 
free to rotate. Before data were recorded, the rig was stopped 
briefly, a cable was attached from the recorder to the slip rings, a 
rope was attached to prevent rotation of the slip-ring housing, 
and the bearings were lubricated. 

During the tests, readings were taken over as large a range of 
weights on bit and rotary speeds as could safely be attained. Con- 
ditions were usually steady, but, in some tests, vibrations of the 
rig and drill pipe could plainly be seen. For example, Fig. 3 shows 7 
a type of drilling-line vibration observed on several occasions. 

For the most part, the discussion that follows will be confined 
to three sets of tests—A, B, and C. 

Test A. This well was being air-drilled at 3962 ft with an 83/,- 
in. Reed YC3 (button) bit. The drill string’ consisted of 3445 ft 
of 4!/2-in., 16.6-lb DP and 535 ft of 68/,-in-OD, 23/,-in-ID DC, 
with stabilizers (three straight blades) at 150 and 240 ft above the 
bit. 

Runs were made at three different rotary speeds and with four 
weights on the bit with a lower sub which had been carefully 
calibrated. 

TestB. This well was being drilled with mud at 3496 ft. The 
drill string consisted of 3353 ft of 6-in. DP,‘ 60 ft of 73/,-in-OD, 
3-in-ID DC, and 95 ft of 11'/,-in-OD, 3-in-ID DC. A 12'/,-in. 
OSC conventional bit was used. 


’The lengths given are those used for calculations; they make 
allowance for subs, kelly, and so forth. 

4 This was actually 6-in. casing which was used as dril! pipe (6.000- 
in. OD, 5 424-in. ID, 18.4 lb/ft without tool joints, 19.64 lb/ft in- 
cluding tool joints). 


» 


Fig. 2 Slip-ring assembly on the top sub. The sub is attached to a kelly 
(in the protective sheath). The whole unit is ready to be lifted. 
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{uns were made at four speeds and with four bit weights. The 
frequency of the rig power supply was very low during the tests 
(about 56 cps), and the data have been corrected for the resulting 
low chart speed of the recorder. An older version of the lower 
sub was used, because the other sub had been destroyed by a 
downhole fire at location A. The resistance of the strain gages 
on this sub to ground was not as high as would be desirable. In 
addition, one of the axial gages was damaged, and the remaining 
axial gages had to be rewired as a half-bridge. For this reason, 
the axial-load recordings made in runs B and C may include 
components due to bending. Although the results from this sub 


Fig. 3. Flapping drilling lines 
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may not be reliable for quantitative interpretation, they are 
probably qualitatively satisfactory. 

Test C. These tests were run in the same well as those at loca- 
tion B. The string consisted of 3175 ft of 6-in-DP, 321 ft of 73/,- 
in. DC, and 95 ft of 11!/,in. DC. The bit was a 12!/,-in. Se- 
curity conventional 84. This set of tests was the most complete; 
runs were made at seven rotary speeds and with six bit weights. 

In addition to these three sets of tests, other runs were made. 
At well A, recordings were made as the rotary table was started 
and stopped and while the weight on bit was being changed. In 
test B, runs were made while the speed was changed (i.e., with 
continuously varying speed instead of with discrete values only). 


Discussion of Results 


In vibrating mechanical systems, the frequencies that pre- 
dominate are often the natural frequencies of the system. It is 
difficult to make quantitative strain-gage readings under field 
conditions. In the present tests, as a result of destruction of the 
first sub by fire, the later tests were run with a sub of questionab!e 
accuracy. However, although the amplitudes of the vibra- 
tions may be subject to considerable distortion, the frequencies 
cannot be changed (though their relative prominence can be 
altered and extraneous frequencies may be introduced). There- 
fore the first step in considering the data was an analysis of the 
frequencies that occurred and a comparison of these with calcu- 
lated natural frequencies. 

Figs. 4, 5, and 6 are reduced copies of parts of the recordings. 
The traces show (starting from the top) axial displacement, tor- 
sional displacement, axial load, and torque. Fig. 4 is from test A, 
with 5000 Ib on the bit at 96 rpm. Slight shaking of the kelly and 
of the rig floor was observed. This chart is a fairly typical ex- 
ample. Fig. 5 is from test C, with 40,000 lb on the bit at 40 rpm. 
The observation made during the test was “‘intermittent unsteadi- 
ness on floor—not violent—derrick top shakes intermittently.’’ 
This figure shows certain features of special interest which will be 
discussed later. Fig. 6 was taken at location A while the rotary 
table was being started with 6000 lb on the bit. 
the three traces are as follows: 


The scales of 


Axial Torsional 
displacement, displacement 
in./em radians/em lb/em 
0.005 4000 
0.005 
0.005 16,000 


Load, Torque, 

lb-ft/em 
1000 
2000 2 
1000 


Fig. 


Calibration questionable. 


= 
Fig. 4 Records from test A, 96 rpm, 5000 Ib on bit 
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Fig. 4, the more or less typical chart, shows that the traces, al- 
though periodic, are not regular—that is, the signals contain im- 
portant periodic components, but the level of random signals is 
relatively high, or the periodie signals are of quite rapidly chang- 
ing amplitude (the latter can be seen very clearly in Fig. 5). 
The frequencies in the traces were analyzed by inspection with 
standard techniques.’ They are probably correct within 10 per 
cent and are generally much more accurate. In addition, an 
attempt was made to read the greatest amplitude of each fre- 
quency component as well as that of the whole trace. 

Before comparing the observed and predicted natural fre- 
quencies, it will be necessary to discuss briefly the method of 
calculating natural frequencies and the manner in which the re- 
sults are plotted. A detailed derivation of the technique is given 
in another paper.? 

Consider first the case of torsional vibration of a stepped drill 
string. If the boundary conditions are known (e.g., the drill string 
is free to vibrate at the bottom end and is clamped at the top end), 
the natural frequencies can be determined by trial and error. 


* See, for example, R. G. Manley, ‘‘Wave Form Analysis,’’ Chap- 
man & Hall Ltd., London, England, 1950. 
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Fig. 5 Records from test C, 40 rpm, 40,000 [5 on bit 


Fig. 6 Records for start-up of rotary at well A, 6000 Ib on bit 


The relative amplitude of the vibration within each section of 
the string is represented by sin Q, say, where 2 is an angle as 
shown in Fig. 7. In addition, the ratio of the angular displace- 
ment to the torque at any point is directly proportional to tan Q. 
Therefore for a free end, Q = 90 deg (or 270 deg, 450 deg, and so 
forth), and for a clamped end, 2 = 0 deg (or 180 deg, 360 deg, and 
so forth). Free and fixed ends are the simplest types of boundary 
conditions; intermediate conditions may also occur. 

In practice, the most convenient method of determining natural 
frequencies is to use the known boundary condition at one end 
and to choose a frequency; the calculations are then carried out 
from the known end until the condition at the other end is deter- 
mined. If the calculated boundary condition agrees with the 
known boundary condition, the assumed frequency is a natural 
frequency of the system. If it does not agree, another frequency 
is chosen and the process is repeated. For example, if this pro- 
cedure were followed for the stepped string shown in Fig. 7, 
starting at the bottom end, and too high a frequency were chosen, 
the result at the top end might be Q = 390 deg, say, instead of the 
required value of 2 = 360 deg. This approach is more con- 
venient than starting with both boundary conditions and solving 
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Fig.7 Variation of torsional displacement with length for a stepped drill 
string with free bottom end and clamped top end. The solid line shows 
displacements for vibration in the second mode. The broken line shows 
the displacements predicted if a somewhat higher frequency is chosen. 


for frequency, especially since, in practice, the boundary condi- 
tions may not be accurately known or may be frequency de- 
pendent. A similar discussion can be presented for the case of 
longitudinal vibrations. 

Figs. 8 and 9 show a comparison of predicted and measured 
frequencies for torsional and longitudinal vibration, respec- 
tively, at well A. The experimentally obtained frequencies are 
shown as short vertical lines. No attempt has been made to dis- 
tinguish between displacement and force recordings. 
plots were constructed for the other two locations. 

In the case of torsional vibration, is appears likely that the lower 
end is free, Q ~ 90 deg. When this condition is assumed, the 
slanting lines in Fig. 8 show the corresponding values of 
boundary condition at the top of the string for a range of fre- 
quencies from 0 to 10 eps. Thus if the top boundary condition 
is a fixed end (Q2 = 0), the first six natural frequencies are 0.36, 
1.58, 3.00, 4.45, 5.94, and 7.37 eps. If the equipment at the top 
end of the string can be idealized as a large inertia restrained by a 
very weak spring, it will be close to a fixed end for all but the low- 
est frequencies. For this reason, the regions in Fig. 8 in which 
natural frequencies would be predicted are near Q = 0 deg, except, 
perhaps, at very low frequencies (less than 0.15 eps, say). 

In the case of longitudinal vibration, it appears likely that the 
bottom end of the string is fixed. As in Fig. 8, the diagonal lines 
in Fig. 9 represent the relation between top boundary condition 
and frequency. The equipment at the top end of the string is 
idealized as a spring-mass system, in which the drilling lines and 
derrick are the spring and the traveling block and swivel are the 
mass. A certain force amplitude is required to vibrate this sys- 


Similar 
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tem with a given displacement amplitude at a specified frequency. 
Thus a certain frequency of vibration implies a certain ratio of 
displacement amplitude to force amplitude and hence implies a 
certain value of Q, since this ratio can be shown to be propor- 
tional to tan 2. The top boundary condition that is, in effect, re- 
quired by the spring-mass system is described by the S-shaped 
curves in Fig. 9. The diagonal curves show the top boundary 
condition provided by the drill string. Where the provided 
boundary conditions coincide with the required boundary condi- 
tions, i.e., where the diagonal and S-shaped curves intersect, the 
frequencies are natural frequencies. 

In practice, there was considerable uncertainty in the calcula- 
tion of the top boundary condition. For this reason, two curves 
are shown on each plot—one depending on a calculated estimate 
of the flexibility of the derrick and drilling lines, and the other 
depending on an experimental estimate. 

From an examination of the torsional plot, Fig. 8, it can be seen 
that the natural frequencies of the system are reasonably well 
evident. The frequency of the first mode is probably too low to 
have been observed by the methods used; however, a cluster of 
points presumably represents the second mode, and other ob- 
servations can easily be interpreted as the third, fourth, sixth, 
seventh, eight, and ninth modes. Errors of observation (in cps) 
are probably more or less proportional to the observed frequency 
and are generally less than 10 per cent. On the other hand, errors 
due to incorrect choice of the bottom boundary condition are 
likely to be more or less independent of frequency. If this is true, 
the variation of the bottom boundary condition can be estimated 
from the spread of the second-mode points to be about 30 deg. 
This figure is about the same as those observed for the other sets 
of tests. When the variations in frequency of the lowest observed 
torsional mode were examined in relation to bit weight and 
rotary speed, no distinct trends were found. The single point at 
5.70 eps and the spread of points between 6.15 and 6.90 cps are 
not easily explained. 

The longitudinal vibrations give a less satisfactory plot, Fig. 9. 
Although some points could be interpreted as representing longi- 
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Fig. 8 Top end condition, ‘2°, versus frequency, n cps, for torsional 
vibration at well A 
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Fig. 9 Top end condition, (2°, versus frequency, n cps, for longitudinal 
vibration at well A 
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tudinal modes, many could not. The grouping of points between 
7 and 8 eps is similar to that on the torsional plot between 6 and 7 
cps. The general interpretation of the results of the other two 
sets of tests is somewhat similar. 

A number of observed torsional frequencies show no correlation 
on the torsional plot but agree with the predicted longitudinal fre- 
quencies, and vice versa. Whether this is significant is not clear 
from the data obtained. 

The recognition of natural frequencies has then been only 
moderately successful. The second torsional mode is clear in all 
cases, and other torsional modes can probably be discerned. In 
the case of longitudinal vibration, the natural frequencies are 
not at all predominant. 

Why this is so and what causes the other observed frequencies 
are of considerable interest. It seems unlikely that the readings 
represent natural frequencies that have been scattered beyond 
recognition either by errors of observation or by changes in the 
boundary conditions that are not accounted for. Even if the pre- 
dicted natural frequencies were incorrect, some characteristic 
grouping of the experimental points, which would indicate other 
natural frequencies, would be expected. Although certain group- 
ings do occur, they do not appear to be indicative of the natural 
frequencies of the svstem, since they are very much spread out and 
are not regularly disposed. As regards distribution of the points, 
the experiment favored the higher part of the range shown. No 
points at less than 0.25 eps are shown, and it is unlikely that fre- 
quencies much lower than this could have been observed by the 
methods used. In the plots, frequencies for both stress and dis- 
placement traces are shown together. However, the sensitivity 
of the displacement transducers is greatly reduced at the lower 
frequencies (e.g., 0.38 at 4 eps and 0.20 at 3 eps, if 1.00 at 10 eps 
is assumed). The sensitivity of the stress measurements does not 
vary much over the frequency range shown, so that the smaller 
number of points at the lower frequencies is due, at least partly, 
to the instrumentation 

\n attempt was made to discard possible extraneous data by 
considering only those oscillations with amplitudes greater than 
some arbitrary quantity (to halve the total number of points, for 


example). However, there was no significant change in the over- 
all picture as a result of this selection. All observed points are 


plotted in Figs. 8 and 9. 

If the damping of the motion of the drill string were sufficiently 
great for disturbances to become small (compared with their start- 
ing amplitude) after propagating the double length of the string, 
the coneept of natural frequencies would become invalid. In the 
presence of heavy damping, it would be possible to obtain fairly 
clear traces at frequencies other than natural frequencies, but this 
would require heavy driving. 

In the presence of heavy damping, the motion of the top end 
of the string could be determined by the natural frequencies of the 
surface equipment.- In the case of the derrick and drilling lines, 
the most obvious motion to consider, perhaps, is that in which 
the derrick and lines act as a spring in the manner previously 
discussed. Their natural frequency in the absence of a drill string 
is given in Fig. 9 by the frequeney at which their characteristic 
line crosses the line Q = 90 deg. In the presence of a heavily 
damped string, the frequency could be lower. As there is no 
marked concentration of points at such frequencies (or rather 
lower), this is probably not the important mechanism. 

Torsional and swaying modes could be observed in the derrick 
during tests, but it is unlikely that they would affect the readings 
appreciably 

\nother possibility seems to lie in the flapping motion of the 
drilling lines. On some occasions, this was very pronounced (see 
Fig. 3, in which one of the lines is shown flapping with a double 
amplitude of several feet). In general, all the lines did not flap at 
any one time; this suggests that there are appreciably different 
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tensions in the lines and hence a number of different possible 
natural frequencies. The observed frequency would be twice the 
flapping frequency. Calculated flapping frequencies are in such 
a range (possibly a little too low) as would account for the mass 
of frequeney observations between about 6 and 8 eps. There is 
also some statistical correlation in this range between higher 
observed frequencies and lower bit weights, which would be 
anticipated if flapping were the source of vibration. This mech- 
anism may be connected with the cluster of readings between 6 
and 8 eps. It is unlikely, however, that it would affect the tor- 
sional observations. Indeed, the general grouping of the observed 
torsional frequencies is not the same as that of the longitudinal 
frequencies for the same set of tests, although the observations 
were made simultaneously. 

Another possible source of driving is the mud system. A 
familiar sight on a rig is that of the rotary hose jerking convul- 
sively as pressure surges from the mud pumps pass through it. 
When the pumps are running light, the pulsation dampeners may 
not be in operation, and, as a result, the condition may be par- 
ticularly severe. In the second series of tests, the pump speed 
was about 3.85 strokes per see. The frequency plot shows no 
concentration near this frequency. Harmonies might be expected 
to appear, and some higher frequencies might be construed as 
harmonics. Although this appears to be a rather promising driv- 
ing mechanism, there is no convincing evidence that it is impor- 
tant. The frequency groupings in the first set of tests are in 
many respects very similar to those in the second set of tests, 
although well A was drilled with air, and the fluid pulsations could 
not be important. 

Since the instrumentation and subs were sufficiently well 
tested, it is not likely that they were responsible for the special 
distribution of points (except for the relative scarcity of points at 
the lower frequencies). 

An interesting phenomenon that suggests the complexity of the 
problem can be seen in Fig. 5, which shows about 6 see of a run 
in test C at 40 rpm with a 40,000-Ib bit weight. The second tor- 
sional mode, which has a period of about 3 sec, can be seen clearly 
in the torque (bottom) trace. The displacement traces are not 
sensitive to such low frequencies. An oscillation at about 5 cps 
builds up rapidly on the longitudinal-displacement (top) trace as 
the torque approaches its least values, and dies out, though not 
quite as rapidly, as the torque increases. On some cycles (possibly 
those distinguished by longitudinal oscillations of rather great 
amplitudes), a corresponding oscillation appears on the longi- 
tudinal-stress (second) trace, but this does not always occur (as 
in the first part of the figure). This oscillation builds up and dies 
out with extreme rapidity during the period of least torque, and 
it contains quite large higher-frequency components (about 10 
cps) that do not appear on the displacement trace. These same 
basic features can be seen on a number of charts from both wells 
tested. Fig. 6 shows a start-up test run at well A (the rotary table 
was started from rest at the left-hand end of the figure with 6000 
lb on the bit). Once more, the nature of these phenomena is not 
clear; without further knowledge of them, their importance can- 
not be evaluated. 


Conclusions 

Vibrations at the top of the string were recorded for a wide 
range of conditions. The records were analyzed-in terms of the 
observed frequencies and amplitudes. Some torsional modes 
were observed; the second mode was very clear. No well-defined 
longitudinal modes were detected. Many vibrations were ob- 
served at frequencies that could not reasonably be construed as 
natural frequencies of the whole system, but their causes are not 
understood. Therefore, although the motions and stresses at the 
top of the string were observed, conditions further down the hole 
cannot be estimated with certainty. 
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Perhaps the most interesting phenomenon was the unexplained 
interaction of the torsional and longitudinal motions, as shown in 
Figs. 5 and 6. This occurred on all start-up tests and on most of 
those exhibiting the low-frequency (about 0.3 cps) torsional 
vibration. This effect was observed at both wells. Also of in- 
terest were the groupings of observed frequencies in the ranges of 
5-8 cps and 40-50 eps. Again, these are unexplained. 

Future testing should be directed toward discovering the 
nature of these unexplained frequencies. This is essential for a 
basic understanding of drill-string vibration. The frequencies are 
presumably controlled by some powerful driving mechanism or by 
the natural frequency of some (perhaps unexpectedly impor- 
tant) part of the system. In order to discover the nature of the 
mechanism, special attention must be given to all likely possibili- 
ties (drilling lines, parts of the rig structure, etc.). 
study, separate vibration pickups 
movie cameras—-may be useful. 


For such a 
perhaps even high-speed 


It is also most important to understand the nature of the 
damping in the system and the boundary conditions at the ends of 
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the string. The damping could be investigated by the transient 
decay of sudden changes at the top of the string. This should be 
done with the string rotating, if possible, and with different lengths 
of pipe in the hole. Another interesting experiment would be to 
drive a known oscillation at the bottom of the hole (by using 
a special footstep bearing, for example) and examine the resulting 
motion at the top. If rotary speed were continuously variable, 
driving could be studied at a complete range of frequencies. It 
would be very desirable if the tests could be run in holes being 
drilled with mud and with air. 


Acknowledgment 


The authors wish to acknowledge the participation of Dr. E. G. 
Chilton who initiated the test program and made many contribu- 
tions throughout its development. The field tests were carried 
out with the co-operation of the Shell Oil Company of Canada. 
Assembly of the test equipment was carried out by Mr. D. D. 
Olson who was also responsible for its successful operation in the 


field. 


may 1960 / 135 


, 
h S 
q 
3 
j 
: 
| 
i 
4 he 
Ue 3 
: 
4 
‘ 
ere, 
3 
‘x 


S. J. BECKER? 


Senior Engineer, Westinghouse 
Atomic Power Division, Bettis Site, 


Pittsburgh, Pa. 


L. MOLLICK 


Chief Stress Analyst, 
Loewy-Hydropress Division 
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A general theory for the elastic design of a compound cylinder made by shrinking to- 
gether many concentric cylinders is described. Each cylinder may have entirely dif- 


ferent strength and operate at any temperature under any pressure, internal or ex- 
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comes a problem. 


|. ORDER to design with a factor of safety based on 
yield, it is necessary to have a suitable vield condition. The con- 
dition selected for use is the maximum-shear-yield stress criterion 
because: (a) It has been proven to be reliable for the behavior of 
evlinders under internal pressure [1},? (b) it is simple to apply in 
this case. 

The vield condition, however, is somewhat modified to the 
form: Shear vield in simple shear 1/V 3 &X tensile vield in 
57.7 per cent isa physically 
observable fact, and also can be corroborated theoretically by 


simple tension. The factor 1/4/ 3 


the equally acceptable yield 
eriterion [2], at least for incompressible material. This last 
criterion, however, is more difficult to apply in this case [7] 
The remaining theoretical material for elastic design can be ob- 
tained from [3] (see Appendix 1). It is to be noted that Lamé’s 
equations developed in this Reference are theoretically correct. 
Other thick-eylinder equations in common use by designers are 
either semiempirical or else they already include some yield 
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Nomenclature 


ternal. The only assumptions, in addition to the usual homogeneity and isotropy in 
each member, are that they have the same Young’s modulus, although even different 
elastic moduli could be accommodated, and that they are not so thin that stability be- 


criterion and some end condition and are thus directly derivable 
from Lamé’s equations. 


This must be emphasized, for there is much misinterpretation - 
in handbooks and texts on machine design on the use of thick- 
eylinder relationships. Usually this makes little difference, but 
for the compound tube it is important. 


It is also to be noted that the elastic solution is valid for both 
conditions of plane strain and of plane stress [4]. 
From Eq. (k) of [3], 


9, 
Tmax = 9 
which by using Eqs. (171) becomes 
(a/r)? 
Tmax = (pj — Po) 
which is maximum at r = a. 


Consider two successive boundary radii, k and l, k <7. Then, 
for incipient yielding, 


1 
Tmax = Ty, = (Pe — Py) X 1) 


1 — (k/l)? 


For maximum utilization of material, all successive rings should 
reach their yield strengths simultaneously. Consider boundary 
radii,a<b<e<...g<k<l. Then from Equation (1) 


tangential normal stress in a 


evlinder, positive if ten- — e,f,g,h, member, the latter in the ference at radius of sub- 
sion alphabetical series being script due to method of 
; the larger; also used as subassembly 
¢, = radial normal lvancadliaeelia subscripts to designate the a mean thermal coefficient of 
radial location expansion for range of 
oe Pa = acting radial pressures at temperature change for 
C. axial normal stress in a Pu radii designated by sub- cylinder between radii 
script; also with prime designated by remainder 
and double-prime nota- of expression 
tion to indicate partial as- AT = temperature change in pass- 
cE, € € normal strains in the radial, sembly effects ing from one temperature 
tangential, and axial direc- ro. shear yield stresses in cylin- datum to another 
tions, respectively Teer s+ der between radii desig- r = aradius variable 
nated by subscripts u = radial deflection 
E, v elastic constants—modulus  6,, 6,,... diametral interferences at 
of elasticity and Poisson’s radius designated by sub- —, —- = first and second derivatives 
ratio, respectively script dr dr? with respect to r 
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radii of boundaries of each 


change in measured inter- 
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Since p, is generally known, usually zero, Equations (2) furnish 
the necessary equations for the determination of all the active in- 
terference pressures for maximum utilization of material as well 


as the internal pressure. Thus 


Let f and 
h be fixed and let g be determined such as to make p,; a maximum 
when both rings are at incipient vielding. 

From Equations (2), 


Consider three successive boundary radii, f < g < h. 


Therefore, 


Pr = Th ( 


For maximum p,, 


Op 
PL = Or 
Og 


This yields the condition 


f 
T fo q 


That this is a maximizing condition can be seen by taking a 
second derivative 


(4) 


for all positive values of the arguments. 

Consider three successive boundary radii, f < g < h and con- 
sider the interference 6, between the two cylinders. This inter- 
ference is equal to the difference between the diametral deflec- 
tions at g of the inner and outer cylinders under the action of 
pressures p,, p,, and p,, whether considered as active pressures or 
passive interference pressures (see Appendix 2). 


Compare with [6]. 
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From Eq. (178), [3] (see also Appendix 1), diametral interfer- 
ence 


1—» 4 (Dy pn) 
1 — (g/h? E 


1 — (g/h)? 


g 


2 
(7) 1+, (2 (py — Py) 
1 — (f/g)* 1 — (f/g)? 


(f/g)* py — Py) 


(Pg — Pr) 
g)? 


— (f/g) 


4 -- (pr — Dz) 
E 1 — (g/h)? ; 1 — (f/g)? 


This last equation is quite general and independent of how the 
vessel is designed. However, if the vessel is designed so that 
cylinder fg and evlinder gh vield simultaneously, then (independ- 
ent of what the boundary radii are), from Equation (1) 


4g 


diametral interference, g = radius. 

Equations (3), (4), and (6) are all that are necessary for the best 
elastic design of a compound cylinder. However, in some cases 
the operating design temperatures may be different from those at 
room temperatures, or else interferences are to be measured at 
each stage of shrink rather than free. These additional condi- 
tions can be easily accommodated as follows. 

If the operating temperature is AJ’ above ambient, then, at 
ambient, all dimensions including the interferences are reduced by 
a factor (1 — @AT7’) as a multiplier. Since the pressures and 
stresses are proportional to the interferences, these are also re- 
duced by this factor, but, in addition, the effect of the change of 
elastic medulus must be considered. The elastic modulus in Equa- 
Let E’ be 
Then, to find the pres- 
sures and stresses at ambient temperatures, multiply those at 


tions (5) and (6) corresponds to operating conditions. 
the elastic modulus in the ambient state. 


the operating temperature by - (1 — a@AT). 


4 


Assuming that the interferences have been corrected to ambient 
temperature and measured at this temperature, suppose a sub- 
assembly is to be performed and then the next cylinder (inside or 
outside) is to be machined to the proper interference. Then, the 
proper interference will be different from that measured prior to 
subassembly, for the interference pressure due to subassembly 
will change the inside and outside radii of the subassembly. Thus 
suppose there is a cylinder with boundary radii e << f<g<h<j 
and suppose that subassembly fgh is performed first with interfer- 
ence as calculated by Equation (4) and corrected to ambient condi- 
tions. Then, from Eq. (178) [3], the additional interference 

necessary at f is(1 — aAT) 
1 — (f/g) 


g\? , 
h Shp,’ /E 


——— = first interference pressure and 
1 — (g/h)? 
f, 9g, h are radii. Suppose now subassembly efgh is performed. 
Then, additional interference necessary at A over that in the free 
state is 


, (see Appendix 2) while at 


his (1 — aAT) 


137 


MAY 1960 


a 
b\2 
c 6 = 2 4 
gE’ 
‘ 
- ( =" 
| — (g/h)? 
or 
k \? 
Px kl Pi 0. = (5) 
a \? b\? 
Py — Pg = 1 — 
§ 
g 
tlt 
h 
f\’ 
\ + Pr 
2 2 
g 
og? h? 
i4 


| 


l1—(g h)? (f h)? 


where p, first interference pressure at g, p,” = subassembly 
interference pressure at f 

The general method for any order of assembly can easily be ex- 
tended in either direction. 

The subassembly passive interference pressures (calculated at 
operating temperature) can be found from Equation (5) but not 
Equation (6). Thus for the peculiar assembly order that has been 


used for the first subassembly p,’ = p,’ = 0. Therefore, 
GY) 
_ £6, h 
4g 1 — (f/h)? 


For the second subassembly pe” = pr,” = 0. Then, 


E(6, + A6,) — (e/f)*) X (1 — (f/h 


4f 1 — (e/h)? 


Furthermore, for interest, although not included in the formula 
as written prevjously, 


(f/a)* (f/h)? 
1 — (f/h)? 


” ” 


Thus, an alternative expression for interference at h is 


” 
On + (1 — aAT 
E 1 — (g/h)? 

after the second subassembly 

It can easily be seen that the method can be extended in- 
definitely 

Example: A three-tube container is to be designed with 11!/, 
in. [ID to sustain an internal pressure of 75,000 psi with a factor 
of safety based on yield of one and one third. The shear yield 
stresses for the operating temperature of 850 F above ambient 
are given as 


Inner tube 75,600 psi 
Middle tube 56,600 psi 


Outer tube 49,500 psi 


The maximum diameters are to be 17, 26!/2, and 35 in. 

The assembly order is to be inside tube with middle tube, sub- 
assembly with outer tube and outer tube is to be machined for 
interference after subassembly. Using acoustical data,‘ the 
moduli of elast icity are 


24 X 10° psi at operating temperature 
30 X 10® psi at ambient. 


Design the container, find the proper interferences in order of 
assembly, state the approximate shrink temperatures, and find 
the proper value of the LD so that it will be 11!/, in. after assembly 
without further machining. 

To take care of the faetor of safety it is merely necessary to 
caleulate with a hypothetical operating pressure of 75,000 x 
100,000 psi 


From Equation (3) for simultaneous yielding 


2c 
49, 500 (=) | 
awl 


*Great care must be taken in the selection of the proper value of 
elastic modulus; acoustical data appear to be most reliable. 


138 MAY 1960 


2b \? 2c \?} 
p, = 56,600] 1 — + 49,500] 1 — 
2d 
2a \? 2b \? 
P, = 75,600} 1 — + 56,600 | 
2c 
+ 49,500] 1 — 
2d 


From Equation (4), for maximum p, under simultaneous yield 


2a \2 2b \2 2c \2 
75,600 = 56,600 = 49.500 [ - 
2b 2c 2d 


These give 


2a \? 
DP, = 75,600 + 56,600 + 49,500 — 3 x 75,600( — ). 


2h 
With p, = 100,000, the solution of this equation makes 
2a © 11.25 
— = 0.6002 or 2b = = 18.74 in. 
2b 0.6002 
But, this is larger than permitted. Therefore, make 2b 17 
in. 
2b 
Then, 


Za 

P, = De raf = = 100,000 — 75,600 0.5621 
Ih 

or 


2b \2 2c \? 
Pp, = 57,505 56,600 | — + 49,500} 1 — 
2c 2d 


This guarantees simultaneous vield of all cylinders. In addition 
for maximum p,, 


2h \2 2c 
56,600 = ) 49,500 


Hence, 


2b 
57,505 = 56,600 + 19,500 — 2 * 56,600 x (¢ ) ‘ 


This has the solution 


2b 17 
—— = 0.6502, or 2¢ = ——— = 26 1n. 
2c 0.6552 


This is permissible. Then, from either the maximizing pressure 
equation, or the simultaneous yield equation, 2d may be found. 


Thus, 
2c \2 56,600 0.6552 0.490: 
= — .6552? = 0.4909. 
2d 49,500 


Then, 2d = 37 in. But, this is larger than permitted. Then 
make 2d = 35 in. and adjust 2c such as to satisfy the condition of 
simultaneous yielding with an internal pressure of 100,000 psi. 
As before, p, is known to be 57,505 psi. Therefore, 


2 Ie \2 
57,505 = 56,600] 1 — >) + 49,500] 1 — 
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56.6 + 49.5 = 48.595 
2¢ 35 
This can be solved as a quadratic in (2c)? and yields imaginary 
The conclusion is that an OD of 35 in. is not sufficient to 
sustain 100,000 psi internal pressure with the materials used, and 
with other limitations stated, under the best design conditions. 
Therefore, as a best design compromise, make OD 37 in. 
the nominal design diameters are 11!/,, 17, 26, and 37 in. 
The diameters have now been selected so that no tube will vield 
before any other does and so that the assembly can sustain an 
internal pressure of 100,000 psi before any yielding. For the 
Therefore, 


roots. 


Hence 


condition of simultaneous vield, Equation (6) applies. 


2x 
6, [56,600 — 0.4379 * 75,600] = 0.0333 in. 
24 X 108 
and, 
2X2 


{49,500 — 0.4275 & 56,600) 0.0548 in. 


24 108 


For this steel and the temperature range between operating 
and ambient temperature, @ = 7.8 X 10~¢ in. per in. per deg F. 
Therefore, the temperature correction factor on dimensions is 1 — 
7.8 X 10° &K (850 — 70) = 0.99392. This need not be ap- 
plied to the nominal diameters, for a slight change in nominal 
diameters will not change any of the previous figures appreciably, 
while it is very convenient to retain the rounded figures for 
diameters. However, it should be applied to the interferences, 
giving, for the interferences in the ambient state, 6, = 0.0331 in., 
6. 0.0545 in. 

For convenience, the first interference pressure is calculated in 
the cold state using Equation (5), 


2x17 p,’ — 0 (4 
0.0331 
30 x 108 ( 17 \2 17 
26 


34 
30 x 1087 


= ’ (1.7467 + 0.7790] 
Therefore 
p,’ (eold) 11,563 psi 


It is obvious that this low interference pressure will not cause 
reversed vielding in the inner tube at zero, internal pressure. At 
operating temperature it is even lower, being 
24x10 


X 11,563 = 9,307 psi 


The additional cold-state interference at cis [from [3], Eq. 


(178)], 
17 \? 
2X 26 11,563 


- | = 00150 in. 
30 10°, — (17/26)? 


c 


Theretore 
é. (cold) = 0.0545 + 0.0150 = 0.0695 in. 


In practice, these last figures would be rounded to the nearest 
thousandth of an inch. 
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For the shrink temperatures many sources are available. Merely 
as a matter of personal convenience, the source of data chosen 
here is [5]. 

0.033 
17 
Add to this, datum value of 0.51 at 70 F taken as ambient, to 
give 2.84 in. per 100 ft. By linear interpolation in the table, the 

first shrink temperature is 360 F, approximately. 


0.070 
For the second shrink, e = - - X 1200 = 3.23 in. per 100 ft. 


For the first shrink, e = % 1200 = 2.33 in. per 100 ft. 


Add to this 0.51 as before to give 3.74 in. per 100 ft. Then, by 
linear interpolation, second shrink temperature is 450 F, approxi- 
mately. 

After assembly, the interference pressures (cold) can be found 
from Equation (5) and from [3], Eq. (172) or (175), depending on 
the order of assembly. Equation (175) will be used in this ex- 
Thus 


ample. 


2X 26 
30 108 y (= ) 
2d 2c 


Note the implied use of the principle of superposition by ignor- 
ing interference pressure at radius b. 

From this p,” = 18,310 psi cold 
and 


0.070 = 


0.5621 
11,563 + 18,310 
O.8128 


= 24,225 psi cold. 


The reversed shear stress on the inner tube due to shrinkage is 
24,225 /0.5621 = 43,100 psi in the cold state, perfectly safe! The 
other shrink stresses, of course, must be safe for they would neces- 
sarily have to be much less than at operating conditions. 

From [3] Eq. (178), the diametral contraction due to p,” at 
p= @is 


1 — (a/b)? 


2x 11.25 x 24 295 x 
~ 30 x 100 * 9.5621 


= (0.0323 in. 


pr” X 


Therefore the proper machining ID dimension is 
11.25 + 0.032 = 11.282 in. 


Hence the complete solution to the example given could be 
stated as: 


1 Finish machine inner tube to a bore diameter of 11.282 in. 
and a nominal OD of 17 in. 

2 Measure the OD of the inner tube carefully and finish 
machine bore of middle tube to 0.033 in. smaller than this dimen- 
sion. Finish machine OD of middle tube to a nominal dimension 
of 26 in., heat middle tube to 370 F, and slip over inner tube. 
Allow to cool to room temperature. 

3 Measure the OD of the subassembly carefully and finish 
machine bore of outer tube to 0.070 in. less than this dimension. 
Finish machine OD of outer tube to a nominal dimension of 37 in., 
heat outer tube to 460 F and slip over subassembly. Allow to 
cool. 


The following points of the problem should be noted: 


1 The formulas give diametral interferences in terms of radial 
dimensions, but calculations were in terms of diametral dimen- 
sions. A little caution will avoid confusion between diameters and 
radii. 
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2 There is no need to calculate stresses in the operating 


condition, for the design is such as to guarantee safety. 

3 There is no necessity to follow this assembly or machining 
procedure. For example, all the intérferences could be measured 
prior to any subassembly, and could be applied to either mating 
surface. If a different assembly order is followed, however, the 
shrink temperatures ¢and interferences if measured in assembly 
order) will be different. 

4 Many variations of this problem could easily arise. For 
instance, it might be required to maintain 11!/,-in. bore at operat- 
ing conditions. This would require a machined bore prior to sub- 
assembly of 0.99392 X 11.282 — 75,000 
24x 10°\1 — (11.25/37)? 


033) [from [3], Eq. (179), corrected for diameter} 11.161 


5 The shrink temperatures are relatively small as compared to 
the operating temperature. Therefore, to avoid trouble, the ves- 
sel should be fairly uniformly heated for operation at 850 F. 

6 The values of 7,,,, 7,,, ete., need not be actual shear vields, 
but ean be working shear stresses if desired. This enables one to 
place a factor of safety on the stress rather than the pressure, and 
in fact, to make the safety factors for all tubes entirely different. 
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APPENDIX 1 


From [4], the following equations are valid for elastic action: 
Hooke’s law 


Mquilibrium 


Ceneralized plane strain 


€. constant (9) 


u du 
€, = and by definition €, = , hence 
r 


dr 


Compatibilit 


€.-—€, (10) 


From (7) and (9), 
o, = Wo, + o,) + Ee, 


Therefore, (7) becomes 


Ke (1 — — + — 
> (7a 
Ee, (1 — v2)o, — + v)o, — | 
Inserting this in (10) and using (9) 
do do, 
(1 — v?)r — wl + p)r (1+ (11) 


But from (8), 


and 


do, do, do, do, 
(1 — v?) r? + 27 + —(1 + pv)r- 
dr? dy dy dr 
which simplifies to 
3 do, 
dr? r adr 


ld ( _ do, 
r3 (12) 
r3 dr dr 


which integrates immediately to 


Use is made of the following boundary conditions: 
When r = a, a, —p,;; whenr = b, a, - Po 
Then 


b? — a? 


par 


Hence, 


and from (8), 


pya* pob? 


b? — a? } 


These equations are identical with Eqs. (171) of Reference [3}. 
However, they are redeveloped here in a manner to show con- 
clusively that no indeterminate form or singularity arises at any 
step if the material becomes incompressible (v = 0.5). 

From the second of Equations (7a) 


— — + v)o,—vEe) (15) 
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de, 
a 
j : 
in 
do, 
sf 
dr 
r? +. Dy 
at Inserting these into (11) 
| 
or 
] 
€, 10 vo, + Oo, 
E 
E pa? pob? Pi — Po a*h? | 
1 b? — a? b? — a? | 
€, = lo vio. + (14) 
| 
> | 
do. 
o = 0 (8 
dy 
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The substitution of Equations (14) provides the desired solution. 
Usually €, is not known directly but rather only the axial force, 


b b 
f = f o,rdr 
0 a a 


However, the sum of Equations (14) is aconstant. Hence, since 
o, = Vo, + o,) + Ee, and €, is constant, then a, is also constant 
under generalized plane strain. (This was not assumed pre- 
viously.) Therefore, a, is known once the axial force is known. 
Thus, from the second of Eqs. (7) and using Eqs. (14), 


1—vp p,a? — pob? 


b? — a? 
l+v pj — po a*b? vo, 
E r b* — a? E 


Thus Eq. (178) [3] holds only for an open-ended cylinder 
(o, = 0) if in a state of generalized plane strain. The equation in 
the reference was actually developed for plane stress. 

For plane strain, itself, €, = 0, and from Equations (14) and (15) 


(1 + v)\(1 — 2v) — pod? 
= 
E b? — a? 


i+? OF. 
(47) 
E r b? — a? 

In a long compound cylinder, just as in a long simple cylinder, 
away from the ends a state of generalized plane strain may be 
assumed for each component. Then Equations (14) and (15) re- 
sult in Equation (17) with vre, added to the right-hand side. If 
it is assumed that €, is the same for both sides of an interface 
then, used in a development parallel to Equations (5) and (6), 
vreé, immediately vanishes. The resulting equations are identical 
with Equations (5) and (6) except for the multiplicative factor 
(1 — v?). This factor represents the presence of o, stresses which 
are different on two sides of an interface and which are intro- 
duced so that €, be the same on the two sides. In general design, 
however, this difference usually represents the presence of axial 
shear stresses which thereby violate the strict condition of gener- 
alized plane strain at all sections. 

On the other hand, if it is required that o, be the same on two 
sides of an interface, then using Equation (16) will result in 
Equations (5) and (6) as before. €, will be constant in each com- 
ponent but not the same constant for both. 


APPENDIX2 
Proposition 1 


In a compound Lamé cylinder, the relative deflections at 
mating surfaces (the interferences) are unchanged by changes in 
the internal or external pressure on the assembly. 

Proof: It is sufficient to consider only the two-tube assembly, 
for the change in internal pressure on this assembly can be con- 
sidered as a change in an interference pressure at the same diame- 
ter due to a change in internal pressure in a third inner tube. 
Similarly for the external pressure. 

Thus in Equation (5) let p, and p, be zero and let the corre- 
sponding passive interference, pressure at g be p,’. Then Equa- 
tion reads 

5. = 


1 — (f/h)? 
E 


(1 — (f/g)*]-[1 — (g/h)*] 


Then introduce internal and external pressures p, and pr. By 
Lamé’s equations [3], and the principal of superposition, these 

& The proof for this and the following propositions assume condi- 
tions of generalized plane stress. From Appendix 1, it can be seen 
how analogous proofs hold for generalized plane strain of the as- 
sembly. 
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give rise to an additional interference pressure at g whose value is 


2 
— Pr (4) (py 
L 1 — (f/h)* 


Thus again by Lamé’s equations and the principle of super- 
position, the increase in OD of the inner cylinder is 


jot 


1 — (f/g)? 


+2——g(-—) 
E 1—-(f/g? E 


while the increase in the ID of the outer cylinder is 


(ty 
‘i 1 — (g/h)? E 


BE ° E 


By taking the difference between these last two, the change in 
interference is obtained 


(;) 
4g 9 <9 1 


(y 
g 
(1+ - 
1 — (f/9)* 
f 2 
+ 
(1 — g/h)? 
+p 
1 — (f/9)? 
g \? | 
1 — (g/h)® + 


The last step is not immediately obvious but requires some alge- 
braic manipulation of the terms within the braces. 


Proposition 2 


For a compound Lamé’s cylinder, the interferences given by 
Equation (5)* are valid for the free state. If successive outer 
bands are to be shrunk on, the growth of the OD of the last sub- 
assembly due to the assembly must be added to the free-state in- 
terference at the same diameter in order to maintain the correct 
interference pressures. 
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| 
: = 
a 
( h 
h 
By substituting the value of p,”, Ad, becomes 
; 


Proof: Consider a three-tube container with surface radii a,b, (in that the first subassembly may be treated as one cylinder), re- 


erat c, and d. From Proposition 1 it follows that 6, and 6, may be quires that 
calculated by Equation (5), using only the passive-state interfer- icp, 1 — (a/d)? 
ence pressures. Let these be p, and p,. Then from Equation (5), E a - (c/d)*)-(1 = (a/c)*) 
1} Py p (a b or 
6, = — 
E {1 — (b/c)? 1 — (a/b)? E(6, + Ad.) (1 — (e/d)?)-(A — (a/c)? 
1 — (a/d)? 
and 
a 9 (b/e)Xp, =» Using the value of Ad, in terms of 6, this becomes 

E Fee 1 — (b/c)? | E6, (b\2 (1 — (a/b)*)-(1 — (c/d)? 
4h 1 — (a d)2 
5 On the other hand, if 6, and 6, are given, p, and p, may be KB. (1 ‘a/e)®)-(1 }2 

found from these same equations, and their values are pee Se Se 

4c (1 — (a/d)?) 
b \? a \? 
ia Eb ee d 77 b This is precisely the same as that obtained from 6, and 6, 
4 rr vat D: alone. Furthermore, from [3] Eq. (175), p, gives rise to an addi- 
, =e tional interference pressure at 6 whose value is 


1 — (a/b)? 


h d 


— (a/c)? 


— (a/b)? 


- (a/c)? 


(b/c)?)-(1 — (a/b)?) 
= (a/c? 
fa\? c \? 
Eo, b \? (1 — (a/b)?)-(1 — (e/d)?) 1 — (a/b)? 
Eo. c d 
th 1 — (a/d)? — (a/c)? 
te 1— (a d)? 
E6. (1 — (a/c)?)-(1 — (e/d)? 1 — (a/b)? 
Consider the interference pressure produced by the first sub- 7 ie 1 — (a/d)? i (a/e)? 
assembly, using Equation (5), 
E6, (1 — (a/b)?)-(1 — (b/d)?) 
= 0, = 23 4h (1 — (a d)?) 
(1 — (b/c *)-(1 — (a/b)? 


E6. (1 — (a/b)?)-(1 — (e/d)?) 
1 — (a/d)? 


eS which is exactly the same as obtained from 6, and 6, alone. 
A Therefore, to obtain the correct interference pressures in the 
3 three-tube assembly in which the assembly order is inner sub- 
i From [3], Eq. (178), the growth in OD of the subassembly is assembly first, it is necessary to add the growth of the OD of the 
“ ™ (b/e)2 subassembly to the interference given by Equation (5), at the 
rai E 1— (be? Dy same diameter, in order to obtain the correct final interference 
pressures. 
2 thin he tod, “Then tee Through use of the principle of superposition, this can be ex- 
tended to n tubes. 

Proposition 3 
Aé. ae ti 2 If successive inner bands are to be shrunk into a compound 
Slay c a \* Lamé’s cylinder, the inward growth of the ID of the last sub- 
sii ( c ) assembly is to be added as a positive value to the free-state inter- 
Be ference at the same diameter, as given by Equation (5) in order 
4 If the next subassembly is performed using 6, + Ad, as inter- to maintain the correct interference pressures. ‘ 
a ference, Equation (5), together with the principle of superposition The proof is similar to that of Proposition 2. 
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In April, 1956, the authors presented to The Institution of Mechanical Engineers a 


paper entitled ‘Fatigue Under Triaxial Stress: Development of a Testing Machine 
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and Preliminary Results,’ 
sented to the International Conference on the ‘Fatigue of Metals.” 
reported tests carried out on cylinders made from a 2'/2 per cent Ni-Cr-Mo steel, which 
were subjected to (up to) ten million repetitions of internal oil pressure of (up to) 
Since the publication of these papers a considerable amount of testing has 


’ 


and in September, 1956, a supplementary paper was pre- 
These papers 


been carried out on cylinders made from a carbon steel, a 3 per cent chrome steel, an 
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tion. 


austenitic stainless steel, a light alloy, a nearly pure titanium, the Ni-Cr-Mo steel in 
a harder state, and both the Ni-Cr-Mo steel and the chrome steel in the nitrided condi- 
In addition, many tests of more academic significance have been carried out on 
the Ni-Cr-Mo steel in an attempt to achieve a better understanding of the extraordinary 
results which have been obtained. 
tion of the results (supported, of course, by ancillary tests on each material), which 
are of importance in design. 
are relevant to the practical problem. 


This paper is concerned mainly with the presenta- 


Points of academic interest are discussed only when they 
In order to make the paper reasonably self-con- 


tained, a brief summary of the previous work is given, together with a short description 


” some years the authors have been interested in 
the problem of fatigue under conditions of triaxial stress, and vari- 
ous papers [1, 2, 3, 4, 6]! have described the work which they and 
their colleagues have undertaken in attempting its solution. 
Some of this work is primarily of academic interest, but one of the 
lines of investigation, that indicated in the title of this paper, is 
also of considerable practical importance. The information so 
far published about thick cylinders [1, 2] has related more to the 
experimental techniques than to the results. Since the latter 
paper, however, a large amount of testing has been carried out, 
and this paper is concerned with the presentation of the results, 
particularly in relation to their practical importance, for instance 
in design. It would be illogical and unnecessarily restrictive to 
omit all reference to the academic issues and difficulties which 
these results raise, but to keep the paper to a reasonable length 
the discussion of such topics has been deliberately minimized. 

To make the present paper self-contained and to avoid un- 
necessary reference to the earlier papers, it may be helpful if a 
brief summary is given of the reasons for attempting such tests, 
of the equipment built to carry them out, and of the earlier find- 
ings. 

The elastic thick cylinder? provides an example of a triaxial 
stress system’ which can be regarded as a uniform triaxial tension 
superimposed on a simple shear stress which has a maximum 
value at the bore. The ratio of the triaxial tension to the shear 
changes with the ratio of the external to the internal diameter of 


1 Numbers in brackets designate References at end of paper. 

2? To avoid confusion with tubular members subjected to other 
types of stress, the word ‘‘cylinder’’ is used throughout this paper to 
denote a hollow cylindrical member subjected to internal pressure. 

3 See Appendix. 
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of the machine which has been developed for this work. 


the cylinder (this latter ratio being henceforth referred to as K). 
As K increases the triaxial tension decreases, but can never reach 


zero. Ina bar subjected to torsion, however, the shear stress is, 
of course, unaccompanied by triaxial tension. If repeated 
stresses are used, comparison between the results of these two 
types of test might be expected to show most elegantly the in- 
fluence, if any, of the repeated triaxial tension. 

The resistance of thick cylinders to repeated internal pressure 
is also of considerable practical importance, however, since failures 
of industrial plants subjected to such conditions have proved 
troublesome. The work has indeed received the generous support 
of Imperial Chemical Industries Ltd., and the experimental 
method used has been developed from a suggestion made by Mr. 
W. R. D. Manning of that company. This method required a 
ram to be reciprocated in a test cylinder filled with oil, so produc- 
ing a high pressure each time the oil was compressed; means for 
making up any loss of oil were outlined also. 

The machine built in fulfillment of these suggestions is shown 
diagrammatically in Fig. 1, together with a few details to a larger 
scale. The lower half of this machine consists of a simple but 
robust device for giving the ram a stroke variable in steps between 
zero and 1!/,in. The stroke is adjusted while the machine is at 
rest by altering the phasing of the crankshafts. Above this 
crankcase come, in order, the block holding the main seal for 
the ram, the test cylinder itself, another thicker cylinder which is 
used as a pressure gage, a top pressure-adjusting ram which also 
provides an inlet for the make-up oil, and the adjusting device 
for this ram. 

Items in this machine which merit special attention are: 


The Main Seal. It was from the outset realized that the problem 
of providing a nearly frictionless and wear-free seal to give neg- 
ligible leakage at the design pressure (3000 atm or 45,000 psi) 
would be extremely difficult and might be insuperable. In the 
event, after more complicated designs had failed to withstand 
such arduous conditions for a sufficient length of time, the simple 
“lip seal,” shown in one of the inserts in Fig. 1, was adopted and 
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Fig. 1 


it has proved remarkably successful. Friction is indeed negli- 
gible, wear isso small that a gland will last for 250 million cycles, 
and the leakage, which in the early stages amounts to a few drops 
per minute, may by the end of that time have reached a few cubic 
centimeters per minute. Needless to say this performance is 
achieved only by producing initially an extremely good surface 
finish on, and fit between, ram and seal. 

Core Bar. To reduce the amount of work expended during com- 
pression, the volume of oil in the high-pressure system has to be 
reduced to a minimum without unduly restricting any passages 
through which it has to flow. The design adopted, Fig. 1, uses a 
split core-bar internally tapered and fitting over a tapered exten- 
sion to the ram. 

Pressure Gage and Joints. The test cylinder is connected to the 
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main block by heavy screwed flanges, a wave ring or 4 cone ring 
being used to make the high-pressure joint. The pressure 
gage, which is a cylinder with a higher K than the specimen, is 
connected in a similar manner to the other end of the test cylinder. 
To determine the pressure, the dilatation of the pressure gage 
cylinder is measured by electrical resistance strain gages, using a 
null method [8]. This system is calibrated against a dead-weight 
pressure gage. 

Pressure-Adjusting Ram. To bridge the gaps between the pressure 
steps afforded by the crankshaft phasing, the position of the top 
ram is made adjustable while the machine is running. The effect 
is to alter the clearance volume which is filled with oil. 

Oil Make-Up System. Inevitably there is a slight leakage from 
the gland and this has to be made good. Oil is supplied to a 
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simple spring-loaded nonreturn valve in the pressure-adjusting 
ram. The supply is maintained by a small variable-delivery 
pump; a spring-loaded relief valve limits the maximum pressure; 
an accumulator minimizes the pressure fluctuations; an edge- 
type paper-pack filter insures cleanliness—which is vitally im- 
portant; and a needle valve provides an adjustable restriction 
in the flow to the high-pressure system. The pressure is meas- 
ured before this restriction by a gage and after it by a second gage 
fitted with adjustable electrical contacts on each side of the indi- 
cating hand. Any abnormality in the make-up oil supply or its 
use causes one of these contacts to close; this energizes a relay 
which short-circuits the hold-on coil in the starter on the driving 
motor. This arrangement then insures that the machine stops 
either when a test specimen fails, or when any excessive leakage 
occurs in the high-pressure system, or if the nonreturn valve fails 
to seat perfectly. 

Oil. At the pressures used the viscosity of oil may change a 
hundred-fold or more. It was therefore originally thought es- 
sential to use an extremely thin oil, that chosen being “Shell 
Tellus 15’ with the usual commercial additives. However, for 
particular purposes the same oil without additives, and neu- 
tralized castor oil, have been substituted, and no trouble has been 
experienced on this account. 

Considerable development difficulties were experienced with 
this machine but four units are now in operation—two under the 
control of the authors—and their performance is entirely satis- 
factory. Originally a running speed of 500 rpm was aimed at; 
this has now been doubled with no ill effect, and there is no indi- 
cation that a limit has been reached. 


Materials and Test Specimens 

Six materials have been tested with sufficient thoroughness to 
warrant publication of results. By far the greatest amount of 
work has been done on one of these materials, a 2'/; per cent 


nickel-chromium-molybdenum steel, ‘“Vibrac,”’ in two conditions 
of heat-treatment, with various protective and surface treat- 
ments, originally stress-free and autofrettaged. The other five 
materials are a 3 per cent chromium steel, ““Hykro’’; a 0.15 per 
cent carbon steel; a stabilized austenitic stainless steel; an alu- 
minum alloy DTD 364; and a “commercially pure’ titanium. 
The Vibrac and Hykro steels were supplied in short lengths in- 
dividually heat-treated, the carbon steel in the form of black bar, 
as-rolled, the stainless steel from rolled bar stock, the light alloy 
in the form of extruded bars, solution-treated and aged, and the 
titanium in the form of annealed forged bars. 

Particulars of these materials and their heat-treatments prior 
to machining are given in Table 1. 

To avoid anomalies arising from the stresses set up by the 
machining operations, the preferred procedure when such a treat- 
ment is suitable is to give every finish-machined test piece, 
whether in the form of a thick cylinder or prepared for the ancil- 
lary test program, a final heat-treatment in a high-vacuum fur- 
nace. Since quenching is impossible in this furnace, in the case 
of the alloy steels the temperature selected was 50 deg C below the 
tempering temperature; for the harder (J) Vibrac steel such 
a temperature would have been ineffective and the treatment was 
omitted. Partly for purposes of comparison, a number of tests 
were carried out on the softer (E) Vibrac in which the treatment 
was similarly omitted. The mild steel was given a full normaliz- 
ing treatment at 900 C. No final postmachining treatment was 
given to the stainless steel, the light alloy, or the titanium. 

Preparation of Test Cylinders. A typical test cylinder is shown in 
Fig. 2. The only feature which calls for special comment is the 
method of finishing the bore surface. Grinding was regarded as 
inadvisable, in view of the serious effects of all but the most care- 
fully controlled grinding on fatigue strength. Clearly, however, 
a very smooth finish is required, since this surface is the point at 
which the stresses reach a maximum. Originally (after drilling 
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Fig. 2 Typical test cylinder 


Table 1 


Laboratory 


Materials tested 


Oil 


hardened —Tempered—. 


letter Material name Si Mno 


Percentages 
3 Ni 


0.15 0.66 
0.21 0.60 
0.21 0.60 
0.19 0.45 
0.047 0.66 


0.015 
0.034 
0.034 
0.014 
0.048 


A Vibrac (basic electric) 
E Vibrac (acid open hearth) 
J Vibrac (acid open hearth) 
H Hykro 

M Carbon stcel 


0.013 
0.025 
0.025 
0.015 
0.020 


o 


L Stainless steel 0.70 1.56 0.013 ee 
Solution- 
treated 2 hr 
at 
512 C 

Annealed !/: hr at 700 C 


Cu 
4.20 


Mg Si 
0.49 0.80 
C 0.055 


Ti Aged 5 hr at 
0.012 175 C 
H; 0.002 


G Aluminum alloy 
K Titanium 
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4 
< 
eis & 
Cr Mo (deg C) (hr) (deg C) 
2.55 0.58 0.59 850 1'/, 660 ‘Ab: 
2.68 0.68 0.52 850 4 645 ee, Siar 
2.68 0.68 0.52 850 4 585 Rate 
0.22 3.27 0.51 910 4 625 a 
Ti 


and reaming or boring with a single-point tool) the bores were 
honed to a surface finish of the order 4 microin. However, it was 
necessary to inspect these cylinders most carefully to insure that 
no “‘deep’’ scratches (say, 0.00005 in. in depth) existed. The 
method now used is diamond-lapping as developed by Willcocks 
(Clevedon) Ltd. this produces a mirror-like finish (better than 
1 microin.) in which it is easy to be sure no scratches exist. This 
change in production method has produced no detectable change 
in test results, which suggests that both finishes are satisfactory. 


Results 


The principal results, those from the thick cylinders, are pre- 
sented in the form of S/log-N curves, that is, graphs in which the 
maximum cyclic stress (initially for comparison pressure and hoop 
stress but thereafter shear stress at the bore) is plotted against 
the number of cycles to failure. In the ancillary tests, however, 

large mass of additional information has been accumulated 
about the properties of the materials when subjected to more con- 
ventional forms of stress, for instance static tension and torsion, 
and bending and torsional fatigue; and in the case of the Vibrae 
steel, fatigue in push-pull with varying mean stress of each sign, 
both with and without ambient fluid pressure, and fatigue in shear 
with varying mean stress, both with and without ambient fluid 
pressure. It is clearly impossible in this paper to discuss such 
ancillary tests at length. The values which seem most directly 
relevant to the main tests are mentioned in the text; a compre- 
hensive summary of the normally quoted properties is given in 
Table 2; and further information about the more unusual prop- 
erties is available in references [1-7]. 

Vibrac. To make the story comprehensible, a certain amount of 
repetition is at this stage essential. Thus much of the informa- 
tion in Figs. 3, 4, and 5, which refer to cylinders made of Vibrac 
(A) steel, stress-relieved in vacuo after machining, has been pre- 
sented before: The present figures, however, include a considera- 
ble amount of additional information relating in particular to 


Ib/sq in.) 


Table 2 General summary of results of tests 


the thickest and thinnest cylinders. In each figure, the abscissas 
indicate the number of cycles of pressure applied. In Fig. 3 the 
ordinates are the maximum cyclic pressure (the minimum being 
150 psi, which is regarded as negligible), and there is no reason 
for surprise at the deduction that thick cylinders are stronger than 
thin cylinders. In relation to other fatigue tests, it is fair to say 
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Fig. Maximum repeated pressure in cylinders; Vibrac (U.T.S. 126,000 


(All stresses psi) 


0.15 18/8/Ti 
Vibrac Vibrae RC stainless Al alloy 
Material (soft) (hard) Hykro steel] steel DTD 364 = Titanium 

Lt. of proportionality (tensile) As yield 123000 82000 As yield 11000 63000 3 L000 
0.1 % proof stress (tensile) As yield 130500 97500 As yield 33000 66000 42600 

Yield stress (upper) (tensile) 111300 47500 

Yield stress (lower) (tensile) 109700 32700 
re Ultimate stress (tensile)U 12€000 149000 120600 57100 860004 73000 600002 

Static prop- ‘ 

— Lt. of proportionality (shear) As yield 67000 42500 As yield 9060 28000 17000 
| 0.17 % proof stress’ (shear) As yield 72000 54000 As yield 20600 31400 21700 

| Yield stress (upper) (shear) 25600 

| Yield stress (lower)? (shear) 58700 15200 
Ultimate stress? (shear) 83300 92000 78400 49700 87500 37400 49000 
{ Longitudinal bending(tensile)L +72000 +89600 +72600 +32700 +56500 +31400 +47000 
Fatigue or} Transverse bending (tensile)T +63000 +73500 +58900 + 27600 +55400 +27000 +40300 
endurance } Solid torsion (tensile)S +43700 +52900 +41700 +19700 +26600 +13900 +22800 


(107) lim-) 
its | Thick cylinders (tensile )C 


{Shear/longitudinal tension: 8/L 0.61 0.59 

' Shear /transverse tension: S/T 0.70 0.72 

Patious | Transverse/long. tension: T/L 0.87 0.82 
stress C/I 0.58 0.59 
ratios C/T 0.66 0.72 

| C/S 0.95 0.99 

O/T 0.33 0.35 


@ At astrain rate of 0.0003 per sec. 


(to 31800) 


0 to 41700 0 to 52600 0 to 43000 0 to 19500 0 to 24200 0 to 20000 0 to 25300 


(to 46000) (to 37500) (roughly) (to 28700) 


0.57 0.60 0.47 to 0.56 0.44 0.49 
0.71 0.72 0.48 to 0.57 0.52 0.57 
0.81 0.84 0.98 0.86 0.86 
0.59 0.60 0.43 0.64 0.54 
0.78 0.71 0.44 0.74 0.63 
1.03 0.99 0.91 to 0.76 1.45 1.11 
0.36 0.3 0.28 0.28 0.42 


> These are the true maximum shear stresses in a torsion specimen. From the torque/twist curve a shear stress/strain curve is 
obtained graphically; the ultimate stress is in these cases the stress at failure. 
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Fig. 4 Maximum repeated hoop stress in cylinders; Vibrac (U.T.S. 
126,000 Ib/sq in.) 
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Fig. 5 Maximum repeated shear stress in cylinders; Vibrac (U.T.S. 
126,000 Ib/sq in.) 


that the degree of scatter in these results is satisfactorily small. 
There is every indication of a genuine and abrupt fatigue limit at 
or near the values indicated, and the “knee” of each curve occurs 
surprisingly early, at a number varying in an orderly fashion with 
the thickness from 0.25 X to 0.6 10® cycles. 

In Figs. 4 and 5 these results are replotted with maximum hoop 
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stress and maximum shear stress as ordinates, respectively. 
There is clear prima facie evidence that hoop stress is not the 
criterion of failure, and almost equally clear prima facie evidence 
that shear stress is at least approximately a critical factor. The 
range of shear stress required to cause failure in cylinders sub- 
jected to repeated internal oil pressure is thus defined as 0 
to 41,700 + 1800 psi, and this result holds over the very wide 
range 1.2 << K <3; Le., from cylinders in which the wall thick- 
ness is one tenth of the bore diameter to cylinders in which the 
wall thickness is equal to the bore diameter.‘ 

That this result is not one which can be accepted complacently 
at its face value is amply indicated by the results of other tests 
on this same Vibrac steel presented in Table 3. (As has been 
stated earlier, a very large number of ancillary tests has been 
carried out. The figures in Table 3 should be regarded as a brief 
summary giving average values only.) The values which most 
obviously clash with these for the cylinders are those from tests 
in simple shear stress obtained by subjecting solid or tubular 
specimens to torsion. For reversed shear stress the fatigue limit 
is thus found to be £43,700 or +41,300 psi, respectively; thus 
the range of stress is at least twice that found for the thick eyl- 
inders. For repeated shear stress in thin tubes the fatigue limit 
is 0 to 70,000 psi. It must be emphasized that this last figure 
is obtained by subjecting small tubes only 0.009 in. thick (K = 
1.08) to 10,000,000 repetitions of stress well beyond the primitive 
yield stress of 58,700 psi. However surprising, this is a genuine 
stress value and should not be confused with a nominal stress 
obtained from tests on solid specimens. That the fatigue limit 


4 It should be added that the term shear stress has been used purely 
for convenience in preference to the Maxwell (or Mises-Hencky) 
function. Since these tests involve simple shear and volumetric 
stress the results in themselves have no bearing on any possible 
controversy as to which is the better criterion of failure. 


Table 3 Average properties of Vibrac steels A and E 


(Unless otherwise stated, tests in air at atmospheric pressure) 
] stresses in pounds per square inch 
Statice tensile properties (longitudinal and _ trans- 
verse): 
Modulus of elasticity. . 
Upper yield stress 
Lower yield stress.............. ae 
Nominal ultimate stress................. 
Reduction of area, per cent 
Longitudinal specimens............... 
Transverse specimens............... 48 


29.5 X 106 
111300 
109700 
126006 


Static shear properties (longitudinal and transverse): 
Modulus of elasticity 
Lower yield stress 
Ultimate stress..... 


11.8 108 
58700 


Fatigue limits in tension and compression: 
Longitudinal specimens: 
Rotating bending tests. 
Push-pull tests... .. 
Repeated tension 
Push-pull, in oil bath at 45000 psi, unprotected 
Push-pull, in oil bath at 45000 psi, protected. . 
Repeated tension, in oil bath at 45000 psi, pro- 
tected 
Transverse specimens: 
Rotating bending 


0 to 120000 
+63000 


Fatigue limits in torsion (any orientation of speci- 
men) shear stress: 
Solid specimens: 
Reversed stress................. 
Reversed stress, in oil bath 
Reversed stress, in oil bath at 45000 psi, un- 
protected 
Reversed stress, in oil bath at 45000 psi pro- 
tected 
Thin tubular specimens: 
Reversed stress....... 
Repeated stress.... 
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thus obtained is so much higher than the yield stress makes the 
discrepancy between it and the fatigue limit for the thick cyl- 
inders even more striking. Had the thick cylinders given a 
nominal stress value in fatigue higher than, or even approaching, 
the primitive yield stress, serious thought would have had to be 
given to the redistribution of stress caused by plastic flow—nat- 
ural “autofrettage” so to speak—but certainly in these tests no 
such complications arise. 

These results, then, apparently indicate a most surprising situ- 
ation in relation to the stress conditions giving rise to fatigue fail- 
ure, in that in the absence of triaxial tension the safe range of 
repeated shear stress is 70,000 psi, whereas when accompanied by 
triaxial tension of any magnitude between 4600 and 29,000 psi 
the safe range drops to an almost constant value of 41,700 psi. 
Although the tests were initiated partly to investigate this very 
possibility of an effect of triaxial tension, such a discontinuity is 
so improbable that it is clearly imperative to look for other pos- 
sible factors which might be partially or wholly responsible for 
such results. The search has not been wholly unfruitful, but it 
has not yet been concluded. The story is a lengthy one and, in 
view of the statement in the first paragraph of this paper, in- 
appropriate in full at this juncture; but it would seem undesira- 
ble to make no mention of it. The following paragraphs set out 
in brief the essentials of the situation as it at present appears to 
the authors. 

Investigation of Anomalous Results. (a) It should be re- 
marked that the surface area of the thick cy_nders exposed to the 
maximum value of the stress is very large compared with that of 
the torsion specimens. If, then, relatively rare and randomly 
listributed flaws of varying severity assist in the initiation of 
ailure, the cylinders, on the average, will give a lower result, the 
small specimens showing probably a larger variability. 

(b) An obvious major difference between the types of test is 
that the torsion specimens are tested under normal atmospheric 
conditions whereas the highly stressed surface of the thick cyl- 
inders is in contact with oil, or whatever fluid is used to develop 
the pressure; these oils may, moreover, have contained a small 
amount of dissolved air. Corrosion fatigue is clearly a possi- 
bility which must be considered. In slight support of the cor- 


f 


rosion hypothesis is the fact that, at stress ranges causing early 
failure, multiple cracks are not uncommon. On the other hand 
no trace whatsoever of corrosion has been found by the most care- 
ful microscopic examination; three fluids have been used, ad- 
mittedly all oils, but with noticeably different properties, and no 
difference whatsoever has been found in the results; many differ- 
ent materials which would be expected to have very different de- 
grees of resistance to corrosion show no correspondingly marked 
difference in behavior; and finally, the shape of the S/log-N 
curves is utterly at variance from those normally obtained under 
corrosion-fatigue conditions. 

(c) The presence of oil on the surface of a material subjected 
to repeated stress may have a deleterious effect which is purely 
physical, and this effect may be altered if the oil is at high pres- 
eure. Such an effect would arise from the existence in the 
material, either originally or as a result of repeated stressing, of 
imperfections into which the oil could penetrate. 

The ancillary tests so far undertaken to disentangle these vari- 
ous possibilities have been described in earlier papers [3,4]. A 
number of the findings are relevant to the present discussion: 
(a) The presence of oil, whether under pressure or not, on 
the surfaces of a fatigue specimen has been shown to be definitely 
harmful. (b) This harmful effect can be partially or completely 
eliminated by protecting the surface by a thin layer of rubber ob- 
tained by applying “rubber solution’”’ as manufactured for repair- 
ing inner tubes for tires, and allowing it to dry. (In these cir- 
cumstances, castor oil is used as the pressure fluid.) (c) A speci- 
men so protected and subjected to static triaxial compression is 
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more resistant to additional fatigue stresses than is a specimen 
under normal atmospheric conditions; and it would seem to be a 
reasonable deduction that material subjected to static triaxial 
tension is less resistant to fatigue. 

Protection of Bore Surface. It is a logical step to try the experi- 
ment of protecting the bore of some of the thick cylinders by a 
similar rubber film. The technique is difficult and uncertain and 
has obvious practical limitations. Hence, many other coverings 
have been tried, but the repeated high pressure is quite remarka- 
bly effective in stripping most of these off the metal surface—a 
fact which is itself suggestive. The only satisfactory alternative 
which has been found is a neoprene covering, considerably thicker 
and tougher, and resistant to the “Tellus” oil. Relatively few 
tests have been carried out; even so, the results, which are illus- 
trated in Fig. 6, are of considerable academic interest. 
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Fig. 6 Effect of various surface treatments; Vibrac (U.T.S. 126,000 
Ib/sq in.) 


Neither of the first two factors, the magnitude of the triaxial 
tension or the surface area exposed to the maximum stress, is 
altered by the film; the second two, the physical and chemical 
effects of the fluid, will both be altered. Now the film has a 
marked effect on the strength of the cylinders but does not raise 
it to the value which would be predicted by the torsional fatigue 
strength. It follows that at least one element in each pair of 
factors is operative. In fact all the evidence so far obtained sug- 
gests that the anomalous results which gave rise to this discussion 
and these experiments are not to be explained on the basis of any 
one of the suggestions advanced, but that several are contribu- 
tory; in addition, perhaps, to others which have not occurred to 
the authors. 

Whatever the explanation, the practical problem of designing 
cylinders to withstand repeated internal pressure remains. The 
rest of the work to be discussed is directed towards assisting in its 
solution. 

Effect of Final Honing. ‘The first variant of practical interest 
was discovered accidentally. Some cylinders were scratched dur- 
ing the heat-treatment in vacuo, and the scratches were removed 
by a light final honing. The strength of such cylinders is mark- 
edly higher: Even cylinders with K = 1.4 withstand a pressure 
about 11 per cent greater, while cylinders with K = 1.8 withstand 
a pressure no less than 41 per cent greater, than corresponding 
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cylinders which have been finally heat-treated. These results 
also are shown in Fig. 6. It must be peinted out that the stress 
values indicated in the latter case should be regarded as nominal 
since they are almost certainly above those at which repetition 
will cause plastic flow, which of course implies a certain degree of 
autofrettage. Neoprene covering of these thinner cylinders, 
Fig. 6, is still effective in raising the strength, but not (so far as 
can be judged from barely adequate results) to a value higher than 
that given by rubber-covered, heat-treated cylinders. It should 
be added that optical microscopic examination of the bore surface 
of these cylinders reveals no effect whatsoever of the heat-treat- 
ment; a preliminary examination by electron microscope, how- 
ever, indicates that the effect may be to alter the texture of the 
surface, perhaps even to allow recrystallization of a layer which 
has been ‘‘smeared’’ by the honing operation. Further investiga- 
tions are being pursued. 

Effect of Autofrettage. The second variant investigated was 
the effect of autofrettage; that is, of applying once to cylinders a 
pressure sufficient to cause plastic flow in the bore layers. Two 
K-values were chosen, 1.6 and 2; and the pressure selected was 
that required to cause overstrain to extend to the geometric mean 
radius. A further low-temperature heat-treatment of 250 C 
for 1 hr in vacuo was given to eliminate the consequent “textural” 
stresses without affecting the residual “body’’ stresses. Very 
careful and elaborate preliminary work [6] established the values 
of these body stresses. The results of the fatigue tests are shown 
in Fig. 7, together with an additional curve to indicate the further 
effect of protecting with rubber the bores of the thinner cylinders. 
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Fig. 7 Effect of autofrettage; Vibrac (U.T.S. 126,000 Ib/sq in.) 


From the practical point of view it is perhaps sufficient to re- 
mark that the increase in the pressure required to cause fatigue as 
compared with unprotected nonautofrettaged cylinders is 10 
per cent for unprotected and 20 per cent for protected autofret- 
taged cylinders with K = 1.6, and no less than 26 per cent for 
unprotected autofrettaged cylinders with K = 2. From the 
academic point of view these results are rather more interesting 
and instructive. A major effect of autofrettage is to reduce the 
mean stress in the bore material without affecting the range of 
stress. If this were the only effect the influence on the fatigue 
strength would be calculable from other work on this material 
(4, 5, 7]. On the assumption that shear stress is the criterion of 
failure, the thinner cylinders would derive the maximum possible 
benefit, since they happen to have been overstrained to a degree® 


5 The calculations, which are too long for reproduction here, are 
given in [6]. 
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which corresponds under repeated pressure almost exactly to the 
condition of reversed shear stress (+23,000 psi) while the thicker 
cylinders withstand a stress range of from —30,200 to +22,300 
psi. 

The most plausible explanation would seem to be that an addi- 
tional effect of autofrettage, that of converting the volumetric 
tension to volumetric compression, an effect which increases with 
K, is responsible for the greater strength of the thicker cylinders. 
Certainly (as will appear from later tests) whenever autofrettage 
occurs, either because of deliberate preliminary overstraining as 
in this case, or because the repeated pressure is sufficient in itself 
to cause plastic flow in the bore layers, thicker cylinders seem to 
derive considerably more benefit from the process in comparison 
with thinner cylinders than can be explained on a simple shear- 
stress-range basis. 

Effect of Nitriding. The last variant to be reported in relation 
to this Vibrac steel in the softer state concerns the effect of nitrid- 
ing the bore surface of K = 1.4 cylinders. Again the results of 
the few tests carried out are shown in Fig. 6. It would indeed 
seem that the effect of protecting the bore by rubber or neoprene, 
or by a nitrided layer, is (perhaps fortuitously) to a first approxi- 
mation the same, and the latter treatment is obviously of much 
greater technological importance. However, this is not a steel 
particularly suitable for the nitriding process, and it was decided 
to continue experiments on the effect of nitriding on another 
material to be discussed later. 

Effect of Higher Tensile Strength. It was considered desirable 
to check the behavior of this Vibrac material in a harder state, 
and a tempering treatment of 4 hr at 585 C was selected. The 
nominal ultimate tensile strength rose to 150,000 psi, and as has 
been stated already no heat-treatment was given to the cylinders 
after final machining. The thick-cylinder results are shown in 
Fig. 8 for the only two K-values used, and again a critical shear 
stress is indicated, in this case 52,600 psi. Ancillary tests on 
rotating bending specimens cut longitudinally from the bar give 
a fatigue limit of +89,600 psi and on specimens cut transversely 
from the bar £73,500 psi, while tests in reversed torsion give a 
fatigue limit of +52,900 psi in shear. 

It is tempting to take these stresses at their face value and com- 
pare them directly with those from the softer steel. Again, the 
maximum value, or whole range, in repeated shear stress from the 
thick cylinders is in close accord with the semirange in reversed 
shear stress obtained from torsion tests. Logically, however, 
they should be compared with tests on the softer steel in the “‘not- 
finally-heat-treated’”’ condition and it will be remembered that 
these results were considerably higher and gave different values 
of shear stress for different thicknesses. It may be that, if the 
tentative explanation of a “smeared’’ layer protecting the crys- 
talline material is correct, no such effective smeared layer is 
caused by honing the harder steel. 
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Hykro. The next material investigated was the 3 per cent Cr 
steel known as Hykro. In brief, the two K-values selected for 
the thick cylinders were 1.4 and 2.0, and the results are shown in 
Fig. 9. It will be seen that the scatter is noticeably greater than 
with the Vibrac and that the fatigue limit is by no means well 
defined; perhaps 45,000 + 2,500 psi could be regarded as a rea- 
sonable estimate. However, the results suggest that this scatter 
is in part due to a real, though small, difference between the two 
K-values. Were an estimate of the fatigue limits for these two 
to be made separately the answers would probably be 43,000 psi 
for K = 1.4 and 46,000 psi for K = 2. Now this material is one 
in which the limit of proportionality is noticeably lower than the 
yield stress; the exact value is difficult to define but it is of the 
order 42,500 psi in shear. It seems fairly certain that in these 
cylinders a slight amount of “natural autofrettage” is occurring 
and that this would account (in terms of the explanation advanced 
before) for the difference in the two nominal stress values. 

The ancillary tests gave less scatter than the cylinder tests and 
the fatigue limits are well defined as +72,600 for longitudinal 
bending specimens, +58,900 for transverse bending specimens, 
and +41,700 psi shear for torsion specimens. On this showing 
the Hykro appears to “lose” a slightly smaller proportion of its 
strength when subjected to the thick-cylinder conditions, perhaps 
because of the autofrettage effect referred to, but the difference is 
not great, particularly if the lower limit is chosen for safety. At 
this limit it is again very nearly accurate to say that the whole 
range of repeated shear stress in the cylinders is equal to the 
semi-range of reversed shear stress in torsion. 

Effect of Nitriding. Since Hykro is a suitable steel for nitrid- 
ing, it was considered desirable to trv the effect of this treatment. 
These results also are shown, for thick evlinders with K = 1.6, in 
Fig. 9. The degree of strengthening in comparison with that 
caused in Vibrae was une xpected, as is evidenced by the test re- 
sults; the first specimen was tested at 58,200 psi. (As the ma- 
chine speed was at this time limited to a million cycles a day, it 
will be seen that the five results obtained represent well over 60 
days of solid running of the machine without any allowance for 
changing specimens The result is a maximum nominal re- 
peated shear stress of 65,000 psi and this could clearly be used as 
a basis for design. However, from the academic point of view it 
should be pointed out that since the 0.17 per cent proof and effec- 
tive yields stress in shear of the basic material is only 54,000 psi 
there can be no doubt that the cylinders are autofrettaging them- 
selves and that the true maximum shear stress, neglecting the 
stress due to nitriding, has a lower value; the stress range will 
not be much if at all less than the nominal value. 

Carbon Steel. Carbon steel is a material of such general impor- 
tance that it clearly deserves a place in the program. As indicated 
earlier, a 0.15 per cent C steel was selected, and again K-values 
ot 1.4 and 2 were chosen for the thick cylinders. The results are 
shown in Fig. 10. Again they support the hypothesis of a critical 
shear stress, with in this case a fatigue limit of approximately 
19,500 psi. It is noticeable that with cylinders of this steel the 
fatigue limit is only fractionally below the limit at which gross 
Vielding occurs. At higher values of the repeated stress the fa- 
tigue failure is preceded by diametral growth of the cylinder, 
accompanied by a quite perceptible deformation of the previously 
mirror-like bore suface. 

Ancillary tests on rotating bending specimens closely define the 
fatigue limits at +32,700 psi for longitudinal specimens and 

+27,600 psi for transverse specimens, and tests in reversed torsion 
the fatigue limit in shear +19,700 Ib/in?. Again the similarity 
between the whole shear-stress range in the thick evlinders and 
the semirange in pure shear will be noticed. 

Austenitic Stainless Steel. The last ferrous material to be tested 
in this series is of a very different type from the others; namely, a 
stabilized austenitic stainless steel. Again the K-values selected 
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for the cylinders were 1.4 and 2.0, and the results are shown in 
Fig. 11. On each curve there is some scatter, but this does not 
seem excessive. What is immediately obvious is that, on the 
basis again used for plotting, i.e., the maximum shear stress, 
there is a large difference between the fatigue limits for the two 
sizes, the thicker cylinders giving a result no less than 55 per 
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cent, higher than the thinner. It is of course possible to argue 
that this indicates that the maximum shear stress is no longer the 
basic criterion of failure; if this is indeed so, much more work 
will be required to establish the correct criterion. For the time 
being, it would seem to be more profitable to suspend judgment 
on such a radical reorientation of ideas, and to look for other ex- 
planations of the anomaly. Two relevant items of evidence al- 
ready adduced in this paper immediately come to mind. 

It has been pointed out several times that whenever cylinders 
are autofrettaged, whether deliberately or because the applied 
pressure causes the stresses to exceed the “natural elastic limit’’ 
of the material, thicker cylinders derive more benefit from the 
process than thinner ones. It has been suggested in explanation 
that this follows as a consequence of the fact that, whereas in an 
initially unstressed cylinder internal pressure causes shear stress 
accompanied inevitably by volumetric tension, the process of auto- 
frettage can set up a state of initial volumetric compression of 
magnitude growing with the thickness. Now it is of course well 
known that austenitic stainless steel has very limited elastic 
properties. The present material gives a limit of proportionality 
of about 9000 psi in shear, and the 0.17 per cent proof stress in 
shear (corresponding to the 0.1 per cent proof stress in tension) is 
only 20,600 psi. Now the cylinders with K = 1.4 and 2.0 with- 
stand, at least for 107 cycles if not indefinitely, shear stresses of 
24,200 and 37,500 psi, respectively. It is therefore obvious that 
a very considerable degree of natural autofrettage must occur. 
Indeed the change of bore diameter, for instance, is quite measura- 
ble with an internal micrometer and amounts to something of the 
order of 0.004 in., even in cylinders which have an indefinitely 
long life, and is much greater in those which fail. If the explana- 
tion just given is accepted, thicker cylinders should therefore 
withstand a higher nominal shear-stress range than thinner cylin- 
ders. 

The second item of corroborative evidence, i.e., that a large 
effect can be produced by an apparently small change, is of course 
that already presented in Fig. 6. Here the effect of finally honing 
Vibrac cylinders after the stress-relieving heat-treatment has 
been shown. Cylinders with K = 1.4 are about 11 per cent 
stronger; cylinders with K = 1.8 no less than 41 per cent stronger. 
Clearly the effect is of the same order of magnitude; and if two 
similar unexplained phenomena can hardly be said to explain 
each other, it is to some extent comforting to find that neither is 
wholly isolated. It has been pointed out, of course, that the stain- 
less-steel cylinders were not heat-treated after being machined. 

Ancillary tests have been carried out on this material. The 
poor elastic properties have been noted already. The bending 
fatigue tests could be presented in the usual form of an S/log-N 
diagram and the fatigue limits deduced therefrom at +56,500 psi 
for longitudinal specimens and +55,300 psi for transverse speci- 
mens; the only obvious abnormality would be the extremely flat 
nature of the diagrams. However, it should be stated that the 
behavior of the specimens under test, both in relation to heat 
evolution and an initial tendency to instability, makes it quite 
clear that much reversed plastic flow is occurring in the early 
stages of the test, and that the ability of the material to withstand 
these stresses is associated with a considerable degree of strain- 
hardening. 

The reversed torsion tests are even more informative if more 
difficult to report fairly by giving a single figure for the fatigue 
limit. The machine used approximates to the constant strain 
type. Typically, when the machine starts, the nominal stress as 
deduced from the torque falls rapidly from its initial value by 
perhaps 10 per cent during the first few thousand cycles, steadies, 
and then rises slowly during the remainder of the life to a value 
perhaps 10 per cent greater than the initial stress. The greatest 
range of torque withstood indefinitely (i.e., beyond 20 x 10° 
cycles) by this material under these conditions corresponds to a 
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nominal stress range of from +26,600 (min) to +31,800 psi 
(asymptote). 

Aluminum Alloy. The material chosen is a relatively high- 
strength aluminum alloy to specification DTD 364, giving an ulti- 
mate tensile strength of 73,000 psi. The thick-cylinder tests were 
carried out on one size of specimen only (K = 2); the results are 
presented in Fig. 12. It will be seen that the scatter is very much 
worse than with any of the previous materials and that there is, 
of course, no true fatigue limit. For the sake of obtaining a value 
to compare with other tests, an approximate endurance limit of 
20,000 psi at 107 cycles might be quoted. Reversed bending tests 
also give a considerable degree of scatter, but the endurance limits 
at 107 cycles are fairly well defined at +31,400 and +27,100 psi 
for longitudinal and transverse specimens, respectively; reversed 
torsion tests show that the endurance limit in shear is +13,900 
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Fig. 12 Aluminum alloy DTD 364 


However, the chief interest to the authors of these tests was not 
the magnitude of the stress at which failure occurred, but the 
manner of the failure in the thick cylinders. By this time they 
had become accustomed to the behavior of the steels, in which a 
longitudinal crack perhaps '/, in. long appears on the outer sur- 
face of the cylinder, allowing a fine oil mist to be sprayed into the 
laboratory. In the case of the light alloy the failure is much 
more catastrophic, the first crack which appears stretching from 
end to end of the parallel part of the cylinder and up the radius 
towards the threaded portion. Fig. 13 shows steel and light-alloy 
cylinders after failure. In each case the first illustration shows 
the outer surface. To obtain the second, the cylinder is sawn 
diametrically opposite the fatigue crack, and in the case of the 
steel only in the plane of the crack also, and then broken open. 

Titanium. The last material to be tested in this series is a com- 
mercially pure titanium. The same two sizes of thick cylinder, 
K = 1.4and K = 2, have been used as in most of the other ma- 
terials; the results are shown in Fig. 14. As so noticeably in the 
case of the stainless steel there is a marked difference between 
these curves which are of course plotted on a basis of maximum 
shear stress. It is natural to suspect that the “explanation”? may 
be similar: That this material is one in which the elastic limit is 
lower than the endurance limit. The ancillary “static’’ tests 
make it amply clear that this is indeed the case; the 0.17 per cent 
proof shear stress is 21,700 psi, whereas the endurance limit at 10° 
cycles in shear for even the thinner cylinders is 25,300 psi. Thus 
this is yet another case where natural autofrettage occurs, and the 
thicker cylinders are strengthened more than the thinner. 

The ancillary fatigue tests give the endurance limits at 10? 
cycles in reversed bending as +47,000 psi for longitudinal speci- 
mens and +40,300 psi for transverse specimens, and in reversed 
torsion as +22,800 psi. 
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(a) Vibrac cylinder K = 


(b) Aluminum alloy cylinder K = 2.0 


Fig. 13 Typical failure of cylinders 


Collation of Results and Conclusions. An attempt is made to 
summarize the main features of these results in Table 2. In- 
evitably in such a summary accuracy in detail is lost, but it is 
hoped that the fuller information given elsewhere in this and 
earlier papers is sufficient to avoid misunderstanding on such a 
score. ‘The virtue of such a presentation is that it allows a broad 
view of the whole picture to be obtained by those who are not in- 
terested in the detail. 

Some explanatory remarks are essential. In general the static 
values given are averages from a number of tests and there is no 
real doubt that they are sufficiently representative. The con- 
ventional fatigue results, that is, those from rotating bending 
and torsion tests, to save space have not been shown as S/log-N 
diagrams. A disadvantage of this is that the reader has to rely 
on the brief remarks in the text to obtain any idea of the degree of 
scatter. It gives a fair picture, however, to say that this scatter 
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Fig. 14 Titanium (commercially pure) 


is of the same order of magnitude as that associated with the 
S/log-N curves for the thick cylinders, which have been repro- 
duced. As in the majority of the latter, a single mean line has 
been drawn as nearly as possible through the points and the 
value quoted in the tables is taken from this mean line. In no 
vase have enough tests been carried out to make it possible to de- 
fine the limits of the scatter; in no case, however, is there any 
doubt as to the substantial validity of the value quoted. For all 
the ferrous materials the limit stated appears to be a true fatigue 
limit; for the two nonferrous materials it is an endurance limit at 
10? cycles. 

A few doublet values in the table must be mentioned again. 
Thus the values for torsional fatigue of stainless steel are the 
lowest and highest stresses arising (in the most highly stressed 
specimen which did not break) as the result of a constant strain 
range; the three doublet values in the thick-cylinder fatigue 
limit are those for the thinnest and the thickest cylinders tested. 
Since the discrepancy is believed to be due to the effect of auto- 
frettage, the lower value has been used for the stress ratios in the 
lower part of the table. 

These stress ratios then are the figures to which the authors 
would wish to direct particular attention. The first ratio, that of 
torsional shear fatigue limit to longitudinal tension (or bending) 
fatigue limit is regarded as conventional. Many such values 
have been published; before any conclusions are drawn from 
them, however, it is worth while to glance at the second ratio, in 
which the transverse tensile limit is used as the denominator. 
Since the shear limit (at least for the Vibrac, the only material 
sufficiently thoroughly tested in this way) does not seem to vary 
much with the direction in which the specimens are cut, there is at 
least an argument for saying that the second is the more important 
for academic discussions. The third ratio is no novelty, but the 
information in the literature on the directional properties of ma- 
terials in fatigue is not so plentiful as to make such additional data 
valueless. 

In the next three lines the whole nominal range of repeated shear 
stress in the thick cylinders is divided by the semirange or 
reversed stress in longitudinal bending, in transverse bending, 
and in torsional shear. In the last line of all this safe range of 
repeated shear stress is expressed as a fraction of the nominal ulti- 
mate tensile stress of the material. 

From the point of view of the authors when they started these 
tests it appeared probable that the penultimate ratio might ap- 
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proach, though not closely, the value 2. It isimmediately obvious 
from the table that the value 1 is a great deal nearer the mark, 
and that in the case of the ferritic steels the value 1 is a very 
reasonable approximation which might well be used in design. 
This value would be on the safe side for titanium and on the 
dangerous side for the stainless steel. However, the table brings 
out additional facts; namely, that for all these ferritic steels the 
ratios of the safe repeated shear stress in the cylinders to the 
fatigue limits in reversed bending and even to the nominal ulti- 
mate tensile stress are equally nearly constant. For such ma- 
terials all the designer need do, apparently, to obtain the limiting 
repeated shear stress which may be applied to a cylinder is to 
divide the ultimate tensile stress by 3. Such a procedure has at 
least the merit of simplicity; it would be even more attractive 
were it more easily explicable. The other entries in the table 
somewhat spoil this simplest picture: Clearly it must not be ap- 
plied too indiscriminately. 

Again, since the table does not include all the results which 
have been obtained, it should be mentioned that autofrettage 
is more effective than might have been expected in strengthening 
cylinders exposed to repeated pressure. 

Finally, from the practical point of view, and for those who 
have found the paper too long to read, the effect of nitriding the 
chromium steel is worth recalling. This treatment increases by a 
factor of no less than 50 per cent the pressure which can be applied 
safely to a cylinder of 1.6-diameter ratio. 
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APPENDIX 


Lamé’s equations for the stress in the walls of a thick cylinder 
supporting its own end load reduce to 


P K?R? 

(1) 

P 

=) (2) 
(3) 


where 77 is the radial stress at radius r, 00 the hoop stress at 
radius r, 22 the axial stress which is constant throughout the wall, 
P the internal pressure, R the radius of the bore, and K the ratio 
of outer diameter to bore diameter. 

In these equations tensile stresses are considered as positive 
and compressive stresses negative. 

If a volumetric tensile stress equal to 22 is subtracted from the 
three stress components, then the axial stress disappears and the 
radial and circumferential components or ‘deviatoric stresses” 
become 


K*R? 


(4) 
P K?k? 

(5) 
66 — zz 


These two stresses are equivalent to a shear stress of 


P K2R2 
= — 
(K?—1) r 
which acts on spiral planes parallel to the axis of the cylinder and 
inclined at 45 deg to successive radii. 
Thus the stress system at radius r is equivalent to a shear stress 


r@ with a superimposed uniform triaxial tensile stress equal to the 


axial stress 2z which is uniform throughout the cylinder. The 
Maximum stress occurs at the bore and is given by 
P(K?/(K* — 1)] (6) 
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The Interaction of Textile Structures 
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and High-Speed Textile Processes 


A study has been made of the interaction which takes place in certain high-speed textile 


processes at points of machine fiber contact. 


have been analyzed. 


Bai design of new processing equipment and 


efficient use of new fibers on existing machines require that we 
understand the dynamic interaction between textile materials 
and textile processes. During processing, textile fibers are sub- 
jected to varying applied stresses and strains, caused by machine- 
material interactions. For example, the rubbing of yarn over a 
guide produces stresses at the point of the frictional interaction. 
The development of a whirling balloon [1]? during yarn spinning 
is a result of many dynamic interactions including aerodynamic 
drag. In another case, abnormal strains may occur if a high- 
speed flow of fiber or yarn is suddenly halted because of material 
faults or machine failure. 

A recent survey |2] has pointed up the importance and probable 
magnitude of transient strains produced by rather minor malfunec- 
tions during processing, such as the presence of a slub or tem- 
porary interruption of material flow. strains 

How- 
ever, their role can be completely assessed only when more is 
known about the normal or steady-state operating conditions 
upon which they are superimposed. 


These transient 
increase in magnitude as the process velocity increases 


The concept of machine-material interaction which is developed 
represents a departure from the traditional studies of machine 
actions. Previously, the studies of were concerned 
strictly with the kinematics and dynamics of the machine itself. 


actions 


Cam motions, rotary motions, reciprocating actions, feed and 
take-up mechanisms have all been studied and analyzed by 
These studies 
usually neglect, and rightfully so, the presence or absence of the 
textile material which is being processed. 


mechanical engineers for several generations. 


The sheer mass of 
past and present textile machines is so out of balance with the few 
grams of fiber which are being operated upon that there is seldom 
For this 
reason, this study deals exclusively with what happens to the 
textile material as a result of the interaction. 


any significant reaction by the fiber on the machine. 


It is anticipated that future machinery, designed for operation 
at higher and higher speeds, will be less massive; and a time will 
come when the action of the material on the machine will be 


more important. Of course, the action of material on machine 


becomes important as the mass of the structure is increased. As 
examples, consider the massive drum of a carding engine whose 
rotational inertia is so great that it is oblivious to the presence 


of fiber even when loaded, and then consider a woven fabric 


' This paper represents a portion of an ScD thesis submitted to the 
Mechanical Engineering Department, Massachusetts Institute of 
Technology, Cambridge, Mass. 

2 Numbers in brackets designate References at end of paper. 

Contributed by the Textile Engineering Division and presented 
at the Annual Meeting, Atlantic City, N. J., November 29-December 
4, 1959, of Toe AMERICAN Society OF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be un- 
derstood as individual expressions of their authors and not those of the 
Society. Manuscript received at ASME Headquarters, August 12, 
1959. Paper No. 59-—A-115. 
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The specific actions of weaving and sewing 


which must be held to dimension on a tenter frame. In the latter 
case, a significant reaction is induced in the machine, but the 
process is no longer one which can be called high speed. In 
general, as the textile structure becomes more organized, the 
mass of material being processed increases. The mass flow re- 
mains reasonably constant and hence the absolute flow velocity 
decreases. In addition, as the structure becomes larger it can 
exert a greater reaction force on the machine. 

Schwarz [3] has introduced the use of symbols to portray 
the multiplicity of actions which occur in all textile processes. 
These action symbols are confined to actions which directly 
By the 


use of action symbols, processes, operations, and subprocesses 


involve the fiber, yarn, or fabric undergoing processing. 
are reduced to their essentials. The clarity and simplicity of the 
approach as exemplified by symbols constitutes the 
greatest aid for the textile-machinery designer to reappraise the 
evolutionary machine monsters and replace them with revolu- 
tionary machine midgets. 

If one were to consider the effect of these basic machine actions 
on the fiber, it could be seen that the mechanical actions are all 
variations on the theme of force transmission to the fiber—struck 
by, striking gripping, restraining, pulling. These 
particular mechanical actions produce the interactions which 
will be considered in the following text. The interactions pro- 
duce surface and body forces which are a consequence of machine 
velocity, fiber inertia, contact areas, surface roughness, bulk 
properties, and so on. 


“setion”’ 


against, 


Some of these factors are fiber properties 
and some are machine and process properties. The combination 
of all factors must be considered in the machine-material inter- 
action. In order to define quantitatively and analyze the effects 
of these interactions, we must be able to delineate both the fiber 
and machine behavior at the moment of the interaction. 

Friction is the means of the development of normal longitudinal 
forces by lateral and shear stresses at the point of action. Fric- 
tional or surface forces are again the consequence of both machine 
and fiber properties. The surface forces can be considered to 
be local forces confined to the contact points. In addition, one 
must consider the bulk forces which exist some distance from 
the local action point. The well-known principle of Saint Venant 
allows one to neglect the particulars of the contact area and to 
treat the bulk forces away from the action point as an independ- 
ent problem. 

The mechanical properties of fibers must be related to the 
times of the measurements of these properties. A general dis- 
cussion of fiber properties must include a general description of 
the time-dependent variations of properties. The properties 
which correspond with the time scale of the process action are 
those which are involved in the machine-fiber interaction. 

A quantitative analysis of an interaction cannot be made from 
a knowledge of fiber properties alone. The machine actions 
must be defined both qualitatively and quantitatively. The 
velocities and other details of striking, pulling, sliding, and the 
like, must be known at each action point of a process. 
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Weaving 

One of the few cases in which the concepts of viscoelasticity 
have been applied to a textile process is the work of Greenwood 
and Cowhig [4] in their study of weaving. These workers have 
related the stress relaxation of warp and woven fabric tensions to 
the movement of the cloth fell during a loom stoppage. In 
addition, Greenwood and Vaughan [5, 6] have made a study 
of the beat-up force as a function of several loom variables; 
e.g., loom speed and shed timing. They measured the force- 
time relationship for the beat-up cycle. These are examples 
of investigations into machine-fiber interactions. —, 

First let us consider Greenwood and Vaughan’s study of the 
beat-up force. They applied strain gages to the back of the reed 
in order to obtain measured values of beat-up force. This 
ingenious use of strain gages does not introduce any new variable 
The 
force-time measurements allow us to analyze the amplitude of the 
forcing function as well as its frequency. 


and allows actual in-process measurements to be made. 


Warprension 
Fabric tension 


Open shed 

tension 

Beat up 

peak 
Basic warp 

tension 
Closed shed 

tension 

Beat upminimum 


Closed shed 


open shed 
Beat up Closed shed 


Fig. 1 Typical warp-tension cycle 


The action of the beat-up in relation to the over-all warp- 
tension cycle can be seen in Fig. 1, taken from Greenwood and 
Cowhig. At a loom speed of 190 picks per min the frequency of 
the complete tension cycle is 3.1 cycles per sec. However, the 
time duration of the beat-up force pulse is much shorter than the 
complete tension cycle, and has been found to be approximately 
30 milliseconds [6]. Thus, the beat-up cycle can be thought of as 
having a frequency of approximately 33 cycles per sec, while the 
over-all warp tension cycle is approximately ten times slower than 
the beat-up cycle. 

For the fabric which Greenwood and Vaughn [6] studied, it 
was found that the beat-up force had a peak of about 15 to 20 
grams per 140 denier warp end. Thus the beat-up cycle can be 
considered to have a force amplitude of approximately 0.15 grams 
per denier at an intermittent frequency of 33 cycles per sec. This 
beat-up cycle is superimposed upon the entire warp-tension cycle 
produced by the shedding. ‘The measured values of the warp ten- 
sion for the shed cycle were between 20 to 27 grams per 140 denier 
end. It can be seen from Fig. 1 that the tensions in the unwoven 
warp and in the woven fabric are identical with the exception 
of the short time during which the beat-up force is active. Dur- 
ing this time the warp tension is increased to a value greater than 
the fabric tension. The total beat-up force is, as indicated, 
equal to the difference between the fabric tension and the warp 
tension. While the warp tension varied between 20 to 27 grams, as 
indicated, the corresponding fabric tension on the other side of 
the beat-up reed varied from 5 to 25 grams, approximately. Thus 
the tension variations on the fabric side are more severe than on 
the warp side of the fell. 
when bumping occurs. 


This difference is particularly severe 


Since Greenwood and Vaughan’s measurements are so complete, 
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the stress cycles of the warp fabric tensions are well defined. 
Simple experimental data can be made on yarns in which the 
weaving-stress cycles are reproduced. In other words, the action 
of the machine has been defined quite explicitly and this machine 
action can be applied to hypothetical or real warp and fabric 
materisls in order to reproduce the macbine-fiber interaction. 
The stress-cycle frequency of interest is in the range of 3 to 33 
cycles per sec. This range has been quite adequately covered 
experimentally by Tipton [7]. Thus a great deal of fiber data is 
currently available in the literature and can be applied directly 
to this problem. 

Another interesting application of the principles of viscoelastic 
behavior has also been made by Greenwood and Cowhig [4]. 
In this case, they have applied the principle of stress relaxation 
to the warp and fabric in a loom when the loom is temporarily 
shut down for one reason or another. If the stop takes place 
with a closed shed, the warp tension and the fabric tension are 
equal and both are equal to a tension value considerably above 
the minimum fabric tension which occurs during beat-up. In 
the stopped condition, the tension in both fabric and warp will 
relax to some value less than that which occurred when the stop- 
page was made. The exact nature of the stress relaxation and/or 
creep depends on the nature of the stoppage as well as on the 
If constant tension 
is maintained during a stoppage, then creep will occur, and the 
position of the fell is apt to move. Fell position can move either 
away from or toward the weaver. This motion of the fell will 
cause a setting-on streak when weaving is resumed. 

Let us consider a hypothetical case of loom stoppage, where a 
constant tension is maintained. It is assumed that the stoppage 
occurred with a closed shed in such a way that the fabric tension 
equaled the warp tension. It is important to mention at this 
time that equal fabric and warp tensions do not necessarily mean 
the same yarn tension in both warp and fabric, since the presence 
of any crimp in the fabric will in fact produce a higher axial 
tension in the yarn as it lies in the fabric than in the same yarn 
as it lies straight in the warp. Thus, the tension and hence 
strain condition of the yarns in the woven fabric are apt to be 
higher than the tension and/or strain in the warp. The motion 
of the fell, or weave point, will now depend on whether the fabric 
grows at the expense of the warp or vice versa. If the fell posi- 
tion moves away from the weaver, the fabric has grown at the 
expense of the warp. This is, in fact, the usual case in practice. 

If, on the other hand, the tensioning device maintains a con- 
stant length during stoppage and the stress is allowed to relax, a 
different type of action can take place. The over-all process 
must follow a line of constant total length. The fell will move 
only if the fabric or warp grows at the expense of a contraction 
in the other member. As indicated previously, the yarn tensions 
in the fabric are apt to be higher than the yarn tensions in the 
warp. In the low-strain region, which is the region which is 
important during weaving, the higher the tension the faster the 
relaxation. Thus, it might be expected that greater relaxation 
would occur in the fabric than in the warp. The fabric will then 
extend as the warp contracts, and the fell will again move away 
from the weaver. 


type of tensioning device used on the loom. 


Waesterberg and Nordhammer [8] have questioned the analysis 
made by Greenwood and Vaughan [5], and they suggest that the 
tension might actually increase during a stoppage because of 
so-called stress-restoration effects. They have analyzed warp- 
tension cycles produced on an Instron tensile tester, and have 
shown that stress restoration can in fact occur if a stoppage is 
made at a low force point in the force-time cycle. While this is 
an interesting and worth-while observation, it does not occur in 
practice according to Greenwood and Vaughan. The reason 
given is that a stoppage usually occurs at higher than minimum 
values of tension. An interesting observation by these authors 
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is the fact that the position of the cloth fell might be maintained 
by leaving the beat-up reed against the fabric in such a way that 
the tendency of the fabric to relax more than the warp will be 
overcome. In other words, the beat-up can be used to produce 
an unbalance in the fabric and warp tensions. Under proper 
conditions, no growth or motion of the fell will then occur. 

In all of these analyses, if a complete understanding of the be- 
havior of the fabric and the warp is to be achieved, it is important 
to know the stress-strain-time surface of the warp after it has 
gone through the multiple-load cycling which occurs in practice. 
Also, we must know the same facts for the yarns in the fabric 
as they exist in the loom state; that is, as the fabric is held in 
width with a minimum fabric crimp. In this case the modulus 
of the fabric may not be as different from the modulus of the warp 
as has been indicated by Waesterberg and Nordhammer [8]. In 
spite of this reservation, which is in the nature of a refinement, 
the application of the viscoelastic phenomenon of stress relaxa- 
tion and/or stress restoration is a good example of machine-fiber 
interaction. Thus we see that two types of loading cycles can 
be important in weaving; the one involving the dynamic loading 
which we commonly think of when considering weaving, and the 
other the stress-relaxation type of behavior occurring when the 
loom is shut down. 


Sewing 


The process of high-speed sewing offers another interesting 
example of machine-fiber interaction. In this process, the sewing 
thread undergoes a complex repetitive loading cycle. An at- 
tempt will be made in this section to detail the machine actions 
which are involved in this important process. This will be done 
both statically and dynamically; that is, at very low speeds as 
well as at high speeds with the use of high-speed motion-picture 
equipment. While the machine actions have been defined, no 
force measurements have been made on the sewing thread as it is 
being processed. Some of the difficulties encountered in making 
such measurements will be outlined. 

While the process of cold drawing is a steady tension, flow 
problem, the sewing action consists of severe cyclical applications 
and removals of force. The sewing action is anything but steady 
and continuous. It can be represented in either of two ways: 
(1) By considering the sewing thread to be under zero static 
stress with a loading and unloading cycle applied, or (2) by con- 
sidering the thread to be under a given static stress with an oscil- 
lating stress superimposed in such a way that the stress varies 
between a maximum value and zero. Observation of the sewing 
action of high-speed machines by high-speed motion pictures 
shows that the thread is actually slack over a large part of the 
cycle. Hence, it is probably better to describe the stress cycle in 
sewing by the first alternative just given. 

The problem of sewing has been treated infrequently in the 
literature. In fact, only two references have been found which 
deal with this problem at all. One is the master’s thesis by 
Schranze [9], who studied the actions of the various machine 
parts as well as described some of their functions in forming a 
stitch. Schranze’s work contributes little that can be used in this 
study of machine-fiber interaction. The second reference ap- 
peared as a paper presented by Hardy [10] to the Thread Insti- 
tute and hence has received only narrow circulation. His work 
is an attempt to actually measure sewing tensions and the cyclical 
forces which are developed. Unfortunately, because of the 
limitations of the tension-measuring equipment which was used, 
the measurements were made at !/19 of the top speed of a commer- 
cial industrial sewing machine. A determination of the forces 
developed at 500 stitches per min does not furnish information 
which is directly applicable to the forces which will be developed 
at 5000 stitches per min. Thus, while Hardy’s data are inter- 
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esting, they have little value in the analysis of tensions in high- 
speed sewing. 

Because of the complexity of the actions in a single sewing cy- 
cle, the first study of the thread motion was made by slowly turn- 
ing the machine through a cycle of operation. This was done on 
a Singer Sewing Machine Model Number 600-1. The motion of 
the needle was measured as a function of the angle of shaft rota- 
tion. In addition, the motion of a point on the thread was co- 
ordinated with the needle motion, Fig. 2. There are three lines on 
the figure, two representing the position of the thread point, and 
one the needle position as a function of angular rotation. The 
position of the lowest point of the needle travel was arbitrarily 
defined as the origin. It can be shown that the needle motion is 
very close to a sine curve with a frequency of 83 cps when sewing 
at 5000 stitches per minute. The maximum needle velocity is 
about 25 fps. 
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Fig. 2. Motion of needle and thread from static motion of Singer Model 
600-1 


Let us follow the motion of the thread point, represented by 
the upper curve of the two. At the start of the cycle the thread 
has just been carried through the fabric by the needle. As the 
needle retracts, the thread forms a loop which is hooked by the 
bobbin carriage beneath the bed of the machine. The thread is 
carried down by the rotation of the bobbin hook. The thread is 
moving down ata velocity estimated at 100 fps, while the needle is 
moving up with a maximum velocity of 25 fps. Thus there is a 
differential velocity of approximately 125 fps between the thread 
and the needle during this part of the cycle. The thread reaches 
its lowermost point when it is carried around the bobbin carriage. 
At this point, it appears that a tension peak should develop. Im- 
mediately following the tension peak, the thread slips off the bob- 
bin carriage, but is still engaged by the bobbin hook. An in- 
stantaneous retraction takes place before the take-up mechanism 
in the upper part of the machine begins to control the motion of 
the thread. During this free contraction the lower portion of the 
thread can move faster than the upper portion and the thread will 
go into compression, buckle, and goslack. Because of the uncon- 
trolled free retraction, the actual velocity of the thread during this 
part of the cycle cannot be estimated from the static method 
used. It can be said, however, that the actual velocity might be 
higher than the 190-fps figure indicated. 

In the ascending part of the cycle, thread motion is controlled 
by the upper, or rotary take-up mechanism above the needle. 
Again the thread motion and the needle motion are out of phase. 
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When the thread reaches its highest point, the newly formed 
stitch is pulled into the fabric. It might be expected that a sec- 
ond tension peak will develop. At this point the needle is al- 
ready moving down and shortly thereafter the needle carries the 
thread back into the fabric for a new stitch. It is only over this 
small portion of the cycle that the needle and thread actually 
move together. During the rest of the cycle, the thread is pulled 
either down by the bobbin hook or up by the take-up mechanism. 
In both cases, the thread motion is considerably different from 
and usually opposite to the direction of the needle motion. 

Throughout the sewing cycle, thread is run over many guides 
and machine parts. However, the most severe actions seem to 
take place between (a) the thread and the needle, and (b) the 
thread and the fabric being sewn. No observations were made 
of the interaction between the bobbin hook and the thread nor of 
any interactions which occur beneath the bed of the machine. It 
must be emphasized that these interactions certainly are impor- 
tant and should not be neglected. 

The relative velocities which occur between the thread and the 
needle are higher than those which occur in any other common 
textile process, reaching values close to 200 fps. It should be re- 
membered that very high coefficients of friction occur at high ve- 
locities. This fact is extremely important, since the forces which 
developed in the sewing thread are entirely due to the friction be- 
tween the thread and the machine parts. This is true both in 
the tensioning device, as well as in the take-up mechanism, and 
also at the point of contact between the thread and the bobbin 
hook below the bed of the machine. 

Further studies of the machine action have been made by tak- 
ing high-speed motion pictures of a sewing cycle. In this case a 
Fairchild camera was focused on the upper part of a Singer Sew- 
ing Machine Model 401. The machine was sewing at 2500 
stitches per min and the camera was taking pictures at 1750 
frames per sec. Thus, it was possible to divide one revolution of 
the sewing machine into 42 separate photographs. Again the 
origin for the motion is the point of lowest penetration of the 
needle. This sequence of pictures shows that the sewing thread 
is slack for most of the cycle and is buckled into widely distorted 
forms most of the time. It is possible to estimate tension points 
by observing when the sewing thread appears to grow taut. The 
thread first appears to be taut in Frame 18, which is equivalent to 
approximately 140-deg revolution. 

teference to Fig. 2 (which admittedly describes the thread mo- 
tion of a different machine) indicates that this might correspond 
to the tension peak predicted when the thread is at its lowest 
point after being carried around the bobbin carriage. The thread 
goes slack immediately thereafter, and in the next sequence of 
pictures appears to buckle quite violently. The next tension 
point appears to begin around Frame 30 and continues in the 
next frame. If we look at the left portion of the machine, near 
the tension disks, it is possible to see that the tension check spring 
has been deflected. The significance of this deflection will be 
discussed later in detail. Shortly after Frame 32, the thread is 
seen to go slack again, and it remains slack through the rest of the 
cycle until the tension peak occurring in Frame 17 re-occurs. The 
second tension peak, that corresponding to Frames 30 and 3], 
corresponds to the tension peak shown in Fig. 2, called peak 2. It 
is during this tensioning that new thread is fed into the system 
to take the place of that which has just been put into the stitch. 

From consideration of the thread motion and the tension peaks 
developed as estimated from this series of photographs, it is pos- 
sible to outline some of the difficulties which might be encount- 
ered in introducing a tension probe to try to measure sewing ten- 
sions. In the first place, the tension peaks seem to develop over 
a very short portion of the sewing cycle. In the two instances, 
the tension peaks seem to exist only over two, at the most three, 
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frames of the sequence of pictures. It will be remembered that 
these frames are spaced !/,79 sec apart, and thus the time between 
exposures is slightly more than '/; milliseconds. Therefore, dura- 
tion of the peaks is in the range of 1 to 1'/2 millisecond. Any ten- 
sion-recording device must have a high enough resonant fre- 
quency to measure accurately these very short peaks. 

The second difficulty which must be considered in measuring 
tensions is what will be the effect of placing the tension probe on 
the sewing thread. If the probe is placed between the take-up 
mechanism and the needle, which certainly appears to be the 
most desirable position, it will inhibit the free motion which is 
seen to occur in the normal sewing operation. It is possible that 
inhibition of this free motion will contribute tension levels and 
tension peaks which are different from those which occur when 
the probe is absent. Some workers have been known to place a 
tension-measuring device back near the position of the check 
spring, that is, near the disk-tensioning device. It can be seen in 
Frame 31, and those immediately thereafter, that the tension con- 
dition back at the disk can be considerably different from the 
tension between the take-up wheel and the needle. The addi- 
tional friction which occurs at the different points of contact can 
produce wide variations in tension along the length of the thread. 

The combination of the high-speed motion-picture study, the 
static study, and Hardy’s data gives us an indication of the type 
of stress cycle which is imposed on the thread during sewing. This 
cycle has at least one large peak, which occurs when the thread 
is pulled into the stitch. There can be other minor tension peaks 
which are developed depending upon the setting of the machine 
and the speed at which the sewing is done. No estimate has 
been made of the relative importance of these secondary peaks. 
At extremely high speeds they could reach the magnitude of the 
main peak although no data are available to show this. 

The stress cycle might be considered to be a single or multiple 
peak of about a millisecond duration, which occurs over a small 
part of the sewing cycle. The loading on the thread is far from 
being a sinusoidal cyclical motion, but instead is a pulse loading. 
This must be considered when we try to apply the viscoelastic 
data to the sewing problem. In addition, the velocities at which 
the sewing thread moves over guides is much faster than the veloci- 
ties at which frictional data have currently been measured. In 
order to know what frictional forces are developed in practice, it 
is obvious that tests must be made in the 50 to 200-fps velocity 
range. We can expect the coefficient of friction to be higher at 
the high velocities, but just how much higher we don’t know. 

The design of a tensiometer which will measure accurately the 
magnitude of tensions developed at high speeds of 5000 stitches 
per min is a very difficult problem. The probe must be so de- 
signed that there is a minimum of distortion of the sewing thread. 
In addition, the resonant-frequency characteristics should be 
sufficient to measure the short tension pulses which seem to occur. 
It is possible that a transducer of the piezoelectric type with a 
resonant frequency close to 10 ke could perform satisfactorily. 
This type of gage has a high output and might prove to be quite 
satisfactory in this application. The transducers which have 
been used to date have all been bonded, wire-resistance strain 
gages attached to a beam of some sort. These gages in the opinion 
of the author are inadequate to study sewing tensions quantita- 
tively. 

It is recommended that a study of sewing tensions be made 
with a tensiometer replacing one of the guides between the needle 
and the take-up mechanism. A minimum amount of distortion 
would then occur. This tension measurement should be made in 
conjunction with high-speed motion pictures. It can be seen that 
a complete study of sewing will thus involve complex and expen- 
sive equipment. It is recommended that this type of research be 
conducted in industrial laboratories where more funds are availa- 
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ble. However, industrial research is not published as quickly as 
academic research, and it might perhaps be better that this equip- 
ment be made available to academic institutions so that these 
could be written on data obtained with suitable equipment 
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= influence on  space-temperature-time history 
and heat-exchanger performance of a transient “‘hot-goods’’ load 
suddenly introduced into the refrigerated zone of an insulated 
structure is of particular interest in cold-storage warehouse 
operations. Investigation of this topic was undertaken initially 
via rheoelectric simulation, on a structure with walls of negligible 
heat capacity, described in a paper by the author at the Ninth 
international Congress of Refrigeration at Paris, France, in 1955 
{1),! and earried along further in 1957 [2]. Both papers inter- 
preted the contribution of the hot body to heat load via a 
concept of “transient-body”’ resistance, timewise variable at the 
exposed surface of the heat-flow body. 

Transient effects with cold-storage loads and finite heat- 
capacity wall under constant outdoor-temperature conditions, 
using experimental prototype simulation methods, were further 
explored ina recent ASME paper [3]. The analysis is now further 
extended, as a continuation paper, by treating transient effects of 
an insulated-wall structure having actual heat capacity, and with 
cyclically varying outdoor-air temperature as well as of solar- 
radiation effect on the exposed surfaces. The analysis is again 
coupled to rheoelectric-simulation analogy (RESA) technique [3]. 

In this study, the different rates of energy flow are followed, 
timewise, using specialized electrical-potential transfer processes. 
This is done in conjunction with the resistance-capacitance repre- 

1 Numbers in brackets designate References at end of paper. 

Contributed by the Process Industries Division and presented at 
the Annual Meeting, Atlantic City, N. J., November 29-December 4, 
1959, of THe AMERICAN Society OF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, Sep- 
tember 8, 1959. Paper No. 59—A-181. 


Effect of Transient Loads on the Refrigerated- 
Space Temperature of a Cooling Complex 


sentation for the heat-flow components of both hot-load and of 
wall-structure, previously treated [3] and hence not necessarily 
to be greatly detailed here beyond the needs of continuity. The 
net energy-flow is ultimately sensed in the heat-exchanger simula- 
tion circuit, as set up for the specific exchanger with given rates 
of brine-flow, circulating space-air, and initial brine-supply tem- 
perature. The space-air temperature is the equilibrium result of 
structure heat-inflow versus exchanger-outflow. In the illustra- 
tive examples, space, wall, and goods temperatures, heat-ex- 
changer energy-flow rates, and exchanger-fluid temperatures 
are predicted as a function of elapsed time, for several cases 
covered in the program of experiments. 

The system studied is shown typically in Fig. 1, along with a 
generalized circuit to represent the heat-flow system. The heat- 
load of the refrigerated space is absorbed in the movement of the 
space-air, cooled through fan circulation in a countercurrent type 
of heat-exchanger supplied with constant-temperature cold brine 
pumped from the separate brine-cooler evaporator. 

Studies via rheo-electric analogy, as extensions of the Resist- 
ance Concept of energy-flow, have covered widely-varying 
applications: (a) In steady-state heat conduction both the 
geometrical-sheet and network methods have been used by the 
author [4,5]; (>) heat-exchanger analysis [6] as herein involved, 
and, (c) further, pipe-network fluid-flow distribution and energy- 
transport {7]—all, differing facets of the RESA technique. 

As continually noted in studies such as these, many simplifying 
assumptions are found to be expedient, and among them are: 
(a) The heat-exchanger and the associated fluids are assumed 
to have negligible heat-capacity; (b) values for the thermo- 
physical properties have been assumed, and taken independent of 
temperature; (c) values of air heat-transfer coefficients, the time- 
varying rate of solar-energy received at the external structure- 


Nomenclature 


k/(pe) = thermal diffusivity, s gross thermal resistance, deg 


surrounding air temperature, deg 


ft/hr 

surface area, sq ft 

specific heat, Btu/(b-deg F ) 

thermal capacity per unit of 
frontal area = cpy 

electrical voltage for analog cir- 
cuit 

dimensionless temperature-dif- 
ference ratio 

heat-transfer coefficient, Btu/(hr- 
sq ft-deg F) 

electrical current for analog cir- 
cuit 

thermal conductivity of solid ma- 
terial, Btu/(hr-sq ft-deg F/ft) 

heat-flow rate, Btu/hr 

thermal resistance, referred to 1 
sq ft, deg F/(Btu/(hr-sq ft)) 
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F/(Btu/hr) 

thermal resistance of solid path 
y, per sq ft frontal area = y/k 

thermal resistance due to surface- 
coefficient h, per sq ft = 1/h 

heat-exchange resistance of air 
stream = 1/(W,¢,) 

heat-exchange resistance of brine 
stream = 1/(W,¢,) 

exchanger heat-transfer resistance 
= 1/(UA) 

electrical resistance, ohms, pro- 
portional to RP,’ 

electrical resistance, ohms, pro- 
portional to PR,’ 

over-all exchanger thermal re- 
sistance, deg F/(Btu/hr) 

electrical resistance, ohms, for 
heat-exchanger, proportional to 
R 


tee 


F 

external atmospheric tempera- 
ture, deg F 

original uniform solid tempera- 
ture, deg F 

surface temperature, deg F 

temperature in solid, at any point 
z distance out from center line, 
deg F 

temperature difference, deg F 

weight flow rate for air, lb/hr 

weight flow rate for brine, lb/hr 

exchanger over-all coefficient of 
heat transfer 

half slab thickness, ft 

density, Ib/cu ft 

elapsed time, hr 

mathematical function 


Subscripts to denote special applications 
are used where required. 
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surface, and, on occasion, the outdoor-air temperature cyclic 
variation, have likewise been assumed. 

(In passing, and in deference to international technology, it 
might be noted that one could wish that specific unit-systems 
need not be employed to portray values of physical properties 
directly, i.e. some metric or English system, that some universal 
system could be used. Work on such a universal dimensionless 
system is under way by the author, not yet sufficiently advanced 
to utilize in this presentation, but will be presented and used at 
an early future date.) 


In general, in this paper, rheoelectrical analysis is carried to the 
further development of the transient refrigerating-load problem 
originally studied [1]. To reiterate, in this work pilot resistance- 
capacitance units are used to monitor the changing heat-flows 
growing out of the time-varying temperature differences between 
the transient-surfaces and the cooled-space. 


Heat-Load Relationships 


Transient Elements With Heat Capacity. 
“hot’’ 


These consist of: (a) The 
load, presumed to be a slab-form body initially uniformly 
at outdoor temperatures and then suddenly exposed to air-cooling 
effect on its surfaces, and (b) the walls of finite heat-capacity, sub- 
jected to inward heat-flow effects on the outside surface, and out- 
ward heat flows to the space-air on the inside surface. The 
actual heat flow to the,cooling-air, in both the “‘slab’’ load and the 
wall case, may ultimately be represented in terms of the relevant 
temperature difference of surface versus the surrounding air. 

Considering, illustratively, the case of unidirectional heat flow 
for a solid-slab section of'y half-thickness, of given thermal prop- 
erties and unit frontal area (1 sq ft), suddenly exposed to cooling 
action with air-coefficient h, (R,’', corresponding thermal resist- 
ance), #,’ as the thermal resistance of the solid material y and 
with a as the thermal diffusivity: Then a useful dimensionless- 
ratio / may be set up [3]. Specifically, for sustained tempera- 
ture-difference Ato. between initial uniform solid temperature 
to and the surrounding air ¢,, and for a point z distance out from 
the center line toward the surface, with position ratio z/y, the 
temperature difference F, may be specified and ultimately 
equivalently shown in terms of the Biot number Cg; and the 
Fourier number and the position ratio 


—, = 
R, y 


1960 


160 7. MAY 


| 
OveRALL At 
EXTERNAL - AIR 


vs 
INCOMING - BRINE 
At t.b 
| REFRIGERANT 
REFRIGERATION EVAPORATING REFRIGERANT 


TEMPERATURE 


WARM-AIR 

TEMPERATURE © 

(COOLED - 
SPACE.) 

_ WARM - BRINE 

~ TEMPERATURE 


2 COLD-AIR 
TEMPERATURE 


TEMPERATURE 


HEAT-FLOW RATE 


General arrangement of refrigerated complex and thermal circuit 


REFRIGERANT t 
TEMPERATURE 


It follows, therefore, that actual temperatures may now be speci- 
fied: 


+ F, Ato. 


(2) 

In the rheoelectric simulation technique, the operation is car- 
ried out by analogy; throughout, electrical values in proper scale 
relationship to thermal values are used and are ultimately 
translated back into thermal values as end result. As noted pre- 
viously and detailed in the earlier report [3], the thermal struc- 
tures (hot body of slab form and wall of finite heat capacity) 
are dealt with in terms of a resistance-capacitance electrical proto- 
type, based on unit frontal area, with the full heat flow then dealt 
with in terms of the total exposed area A,, the surface heat- 
transfer coefficient h,, and the surface-to-fluid temperature dif- 
ference At, .. 

The gross heat flow g,; between a transient body with varying 


surface-temperature ¢, and surrounding air at t, may, as indi- 


cated, be stated in terms of the temperature difference At, , with 
an applicable resistance, R, = 1/(A,h,): 
A tl, a , 
¢ = — 
R, (3) 


Wall of Negligible Heat Capacity. In contrast to transiently- 
operating walls of finite heat capacity, walls of negligible heat 
capacity are more readily handled. The heat flow is given pri- 
marily via an equivalent wall-resistance R,, which may also 
provide for ventilation and infiltration-air cooling, and the over- 
all temperature difference At, between outside and inside space- 


air, or alternately as the product of this At, and the gross heat- 
transfer rate UA: 


(4) 

Solar-Radiation Effect. The solar-radiation effect, incidentally, 
will be established via timewise varying unit-area heat input on 
the external-surface side of the wall. Its effect will ultimately be 
apparent in the temperature behavior of the external wall surface. 


Heat-Exchanger Relationships 


The heat exchanger, for air-to-brine heat flow, may be con- 
sidered on a countercurrent basis. As developed in earlier papers 


[6, 3], its performance may be explored and predicted through an 
extension of the Resistance Concept. The exchanger heat flow q, 
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is related to the over-all exchanger temperature-difference At,, 
between the entering air ¢, and the initial brine ¢,, through 
the over-all resistance F,,,: 


q = (5) 


R, - R,e(Re— Re)/Rua 
1 — e(Ra—Re)/RvA (6) 


= 


R,and R, represent the heat-exrchange resistance of the air and 
brine streams, and Ry,, the over-all heat-transfer resistance within 
the exchanger. 

As will be apparent in the specific examples to be cited, pro- 
portional resistance along the total resistance-path R,, will 
represent a measure of the temperature change for the different 
fluids in the exchanger. 


Electrical Simulation of the Thermal Processes 


Solid bodies, including walls, undergoing transient change, 
comprise both thermal-resistance and thermal-storage, i.e., heat- 
capacity, effects. Their internal behavior, therefore, may be 
simulated electrically, through electrical resistance and electrical 
capacity, multiple-slicewise, on a lumped network basis, and 
incorporated in the general energy-flow circuit, through a “‘pilot’’ 
prototype arrangement. 

Fig. 2 represents the general rheoelectric simulation circuit, 
typically for the study of the subsequently cited illustrative cases. 
It is made electrically in proper scale relationship to thermal, for 
example, electrical resistance versus thermal resistance. In addi- 
tion to the resistive components for the structure and for the heat- 
exchanger, load-transfer slide-wires are also included for the pilot 
monitoring operations. Via the pilot units, ‘surface-tempera- 
ture’’ voltages result from the operation of the rC prototype net- 
works, their values being reproduced at the load slide-wires via 
high-impedance elements. In this way it is possible conveniently 
to track the real load change for the main circuit as a function of 
time. Thus the flow of heat off the exposed surfaces of the trans- 
ient loads to the space-air may be simulated, using load-resistance 
values as indicated via Equation (3). The time effect contributed 
through the electrical condensers provides, via appropriate time 
scale, the thermal time basis for the different simulation experi- 
ments. 


Illustrative-Studies Program and Results 


} For the refrigerated space and heat-exchanger complex repre- 


sented by Fig. 1, timewise performance prediction is desired 
for certain imposed conditions. In the illustrative cases herein 
reported, two thermally-different structures are involved, each 
with the same 2400 sq ft of finite-heat-capacity wall, but with 
4800 sq ft of negligible-heat-capacity walls of dissimilar heat-flow 
characteristics, representing Series I and Series II with separate 
schedules of operating conditions. The same heat-exchanger and 
auxiliary services are provided for each series, with brine at the 
fixed supply temperature of 12 F. 
Thermal data are as follows: 


1 Hot load: A tightly packed vertical stack of goods, 2-ft- 
thick, of 1250 sq-ft exposed area; thermal conductivity, k = 
0.35 Btu/(hr-sq ft-deg F/ft); thermal diffusivity, a = 0.017 sq 
ft/hr; air surface-coefficient, h, = 2.0 Btu/(hr-sq ft-deg F); 
load uniformly originally hot at 100 F. 

2 Structure-envelope: Series J—2400 sq ft of concrete-cork 
wall, with 8-in. (0.667-ft) concrete, k = 0.50, a = 0.0139; and 
with 6-in. (0.50-ft) cork, k = 0.025, a = 0.005; —4800 sq ft of 
negligible heat-capacity wall, gross heat-flow rate of 650 Btu/(hr- 
deg F); Series ]I—same concrete-cork wall as in Series I, but with 
4800 sq ft of negligible-heat-capacity wall having 213 Btu/(hr-deg 
F) gross over-all heat-flow rate. Air heat-transfer surface co- 
efficient, h = 1.65 Btu/(hr-sq ft-deg F), as a general value for 
both exposed sides of the walls. 

3 Heat-exchanger: Brine flow-rate, W, = 8000 lb/hr; specific 
heat, c, = 0.70 Btu/(Ib-deg F); circulating air-flow, W, = 
38,500 lb/hr; c, = 0.241; gross heat-transfer rate UA = 11,850 
Btu/(hr-deg F) 

4 Solar radiation: The solar load, presumed to repeat daily, 
and received at the outside concrete surface, with appropriate 
allowance for absorptivity-coefficient, location, season (late 
summer), has been illustratively assumed according to the follow- 
ing symmetrical schedule, in Btu/(hr-sq ft): 


Time, hr: O04 5 6 7 8 9 10 11 12 
Heat flow: 0O 6 27 55 88 110 121 128 132 


Time, hr: 13 14 15 16 17 18 19 20-24 
Heat flow: 128 121 110 88 55 27 6 0 


5 Sinusoidal diurnal outdoor-air-temperature variation has 
been assumed, as evident directly in the appropriate figures. 
It is to be noted that the hot-load, under the programs of 
operation, will be subjected to cooling heat transfer under time- 
wise-variable space-temperature. Accordingly, for reference 
purposes, the general characteristics of temperature-variation 
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under sustained cooling conditions, indicated via the dimension- 

less ratio F, will be of interest, and similarly for the concrete-cork 
finite-heat-capacity wall. These characteristics have been deter- 
mined separately through the electric-analogy resistance-capaci- 
tance study of the prototype models, and have been carried over 
from the earlier study [3]. 

Fig. 3 shows the ratio F, for both surface and core of the hot 
load, versus elapsed time, under constant-level cooling conditions. 
Fig. 4 is a similar study for the finite-capacity wall, and shows 
the value of F as a function of elapsed time; this is for different 
positions within the wall including the concrete-cork interface, 
while ultimately achieving steady-state heat flow between the 
high and the low-temperature air-realms on the two exposed faces 
of the initially uniformly-warm wall. Fig. 5, with dimensionless F’ 
plotted against the traversed wall-resistance, highlights the idea 
of steady-state as compared with the 10-hr intermediate transient 
value derived from Fig. 4. 

Figs. 6 to 9 represent Series I results, and include the effects of 
external solar radiation to the 2400-sq-ft concrete-cork wall sur- 
face on the performance of the structure-complex. Here the out- 
side air-temperature is presumed to stay constant at LOO F. It is 
interesting to trace the different internal temperatures and the 
heat flows 

Figs. 6 and 7, for temperatures and heat-flow rates, respec- 
tively, show the effect of the solar load on the empty structure 
and heat-exchanger complex as a function of time, with the in- 
fluence of the diurnal heat-flow variation very evident throughout 
the various parts of the complex. 

Figs. 8 and 9 concern the complex with the hot load suddenly 
introduced into the cooled space at 9 a.m. on the first day of this 

portion of the study. It is interesting to note that the outside 
concrete surface temperature just begins to decline after the 
third-day peak, with the effect of the hot load within the space 
quite marked, as inherently superimposed on the solar effect. 
Temperatures for the load-material are also shown, and bear com- 


parison with the reference values indicated by Fig. 3. 
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Fig. 3 Dimensionless ratio F for hot load, surface and core, versus 
elapsed time 
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Fig. 7 Effect of diurnal solar heat flow on system energy-flow-rates; 
elapsed time (Series 1) 
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Fig. 9 Effect of combined solar-energy and hot load on system heat flows (Series 1) 


Figs. 10 to 12 represent Series IT results. They cover three peri- 
ods of four days each, showing temperature values as a function 
of time, well worthy of study and particularly illustrative of the 
value of this type of rheoelectric-simulation analog analysis. 
This is especially to be noted in terms of the very involved situa- 
tion of solar-heat variation, external air-temperature variation, 
structure time lags, and transient-load heat flow, all coupled to 
the heat-exchanger performance of the complex. 

Fig. 10 represents Period 1 of the Series II study. It depicts a 
prior steady-state condition for the empty structure, with a con- 
stant outside temperature of 80 F. Two days of the prior steady- 
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state are shown, then followed by two days with cyclically-vary- 
ing outdoor temperature (from 80 to 100 F). It is interesting to 
note the different structure temperatures, especially that of the 
outside surface of the concrete wall, and to follow them through 
their succeeding histories. In Fig. i1, Period 2 next follows 
directly, with Day 1 continuing with merely the outside-air varia- 
tion. Then, with the beginning of Day 2 of this Period 2, the 
solar-radiation cycle is added, in accordance with the given 
schedule. It is of course evident, among other items, that the 
outer, concrete, surface temperature cyclically exceeds the vary- 
ing outdoor temperature, under which circumstances some of the 
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solar energy is given up to the outdoor air, with the remainder 
entering the wall and causing the various temperature disturb- 
ances so evident in the figure. 

Finally, Fig. 12 for Period 3 is to be noted, as a further con- 
tinuation of Period 2 shown in Fig. 11. At about 1 p.m. of Day 1, 
the hot-load of the specified design characteristics and at 100 F 
has been delivered to the warehouse and is promptly charged 
into the empty refrigerated space; then follow the myriad of dis- 
turbances that are so evident and which are continued to be de- 
picted until the close of Day 4 of this Period 3 study! 

The results for Series I and II represented by Figs. 6 to 12 are 
particularly illustrative of the value of this type of rheoelectric- 
simulation analog analysis, a study which, by orthodox computa- 
tional methods, would indeed prove quite cumbersome. To re- 
iterate, this is especially noteworthy in terms of the very in- 
volved situation of solar-heat and outdoor-air temperature varia- 
tions, the different structure time lags, and the transient-load 
heat flow, all coupled to the heat-exchanger performance of the 
total complex. 
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stretching zone. 


flow lines. 


Uniform Two-Way Orientation 
of Plastic Films 


In order for a plastic film to be stretched at equal rates in the longitudinal and the trans- 
verse direction a highly specific geometrical flow pattern must be established in the 
This flow pattern involves diverging radial flow lines; the velocity at 
any point is directly proportional to the distance from the convergences point of the 
These geometrical facts lead to certain design requirements for a machine 
to be used in the production of plastic film with uniform two-way orientation. 


In order 


for new concepts to be expressed in terms of practical working designs for machinery, 
mechanical development is required. The reduction to practice of the foregoing concept 


PART 1 


N THE production of thin plastic films, an important 
factor is the nature and extent of molecular orientation. Orienta- 
tion is achieved by hot-stretching or cold-stretching of the film at 
some stage of manufacture. The physical properties of the final 
plastic film are very sensitive to the molecular orientation. For 
example, unoriented polystyrene film is brittle and fragile, whereas 
the same material when properly oriented is strong and tough. 

If the orientation is in one direction only, the plastic film will 
have different properties in different directions. For example, 
a polystyrene film which is oriented only along the machine 
direction may be very strong and tough in this direction, but will 
fibrillate easily when subjected to transverse tensile stresses. For 
this reason, it is often desirable to impart equal molecular orienta- 
tion in all directions within the plane of the film—which implies 
a film-stretching step in which the film is subjected to equal 
longitudinal and transverse stretching. Part 1 of this paper 
analyzes some of the geometrical and rheological factors involved 
in such a uniform two-way stretching operation. 


The Kinematics of Uniform Two-Way Stretching 


If all regions of a piece of film are being stretched at equal rates 
in all directions of the plane, the velocity distribution must con- 
form to a highly specifie pattern; all particles move in diverging 
straight lines, with velocities proportional to the distance from 
the convergence point of these straight lines. 

This pattern of motion is a two-dimensional analog of the ex- 
panding universe, in which galaxies (in general) are moving away 
from each other at speeds proportional to the separation dis- 
tances. Thus if an observer on any galaxy considers himself to 
be at rest, he observes all other galaxies to be moving away from 
him along diverging radial straight lines, at velocities propor- 
tional to distance from the origin. As far as the relative motion of 
the galaxies is concerned, any point in the universe can be taken 
as the center of co-ordinates; i.e., as the convergence point of the 
lines of motion. 

Likewise, we can imagine that a piece of plastic film which is 
undergoing uniform two-way stretching is populated by a race of 
scientific ants, each of whom observes the motion of the others 
away from him. Each ant can consider himself to be at rest, and 
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is shown by the accompanying illustrations. 


thus at the convergence point of the path of motion of the other 
ants. However, in this case, the earth and the film-production 
equipment provide a unique frame of reference. In terms of this 
reference frame, there will be some particular point from which 
all the velocity vectors diverge. Using this point as the center 
of a polar co-ordinate system (r, @), the velocity pattern corre- 
sponding to equal two-way stretching can be written very simply: 


ve = 0 v, = kr or, v = ke (1) 


Here v is the vectorial velocity of any particle of the film (with 
angular and radial components equal to vg and »,), and r is the 
vectorial location of the particle (with co-ordinates 9 andr). The 
rate of tensile straining of the film is equal (in all directions in the 
plane, and at all locations) to k. 

The foregoing type of velocity pattern is necessary in order for 
a flat film to be stretched at equal rates in all directions. If we 
are dealing with a “‘batch’’ operation on a small, bounded piece of 
film, k can be a function of time. Also, we can superimpose any 
arbitrary motion of the center of co-ordinates, O, as long as the 
motion of the film relative to O follows equation (1). 

However, we are most interested in the steady-state operation 
of a continuous film-stretching device, in which an infinite length 
of film passes through a stretching zone within which it is sub- 
jected to equal rates of stretching in the longitudinal and the 
transverse directions. In such a case, k is a constant (rather than 
a function of time), and the point O is stationary (relative to the 
earth). 

Such a stretching zone can be either of the form shown in Fig. 
1(a) or of that shown in Fig. 1(b). The nature of the continuous 
supply of film to the inner, or starting boundary of the stretching 


\a) 


Fig. 1 Continuous film-stretching unit. 


(a) 360-deg stretching zone. 
(b) Limited sector stretching zone. 
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zone, and the methods of continuously removing the film from 
the outer, or finishing boundary of the zone will be discussed later 
in the paper. For the moment let us confine our attention to the 
stretching zone itself. 

Within a 360-deg stretching zone such as that shows in Fig. 
l(a), the film will be under tension. Let o, and ¢, represent the 
tensile stress in the tangential and radial directions, respectively. 
These stresses will be functions of r, but independent of 6. o, = 
o,(r); o, = 0,(r). Furthermore, as the film moves through the 
zone and stretches in two directions, its thickness, b, will decrease; 
b is a function of r, but independent of @. Finally, the velocity 
pattern which develops in the unsupported interior of the stretch- 
ing zone (ag a result of the stresses 7, and @,) will depend upon the 
rheological properties of the film material. At this point the 
problem becomes one of mechanics, rather than simply kine- 
matics. In the next section we shall analyze the mechanics of 
film-stretching in such a zone, on the assumption that the ma- 


terial behaves as a high-viscosity Newtonian fluid (with sufficient 


tensile strength to support the stresses o, and oa, without rup- 
ture). 


Mechanics of Uniform Two-Way Stretching (Film Exhibiting 
Newtonian Flow) 


Statement of the Mathematical Problem. We are concerned with 
a thin sheet of a highly viscous (Newtonian) material. The 
sheet is flat and of the’ form indicated in Fig. 1. The sheet is 
continuously flowing outward—being fed continuously at the 
inner boundary and removed continuously at the outer boundary. 
The sheet is everywhere in a state of two-dimensional tension. 
The radial tension will be represented by ¢,; the tangential 
tension by o,. Thickness of the sheet will be represented by b. 
The fluid velocity v is everywhere radial in direction. The vis- 
cosity will be represented by 7. 

Only the steady-state condition of the system will be treated. 

The following quantities are, in general, functions of r (but in- 
dependent of @): n, b, o,, 

Ordinarily, we shall start with a certain 7n(r)—assumed 
known—and a set of boundary conditions (v at inner and outer 
boundary), and proceed to compute o(r), b(r), o,(r), o,(r). 

The Differential Equation for the System. The radial velocity dis- 
tribution v,(r) is governed by the following differential equation 
(7 is a known function of r): 


(= 1 dv 1 
dr? ” \ dr \ dr v \@adr 
2r \dr 2nr \ dr a 


If »(r) has been determined by solving (2) subject to definite 


boundary conditions, b(r), ¢,(r), and o;,(r) can be computed as 
follows: 


Q 


2rrv(r) 


b(r) = (3) 


where 


Q = total volumetric feed rate across the inner boundary 


dv v 
(4) 
dr 
dv v 
o, = 2n — |} + 491 — (5) 
dr r 


The derivation of equations (2), (3), (4), and (5) follows: 
Balance of Forces. The condition of force balance on each ele- 
ment of the viscous membrane is 
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d 
— (rbo,) — bo, = 0 (2a) 
dr 
This equation can be expanded as follows: 


b 
+ ro, (4) + bo, — bo, = 0 (2b) 
dr dr 


Steady-State Material Balance. In a steady state of flow, the total 
outward flux of material at each value of r is equal to Q (the total 
volumetric rate of extrusion): 


Q = 2rrb(r)v(r) (3a) 
= (3b) 


Relationship Between Velocity Function v(r) and Rate of Strain Tensor 
é,): 


Relationship Between Velocity Function v(r) and Stress Tensor 


o,=0 (5a) 


Inserting Equations (3b) and (5a) Into (2b) Leads to Basic Differential 
Equction: 


ar, n \dr dr v \dr 


2r \dr 2nr \ dr Qr2 


Solutions of Equation (2). When a particular (7) is inserted into 
equation (2), the resulting differential equation is (in general) too 
complicated to solve analytically. It is quite amenable to ma- 
chine calculation. 

Certain specific conditions lead to simple analytical solutions of 
(2). For example, if n(r) = Cr®, a solution can be found of the 
form v(r) = Ar™. In order to utilize this solution, not only must 
n(r) be of the correct form, but also the boundary conditions must 
be consistent with the solution [v,/v, = (7,/r,)™]. This is thus a 
highly special case. Up to now, no »(r) has been found which 
permits general solution of (2). 

In spite of the highly special nature of the foregoing case, it 
will be discussed briefly. 

If n(r) = Cr* equation (2) becomes 


d’v (2n — 1) f dv 1 f{ dv \? v(n — 3) 

(4) + 2r () v (*) 2r? 

Assuming a solution of the form o(r) = Ar™ and inserting in (6), 


we obtain 


= 7 
2n — 3 (7) 


m 


This relation between n and m is plotted in Fig. 3. Note, how- 
ever, that m must be greater than ( —0.5) or the sheet will buckle 
radially. The left branch of Fig. 3 is eliminated for this physical 
reason. 

The case where 7 = Cr? is of particular interest, since it leads 
tov = Kr. In this case, dv/dr = v/r; hence the radial and tan- 
gential tensions are identical; i-e., we have a condition of equal 
two-way stretching rates at all points. 
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Applicability to Non-Newtonian Fluids and Viscoelastic Films. 
Under the conditions encountered in film production, polyethyl- 
ene is far from being a Newtonian fluid. Its rate of flow is not 
simply proportional to the stress; and in the latter part of the 
process, deformation is largely elastic, in order to achieve molecu- 
lar orientation. 

No attempt has been made to extend the foregoing mathemati- 
cal analysis to non-Newtonian and viscoelastic fluids, although 
this could be done if desired. The reason is that such an extended 
treatment would leave unaltered some of the geometrical con- 
siderations which arose in the simple treatment. Specifically: 
Equations (2a) and (2b) apply in the case of non-Newtonian and 
elastic fluids. Also, equations (3a) and (3b) continue to hold. 
Finally, equations (4a) remain valid in the extended treatment. 
Only equations (5a) must be changed from their original form. 
This of course changes the governing differential equation from 
(2) to some more complicated form. However, the flow pattern 
v(r) = Kr will still represent equal rates of stretching in the two 
directions (see equations 4a). 

The chief change introduced by non-Newtonian and visco- 
elastic behavior is that the femperature distribution needed to 
achieve v = Kr is no longer given by equation (9). 

Thus although equation (2) applies rigorously only to New- 
tonian fluids, certain conclusions derived therefrom can be applied 
with confidence to real plastic films. 


Inner Boundary 


Outer Boundary 
Fig. 2 Radial and tangential stress components Idealized “Pancake” Form of Stretcher 
The foregoing theoretical analysis of the basic kinematics and 
mechanics of uniform two-way stretching leads to the following 


idealized concept of a mechanical device to actually carry out 
such an operation. We have called this idealized machine a 


Fig. 3 Relationship between mand n. n(r) = Cr®; vir) = Ar™; m = 
(3 — n)/(2n — 3). 


r 


If the viscosity of the material varies with temperature ac-  Fig.4 Type of temperature-distribution required in order to achieve equal 
cording to the Arrhenius law, rates of radial and tangential stretching (Newtonian fluid) 


= (8) 


the temperature distribution with respect to r must be 


B 
InC/A+2Inr (9) 


in order to obtain the viscosity distribution Cr?. This tempera- 
ture distribution is of the general form shown in Fig. 4. 

Another case of interest is where 7 = Cr’; here m = 0 which 
means that (7) = K. In this case, 7, = 1/20, The orientation 
in the transverse direction will be much greater than the orienta- 
tion in the longitudinal direction. In order to achieve uniform 
two-way orientation in the film such a stretching step should be 
followed or preceded by a step in which additional longitudinal 
orientation is introduced. 


Fig. 5 Schematic diag of ke” stretcher 


Journal of Engineering for Industry may 1960 / 169 


lg 
| 
‘ 
. | 
! 
-4 
-5 
-6 | 
| 
| | 
| 
| 
= 
= If \ Se - 


pancake stretcher. Its form is indicated in Fig. 5. The polymer 
is extruded as a molten tube of radius Rf, and is immediately 
folded! over and stretched in the horizontal plane, and taken over 
a toroidal roll of radius Re. 

The diverging flow lines provide the transverse stretch; the dif- 
ference in take-up velocity and extrusion velocity provides the 
longitudinal (radial) stretch. 


PART 2 
Mechanical Development 


The statement of a fundamental principle sometimes necessarily 
includes certain assumptions of ideal physical situations. Two 
problems of this nature were apparent in considering the reduc- 
tion to practice of a film-stretching device based on the principle 
as outlined in the first part of this paper. The first problem was 
to get a toroidal roll, driven or otherwise, that would function 
properly, and would be mechanically sound. 

The second problem was to pull a cylinder of thin film uni- 
formiy, and to separate it into flat unmarked sheets. To sur- 
mount. these practical difficulties an approximation of this 
principle was considered to be possible. The second part of this 
paper deals with the reduction to practice of this approximation, 
and the subsequent steps taken to perfect the process. The first 
step toward using the foregoing principle was the concept of 
dividing the pancake into eight separate sheets, thus approximat- 
ing the circular boundary with an octagon, Fig. 6. 


Fig. 6 Eight-sided film stretcher 


In this concept the extrudate issues from a radial flat or pan- 
eake die, is drawn radially outward to point B where suitable 
cutters separate the web into sheets. 

The stretching to be carried on can be described as approximat- 
ing the following: 


1 The lateral stretching ratio is approximately equal to 


Diam at Bor D, R 
=f 
Diam at A or D, 


2 The longitudinal stretching ratio is approximately equal to 

? Actually, this ‘folding’? would be carried out inside the extruder 
die. The die would discharge the polymer in a horizontal peripheral 
stream. 
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Vel at B or V, 
Vel at A or V, 


=R 


If R, is fixed by the physical dimensions of the machine, R, can 
be made to equal #, by regulating the velocity at B with respect 
to A. 

At a fixed rate of extrusion, V, can be increased or decreased by 
lowering or raising the die opening, and V, ean be increased or 
decreased by changing the speed of drawdown. Thus for the 
situation V,/V, = D,/D, this concept can approximately equal 
two-way stretch. 


Fig. 7 Initial trial run of eight-sided film machine, Model 1 


Fig. 7 shows the initial trial run of such a device. Poly- 
ethylene was extruded from a pancake die and started by hand 
through the draw rolls. 

Those experienced in process development will recognize that 
this device was not yet a process for it had many shortcomings at 
this stage. Among these were: 


(a) The edges of the sheets were rough and baggy since the 
cutters were located too far ahead of the draw rolls. 

(b) The film heated up the rubber draw rolls to the extent that 
the film became embossed. 

(c) There were no provisions for gage control. 

(d) Many die lines were present. 


Fig. 8 shows the initial steps taken to correct these short- 
comings. The cutters were moved to give a better edge trim. 
The rolls were given a good finished surface to eliminate the em- 
bossing problem, and arranged alternately up and down to sup- 
port the web at the point of separation. After several trials the 
decision was made to scrap the original machine and build a better 
one, based on our findings. Fig. 9 shows the improved machine 
wherein the rubber draw rolls have been moved up to drive the 
finished metal rolls giving a more positive drawing mechanism. 
An adjustable radial extrusion die was also designed to improve 
gage control, and improved to eliminate die lines. 

Runs were then made to attempt gage control and to wind a 
good roll of film. 

The fact that the circular flat pancake die can be made readily 
adjustable was thought to be one of the most promising features 
of this process. It was found that good gage control was possible 
by following this procedure: 


(a) The eight sheets of film were gaged with micrometers 
and the low and high spots located. 
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Table 1 Elongations of experimental films 
materials 


Sample Percent elongation to break 


Longitudinal | Transverse | Diagonal 


This invention: 
Single sheet 489 492 
Eight sheets (avg.) 477 486 


Commercial products: 


Monufocturer A------------- 
Manufacturer B--------------- 
Monufaocturer 
Manufacturer 
Manufocturer 
Manufacturer 


(b) Using crayon marks, the position on the face of the die was 
determined which corresponded with the imperfection in the 
sheet. 

(c) It was then necessary to adjust the proper screw to effect the 
correction. 


It was observed that the temperature of extrusion had a pro- 
nounced effect on the quality of the sheet in so far as flatness and 
clarity were concerned. 

As the extrusion temperature was lowered below 200 C, using 
polvethylene, the sheets became very flat and entirely free of 
creases or bags, while the opacity of the film increased slightly. 

At higher extrusion temperatures some sagging occurred which 
allowed air pockets to form where the sheet contacted the roll, 
forming creases in the sheet. 

Other observations indicated that when polyethylene remained 
in the extrusion system during shutdowns, excessive gels were 
formed. These were purged out of the system without too much 
difficulty. 
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Fig. 10 Eight-sided film machine, Model No. 4 


In addition to this, there was usually a formation of decomposed 
polymer on the die surfaces which caused die lines and influenced 
the gage. 

An extended shutdown usually required that the dies be re- 
moved, cleaned, and reset. 

The rate that was attained at this point was limited by the 
capacity of an available 3!/:-in. extruder. One hundred and fifty 
pounds an hour was reached with no difficulty. Arrangements 
were made to load a tank from the extruder and then by pumping 
from the tank much higher rates were realized. Rates in excess of 
400 lb per hr were achieved with no difficulty. The top rate re- 
mains unknown as of this writing. 

To prove that a uniform two-way stretch takes place required a 
testing program carried out concurrently with the development 
work. Table 1 shows the per cent elongation required to break 
the film. This is the only table of data included in this paper since 
it is sufficient to illustrate the point that the film properties are 
very uniform. A study of Table 1 shows that this film has far less 
variation in elongation than commercially available films made 
by conventional blowing and flat-sheet-extrusion techniques. 
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Fig. 9 Eight-sided film machine, Model No. 3 
Fig. 8 Eight-sided film machine, Model No. 2 
480 
480 
297 667 --- 
490 705 
25! 629 --- 


Fig. 10 shows the final embodiment of this device. The ac- 
cumulated knowledge and experience gained in the course of the 
development are here expressed in terms of a new machine. This 
machine has an improved drive system, powered cutters, a com- 
pletely new frame, uniformly cooled high quality draw rolls, and 
the final embodiment of an adjustable radial extrusion die. 
Fundamentally sound in principle, and mechanically sound in 
construction, this machine is capable of producing eight sheets of 
quality film continuously at high production rates. 

This paper points out that it is possible to start with a funda- 
mentally new approach to an old problem and end up with a 
totally new type of process. 

Also, it shows that, in order to do this properly, mechanical de- 
velopment work is required to develop suitable equipment that 
can be used to produce the desired results. 

In closing, it should be stated that this paper covers a great 
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deal of work over a long period of time. It hits the high spots, 
and does not, or could not cover this development in detail. We 
hope, however, that it does show how such a problem can be ap- 
proached and solved. 

Norte: The process shown in this paper is covered by U.S. Pat- 
ents: No. 2,769,200, No. 2,852,813, No. 2,779,053. 


Glossary of Terms 
= viscosity 
thickness 
radial velocity 
radial stress 
tangential stress 
radial rate of strain 
tangential rate of strain 
volumetric rate of extrusion 
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